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Abstract 

The stability of quantum systems to perturbations of the Hamiltonian 
is studied. This stabihty is quantified by the fidehty, an overlap of an 
ideal state obtained by the unperturbed evolution, and a perturbed state 
obtained by the perturbed evolution, both starting from the same initial 
state. Dependence of fidelity on the initial state as well as on the dynamical 
properties of the system is considered. In particular, systems having a 
chaotic or regular classical limit are analysed. The fidelity decay rate is 
given by an integral of the correlation function of the perturbation and is 
thus smaller the faster correlation function decays. Quantum systems with 
a chaotic classical limit can therefore be more stable than regular ones. If 
the perturbation can be written as a time derivative of another operator, 
meaning that the time averaged perturbation vanishes, fidelity freezes 
at a constant value and starts to decay only after a long time inversely 
proportional to the perturbation strength. In composite systems stability 
of entanglement to perturbations of the Hamiltonian is analysed in terms 
of purity. For regular systems purity decay is shown to be independent of 
Planck's constant for coherent initial states in the semiclassical limit. The 
accelerated decoherence of macroscopic superpositions is also explained. 
The theory of fidelity decay is applied to the stability of quantum compu- 
tation and an improved quantum Fourier transform algorithm is designed 
and shown to be more stable against random perturbations. 

Keywords: quantum stability, fidelity, purity, quantum chaos, de- 
coherence, entanglement, quantum computation, quantum Fourier 
transformation, the kicked top, Jaynes-Cummings model. 
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Chapter 1 

Introduction 



In science one tries to tell people, in such a way as to 
be understood by everyone, something that no one ever 
knew before. But in poetry, it's the exact opposite. 

— Paul Dirac 



Irreversibility of macroscopic behaviour has a ttracted attention ever since Boltzmann intro- 
duced his H-theoiera in a seminal paper in 1872 (jBoltzmannl . Il872l ) . together with the equation 



bearing his name. The //-theorem is a precursor of what is today known as the second law 
of thermodynamics, stating that entropy can not decrease in the course of time. The problem 
was how to reconcile this "arrow of time" with the underlying reversible microscopic laws. How 
come that macroscopic systems always seem to develop in one direction whether the underlying 
dynamics is symmetric in time? This so called reversibility paradox is usually attributed to Jose f 



Loschmidt. He mentioned it briefly at the end of a paper published in 1876 (jLoschmidtJ . 
discussing the thermal equilibrium of a gas subjected to a gravitational field, in an attempt to 
refute Maxwell's distribution of velocities for a gas at constant temperature. He also questioned 
Boltzmann's monotonic approach towards the equilibrium. Discussing that, if one would reverse 
all the velocities, one would go from equilibrium towards the initial non equilibrium state, he 
concludes with : "...Das beriihmte Problem, Geschehenes ungeschehen zu machen, hat damit 
zwar keine L osunq. . . ". Bolzmann was quick to answer Loschmidt's objections in a paper from 



1877 (Boltz mannl . [12221)) pointing out the crucial importance of the initial conditions and of 



the probabilistic interpretation of the second law. The fact that for macroscopic systems we 
always observe a definite "arrow of time" is a consequence of a vast majority of the initial con- 
ditions representing an equilibrium state. If we choose the initial state at random it will almost 
certainly evolve according to the second law of thermodynamics. Of course this probabilistic 
interpretation still does not solve the problem. Resolution lies in the initial condition, chosen 
to be a macrostate representing a very small part of the phase space, i.e. having a low en- 
tropy. Soj_2rovided we start with a low entropy state, the above probabilistic arguments can be 
used ( Lebowit j . Il999) to explain the second law. The question of the initial low entropy state 



of say the whole u niverse still remains. For a popular account of this "cosmological" subject 



see 



Penrose! (|l989l ^. 



The first written account of the reversibility paradox is actually not due to Loschmidt but 
due to William Thompson (later Lord Kelvin) , although it is possible that Loschmidt mentioned 
the paradox privately to Boltzmann before. Boltzmann and Loschmidt became good friends 
while working at the Institute of Physics in Vienna around 1867 (directed at the time by the 
Slovenian physicist Jozef Stefan). In 1869 Boltzmann moved to Graz, but returned to Vienna in 
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the per iod 1873-1876. Fo r a detailed biography of Boltzmann see ICercignanil (|l998h . In a vivid 
paper ( ThompsonL Il874l l from 1874 Thompson gave a very modern account of irreversibihty. 
When discussing the heat conduction and the equahzation of temperature he says: "...// we 
allowed this equalization to proceed for a certain time, and then reversed the motions of all 
the molecules, we would observe a disequalization. However, if the number of molecules is 
very large, as it is in a gas, any slight deviation from absolute precision in the reversal will 
greatly shorten the time during which disequalization occurs... Furthermore, if we take account 
of the fact that no physical system can he completely isolated from its surroundings hut is in 
principle interacting with all other molecules in the universe, and if we believe that the number 
of these latter molecules is infinite, then we may conclude that it is impossible for temperature- 
differences to arise spontaneously...". The interesting question is then, how short is this short 
disequalization time? The quantity in question is nothing but the fidehty or the Loschmidt 
echo as it is sometimes called. We evolve the system forward in time with the unperturbed 
evolution, then backward in time with the perturbed evolution, and look at the overlap with 
the initial state. It can be considered in classical mechanics as well as in quantum mechanics. 
The decay time of the fidelity will then be the time of disequalization in question. Despite 
its importance for thermodynamics the fidelity was not considered until some years ago, with 
motivation coming from quantum rather than classical theory. 

Quant um theory is arguably the g reatest achievement in 20th century physics. There are 
estimates (JTegmark fc Wheeled . 120011 ) that up to 30% of the gross national product of the US 
relies on quantum devices. Alone the semiconductor industry is of enormous importance. Still, 
all these quantum devices do not manipulate individual quanta but rather exploit macroscopic 
phenomena involving many particles. Experiments involving individual quantum systems be- 
came possible only in the '80 with the progress made in e.g. manipulation of cold atoms in 
traps, single electron devices, entangled photons etc. . This was so to say the e xperimental birth 



of what is now called quantum information theory ((Nielsen &: Chuand . 120011 ) . Quantum infor- 
mation theory married quantum mechanics with information theory and with computer science. 
It deals with means of processing and transmitting information, and by using quantum systems 
can achieve things not feasible in any classical way. For instance, one can teleport a quantum 
state, or perform secure communication over a public channel or do a quantum computation. 
By using quantum resources to do a computation one is able for instance to factorize a number 
in a polynomial time, which is presently not known to be possible by classical computer. Also, 
quantum com puters are very e fficient in simulating other quantum systems, answering the prob- 
lem posed bv iFevnmanI (1983). Namely, he asked whether it is possible to build a computing 
machine whose size will grow only linearly with the size of a quantum system simulated on it. 
With classical computers this is not possible since the size of the Hilbert space needed merely 
to describe the system grows exponentially with the number of particles. We do not know yet 
if it is possible to build a quantum computer that will achieve that goal, but on paper quantum 
computer will be able to do the trick. We say on paper because presently one is able to perform 
laboratory computations only on a few qubit (less than 10) quantum computer. The main ob- 
stacle are errors in the evolution, either due to unwanted coupling with the environment or due 
to internal errors. Therefore, the main goal is to build a stable quantum computer, resistant 
to such perturbations. The usual benchmark for stability is fidelity and therefore one ought 
to understand the behaviour of fidelity in different situations to know hot to maximise it. Yet 
again, the original push to study fidelity came neither from quantum information theory nor 
from thermodynamics, but from the field of quantum chaos. 

The exponential instability of classical systems is a well known and much studied subject. 
As the underlying laws of nature are quantum mechanical the obvious question arises how this 
"chaoticity" manifests itself in quantum systems whose classical limit is chaotic. The field 
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of quantum chaos mainly dealt with stationary properties of classically chaotic systems, like 
spectral and eigenvector statistics. Despite classical chaos being defined in a dynami cal wa' 



I 



It was easier to pinpoi nt tne signatures oi classical cnaos m stationary properties ijirLaaKei . 
199ll : IStockmannL Il999h . There were not as many studies of the dynamical aspects of quantum 



evoluti on in chaotic systern s , some examp l es bei ng studies of the reversibility of quantum evo- 
lution dShepelvanskv , 119831: ICasati et al. 1 119861 ') , the dynami cal localization ( Fishman e^ al. 



1991 



1982 ; iGrempel et al. I . ll^SJ) , energy spreading (|Cohenl . 1200(1 ) , wave-packet evolution ( Hellei , 



). Classical instability is usually defined as an exponential separation of two nearby tra- 
jectories in time. In quantum mechanics the state of a system is completely described by a 
wave function and so one could be tempted to look at the sensitivity of quantum mechanics to 
variations of the initial wave function. But quantum evolution is unitary and therefore preserves 
the dot product (i.e. the distance) between two states and so there is no exponential sensitivity 
with respect to the variation of the initial state. This, at first sight perplexing, conclusion has 
been reached because we compared two different things. Classical mechanics can also be stated 
in terms of a Liouiville propagation of phase space densities and this is also unitary. If we want 
to compare quantum and classical mechanics we have to compare them on the same footing. 
Quantum mechanics is a probabilistic theory, the wave function just gives the probabilities of 
measurement outcomes. Therefore we should also formulate classical mechanics in a proba- 
bil istic way as a propagation of probability densities in phase space. The idea, first proposed 
by IPeresI (|l984l ) , was to study not the sensitivity to the variation of the initial condition but 



with respect to the variation of the evolution. He compared two slightly different evolutions 
starting from the same initial state - the quantum fidelity. For classical systems fidelity gives 
the same exponential sensitivity to perturbations of the evolution as to perturbations of initial 
conditions; the two things are equivalent. The quantum fidelity though can behave in a very 
different way, as we will see in the present work. 

The fidelity lies at the crossroad of three very basic areas of physics: thermodynamics, 
quantum information theory and quantum chaos. The features of the fidelity turned out to be 
very interesting, as one would expect for such a crossroad. 

1.1 Historical Overview 

1.1.1 Quantum Fidelity 

The quantum fidelity F(t), being the square of the overlap of the state |V'(0) = t^(OI'0(O)) 
obtained by the unperturbed evolution Uo{t) and the state \ips{t)) = Us{t)\^p{0)) obtained by 
the perturbed evolution Us{t), 

F{t) := immm' = Kvwic/jwt/^wivw)^ (i.i) 

has b een first used as a measure of quantum stability bv IPeresI (|l984l ). see also his book (jPerea . 
I1995I ). The fidelity can also be interpreted as an overlap of the initial state and the echo state 
obtained after forward unperturbed evolution followed by a backward perturbed evolution, 
i.e. after evolution with Ug{t)Uo{t). Quite generally we can imagine unperturbed evolution 
being governed by an unperturbed Hamiltonian Hq and the perturbed evolution by slightly 
perturbed Hamiltonian Hs = Hq + 5V, with 5 a dimensionless perturbation strength and V 
the perturbation operator. Peres reached non-general conclusion that the decay of fidelity is 
faster and has a lower asymptotic value for chaotic than for regular classical dynamic. As we 
will see the general situation can be exactly the opposite. Non decay of the fidelity for regular 
dynamics in Peres's work was due to a very special choice of the initial condition, namely that a 
coherent wave packet was placed in the centre of a stable island. Such a choice is special in two 
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ways, first the centre of an island is a stationary point and second the number of constituent 
eigenstates of the initial sta te is very small. After P eres's work the subject lay untouched 
for about a decade. In 1996 IBallentine fc ZibinI (|l99fil ) numerically studied a quantity similar 
to fidelity. Instead of perturbing the backward evolution (i.e. the Hamiltonian) , they took 
the same backward evolution but instead perturbed the state after forward evolution by some 
instantaneous perturbation, like shifting the whole state by some 5x. They also looked at the 
corresponding classical quantity. The conclusion they reached was that for chaotic dynamics 
quantum stability was much higher than the classical o ne, while for regular dynami cs the two 



u 



agreed. All these results were left mainly unexplained. iGardiner et al. I (|l997l . Il998l ) proposed 
an experimental sheme for m easuring the fidelity in an ion trap. Somehow related to the studies 
of stability was also a work bv lSchack &: Caved ( 1992I . 119931 . Il99a ^. where they studied how much 
information about the environment is needed to prevent the entropy of the system to increase. 
Fidelity studies received new impetus by a series of NMR experiments carried out by the group 
of Pastawski. 



In NMR echo experiments are a standard tool. The so called spin echo experiment of lHahn 

(|l95ri ) refocuses free induction decay in liquids due to dephasing of the individual spins caused 
by slightly different Larmour frequencies experienced due to magnetic field inhomogeneities. 
By an appropriate electromagnetic pulse the Zeeman term is reversed and thus the dynamics 
of non-interactin g spins is reversed. The first real interacting many-body echo experiment was 



done in solids bv lEhim et a,l. 



(|l97nl ). Time reversal, i.e. changing the sign of the interaction, 
is achieved for a dipolar interaction whose angular dependence can change sign for a certain 
"magic" angle, that causes the m ethod to be ca l led m agic echo. Still, the magic echo showed 



strong irreversibility. Much later, IZhang et al. I (199j) devised as sequence of pulses enabling 
a local detection of polarisation (i.e. magnetic moment). They used a molecular crystal, fer- 
rocene Fe(C5H5)2, in which the naturally abundant isotope ^^C is used as an "injection" point 
and a probe, while a ring of protons ^H constitutes a many-body spin system interacting by 
dipole forces. The experiment proceeds in several steps: first the ^^C is magnetised, then this 
magnetisation is transfered to the neighbouring ^H. We thus have a single polarised spin, while 
others are in "equilibrium" . The system of spins then evolves freely, i.e. spin diffusion takes 
place, until at time t the dipolar interaction is reversed and by this also spin diffusion. After 
time 2t the echo is formed and we transfer the magnetisation back to our probe ^^C enabling 
the detection of the polarisation echo. Note that in the polarisation echo experiments the total 
polarisation is conserved as the dipole interaction has only "fiip-flop" terms like S'\_S'!_ , which 
conserve the total spin. To detect the spin diffusion one therefore needs a local probe. With 
the increase of the reversal time t the polarisation echo - the fidelity - decreases and Zhang et 
al. obtained approximatelv exponent ial decay. The nature of this decay has been furthermore 



elaborated in iPastawski et al. I (|l995j). The group of Pastawski then performed a series of NMR 



experiments whe re thev studied in more detail the dependence of the polarisatioii echo on var- 



ious parameters ( Levstein et al. 1 . 119981 : lUsai et al. 1 . 119981 : JPastawski et al. 1 . 12000|) • They were 



able to control the size of the residual part of the Hamiltonian, which was not reversed in the 
experiment and is assumed to be responsible for the polarisation echo decay. For small residual 
interactions they obtained a Gaussian decay while for a larger ones the decay rate saturated 
and was independent of the perturbation strength, i.e. of the size of the residual interaction. 
While there is still no complete consensus on the interpretation of these experimental results 
they triggered a number of theoretical and even more numerical investigations. We will briefly 
list them in chronological order. 

Using the semiclassical expansion of the quantum propagator Ijalabert Sz Pastawskil ( 200ll ) 
derived a perturbation independent fidelity decay for localised initial states and chaotic dynamics, 



F{t) 



-xt 



also called a Lyapunov decay due to its dependence on the Lyapunov exponent A. 
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Perturbation independent fide lity decay obtained for st rong perturbations has also been studied 



numerically in a Lorentz gas ( Cucchietti et al. I . l2002al ^ . the perturbation being in the mass of 



the particle. In the same paper the authors also studied the asymptotic fidelity saturation value 
which, for the strong perturbations considered, is independent of the perturbation strength, 
but becomes perturbation dependent for smaller perturbations (see Section 12.2.11). Studyin g 
fidelity t urned out to be particular ly fruitful in terms of the correlation function (jProsem |2003) , 
see also (|Prosen fc Znidarid 120011 ). For coherent initial states, the fidelity of regular systems 
decays as a Gaussian while for chaotic systems we have different regimes, the most prominent 
being the perturbation depende nt exponential decay. T he dec ay of quantum fidelity in ch aotic 
systems has also been studied bv l.Tacquod et al. I ()200ll ) and bv lCerruti fc Tomsovid (|2002h . For 
sufficiently small perturbations one gets a Gaussian decay (also called a perturbative regime), 
for intermediate ones the so-called Fermi golden rule regime of exponential, perturbation de- 
pendent decay, while for still stronger perturbations we get the Lyapunov decay. All these 
regimes, including the fidelity decay for regular dynamics, an d for different initia, ! state s, were 



Prosen &: Znidarid (200j). Sev 



carefully discussed using the correlation function approach in 
eral interesting results were obtained, perhaps the most surprising one being that in a certain 
range of parameters we can, by increasing chaoticity of the corresponding classical system, in- 
crease quantum fidelity, i.e. improve the stability of quantum dynamics. Different time and 
perturbation scales were discussed as well as their dependence on the number of degrees of 
freedom. It is well known that the quantization of classic al system is no t unique, i.e. there are 



different quantizations leading to the same semiclassics. Q3 li^ compared the quanti- 
zation ambiguity in chaotic and regular systems, reaching a conclusion th at in chaotic systems 



the qu antization ambiguity is supressed as compared with regular ones. iKarkuszewski et al. 
(2002) connected the fidelity decay with the decay of the off-diagonal matrix elements of the 
reduced density matrix, and therefore with decoherence. They claimed that quantum systems, 
whose classical limit is chaotic, are particularly sensitive to perturbations due to small scale 
structures in their Wigner functions. Actually, what their results show, was just the depen- 
dence of the fidelity decay on the size of the initial state. Fidelity decay is faster for random 
initial states than for coherent ones a nd there are no quan t um e f fects in the regime the y stud- 
ied. For related comments see also ( Jordan fc Srednickil . I2OOII : JJacouod et al. I |200j). The 



transition between Fermi golden rule and the Lyapunov decay h as been fu rther considered in 



a Bunimovich stadium fjWisniac ki et al. I . l2002l ) . see also (Cucch ietti et al. I .l2002b) for a study 



of Lyapunov decay. We should stress that the Lyapunov decay of quantum fidelity is purely a 
consequence of the quantum-classical correspondence. There is nothing "quantum" in it and 
can be explained in terms of t he classical fi delity displaying the same perturbation independent 
decay as the quantum fidelity (|Veble fc Prosen. 2004 ) . The fidelity decay c an also be connected 
with the local density of states, although not in a straightforward manner ( Wisniacki &: CohenL 
l20n2l V 

Fidelity decay in mixed systems, having a coexisting regul ar and chaotic componen ts, is 



not as well studied as in purely chaotic or regular situation. In Jwemsteiu .^ .1 1 . 1^?^ ) they 



studied the fidelity decay for initial states placed at the border between regular and chaotic 
regio ns and observed a pow er-law fidelity decay although over less than an order of magni- 
tude. JJacauod et al. I (2003) studied averaged fidelity decay in regular systems for perturba- 
tions with a zero time average for which they predict a universal power-law decay with a power 
3/2. This does not agree with our findings, see discussion at the end of Section 14.2.21 If the 
packet after an echo in the chaotic system drifts exponentially a way from its position at the 



beginning, the fidelity can decay also in a super-exponential way ( Silvestrov et al. 1 . 12003 ) . In 



the same paper the authors also considered the influence of different averaging procedures on 
the decay of short-time fidelity. From the correlation function approach one easily sees that 
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the decay of the fidehty does not only depend on the unperturbed dynamics, say being regular 
or chaotic, but also on the perturbation. If the correlation function decays to zero sufficiently 
fast, we can get exponential decay also for regular unperturb ed dynamics. Such is t he case for 
instance if the perturbation is modelled by a random matrix ( Emerson et al. LI2002I ). but note 
that such perturbation has no direct classical limit. Instead of t he fid elity one can als o study 
the Fourier transformation of the fidelity, see ( Wang &: Lil . l2002h . In ( Wisniacki l2003h it was 
studied how different perturbations of the billiard border influence the short-time fidelity decay. 
For an attempt to giye a uniform approximation of the fidelity decay in the crossoyer regime 
between the Gaussi an perturbatiye and t he exp onential Fern ii golden rule decays using random 



2nn4h. The asymp- 



matrix theory see ( Cerruti fc Tomsoyid . |2ilQ2ilM) and also iGori n et al. 

totic satu ration leyel of the fidelit 

results of ( Prosen fc Znidarin . 120021 ) . Perturbatiye calculation of the fidelity decay in disordered 

systems by field theoretical me thod, i.e. diagrammatic expansion of Green's function, has been 



2002a ). confirming 



done by lAdamoy et al. I ( 200,'j 



Recently IVanicek Sz Hellej ( 20031 ) deyised an efficient numerical scheme for a semiclassical 
eyaluation of the quantum fidelity. The method consists of transforming intractable (due to 
exponentially many contributing orbits) semiclassical expressions in coordinate space into an 
initial momentum space representation. A surprising quantum phenomenon of a prolonged 



Prosen fc Znidarid (2003b|) for regular 



stability, called freeze of fidelity, has been described in 
syster ns. Later, i t has been generalised to arbitrary dynamics, in particular to chaotic sys- 
tems (IProsen &: Z nidaric. 2 0041 ). The decay of fidelity in regu lar one dimensional systems has 



been stuHieH in jRaukarauaravauau X. T.^shmiuarav t:! ^^) and they also numerically Ob- 
seryed a yery short correspondence between the classical and the quantum fidelity for initial 
states placed in a rotational part of phase space, where the ayerage perturbation is zero. This 
is nothing but the freeze of fidelity, not present in classical fidelity. The Lyapunoy regime, 
bei ng of purely class i cal orig in, and its borders of yalidity have been furthermore elaborated 
in ( Gucchietti et al. V I20031J ) , see also ( Prosen fc Znidario . l2002h for a detailed discussion of 
borders within which different regimes are yalid. They also stressed the importance of non- 
commutativit y of the limits fi — > and 5 — > in recoyering the classical behayiour, as already 
explained in ( Prosen &: Znidarid. 2002^ . A nontriyial question, addressed by Hiller, Kottos, 
Gohen and Geisel (jHiller et al. Ll2004^ concerns the optimal time of the unperturbed eyolution, 
i.e. we fix the duration of the unperturbed eyolution and seek the duration of the perturbed 
eyolution for which the fidelity will be maximal. The quantum fidelity decay in weakly chaotic 
systems when one might not get the Lyapunoy or the Fermi golden rule decay has been ex- 
pl ored in (IWang et al. 1 . 12004^ . T he fidelity for various rand om matrix models has been analysed 
by iGorin et al. I (|2QiJJ) , see also ( Gerruti fc Tomsovid . l2003b fai) . For further random matrix re- 



Dy n.Torm et al. I Ij^uu4i i , see also iji.erruti <v: lomsovid . iztni.-SDiif^). .for lurtner rando m ma trix re- 
sults and their relevance for quantum computation see also (jFrahm et al. L |200J) . In ( lominL 
1200,31 ) the fidelity of a nonlinear time-dependent chaotic oscillator with respect to the time- 
dependent perturbation of its frequency has been studied analytically. Most of the s tudies 
so far focused on a few degrees of freedom systems, the exce ption being ('Prosen, '2002*) . For 
addit ional results on the stability of many-body systems see ( Izrailev fc Gastaheda-Mendozal . 

ET" 



There has also been a large number of papers dealing with the stability of quantum computa- 
tion, i.e. studies of the fidelity for specific quantum algorithms and perturbations, mainly relying 
on numer ical simulations. Some of these include ( Miguel et al. V 119961 . 119971 : iGea-Banaclochd . 
119981 . Il999 . 2000) . The group o f Berman studi e d in d etail t he stability of the Ising quan- 



tum computer, see for in stance (|Berman et al. 



also (JGelardo et al. 



2001 



2002 



and references therein and 
I) . The group of Shepelyansky studied the stability of rnany different 



quantum algorithms (jGeorgeot &: Shepelyanskvll200d : ISong &: Shepelyanskyll200ll : lBenenti et al. . 
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bond : lTRrra,nRn fc ShenRlvajiskvi i2nn.liV srr r1so (JReneTiti p.t al. I . 1200301 and (jRettRllii . l2004f r 
The scaling of errors with various parameters is easily explained by our correlation function 
approach. 

With the advances in experiments manipulating individual quantum systems it has become 
possible to actually measure quantum fidelity. Apart from NMR echo experiments already 
menti oned, one is able to mea sure the quantum fidelity on a few qubit NM R quantum com- 
puter (TWeinstein et al. I . l2002bi l. Particularly pro mising candidate s are ion traps Gardiner et al. I 



juter (IWei 



( 19971 . Il99 8h , ultra cold atoms in optical traps ( Andersen et al. I l200.'J : 



Experiments with microwave resonators in billiards are also under way ( 



Schlunk et oi. 1. 1200, 



Schafer et al. 



EST 

2003). 



1.1.2 Classical Fidelity 

Cl assical fidelity can be defin ed in an analogous way as quantum fidelity and has been first used 



m 



( Prosen fc Znidarid . l2002l ) . It is an overlap integral of two classical densities in phase space, 



obtained by unperturbed and perturbed evolutions, 

-^clas(t) := dxpo{x,t)ps{x,t), (1.2) 

where ps{x,t) is the density in phase space obtained by a perturbed evolution 4>s{x), where 
(f)s{x) is a volume preserving flow in phase space. The density at time t can be obtained 
by backward propagating phase space point x, ps{x,t) = p{4>s (^))0)- Note that in order for 
fidelity to be normalized to 1, the classical density has to be square normalized, Jdxp'^{x, 0) = 1. 

Classical fidelity behaves markedly different from quantum fidelity. For chaotic systems and 
localized initial states the classical fidelity agrees with quantum only up to the short Ehrenfest 
time, logarithmic in Planck's constant, when the quantum-classical correspondence breaks down. 
For regular systems and if the time average perturbation is nonzero though, classical fidelity 
follows quantum fidelity. On the other hand, if the time averaged perturbation is zero, quantum 
fidelity exhibits the so called freezing (Chapter^)), while classical fidelity does not, except in the 
non-generic case of a harmonic oscillator (Section 15.5(1 . The classical fidelity will depend on the 
stability of orbits in phase space. Linear response will therefore depend on the stability of the 
flow and this involves the derivatives of the flow, a derivative being an unbounded operator. For 
chaotic systems for instance, the derivatives grow exponentially in time due to orbit separation, 
a simple consequence of the famous Lyapunov instability. In a phase space picture, the classical 
dynamics can produce structures on an arbitrary small scales, even if we start from a smooth 
density. Resolution of the quantum mechanics on the other hand is limited by a finite Planck 
constant. 

We will give a brief overview of the known results about classical fidelity. This will help us 
to understand the differences with quantum fidelity which will be the object of study in the 
present work. The literature on classical fidelit y is not nearly as extensive as on quantum fidelity. 



Linear response calculation has been done in (|Prosen fc Znidaric. 2002 ). Numerical results on 



the classical fidelity and its correspondence with the quantum fidelity in chao tic systems and in 



systems exhibiting diffusion have been presented in (JBenenti fc Casatf . 2002). Classical fidelity 



in regular and chaotic systems has also been theoretically discussed in (jEckhardi:, 20 03). A de- 



tailed explanation of the asymptotic decay in chaotic systems has been given i n ( Benenti et 



2003bri and a theoretical explanation of the Lyapunov decay for short times in (IVeble &: Prosen , 
20041'). Classically regular systems on the other hand have been worked out in (|Benenti et al. 
2003ah. 
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Chaotic Systems 

The classical fidelity in chaotic systems will go trough different decay regimes as time increases. 
We will consider a localized initial state of size v. St arting f rom t = fidelity stays close to 1 



until time ty (JBenenti fc Casatil . 12002 : Benenti et al. I . l2003bt ) 



u 






(1.3) 



where A is the Lyapunov exponent. The time ty can be thought of as the time in which the initial 
perturbation is "amplified" to the size of the initial packet and ps starts to differ from pQ. For 
a rigorous derivation of ij^ as well as for a discussion of multi-degree of freedom systems, where 
there is a cascade of times t,y, see ( Veble Sz ProsenL l2004 ) . After t^, the so called Lyapunov decay 
F{t) = exp {— A(t — t^)} sets in and lasts until time ts determined by the spreading of the packet 
over the whole phase space. The Lyapunov decay has been explained bv lVeble Sz ProsenI (200J) 
and for systems with more than one degree of freedom there is a whole cascade of Lyapunov 
decays determined by the Lyapunov spectrum*. If we have a classical system with diffusion, 
such that the phase spac e is much larger in one direction, say q,p £ [0, 2tt] x [0, L], with L ^ 2tt, 



time ts will be given by ( Benenti et al. I . l2003bl ) 



ts 



1 , 27r 



(L4) 



After ts fidelity will decay diffusively F{t) ~ l/\/Dt with the diffusion constant D, until the 
diffusive process reaches the phase space boundary also in p-direction, i.e. at ^d ~ L'^/D. 
Note that this diffusive regime is present only if L ^ 27r. After to the asymptotic decay of 



classical fi delitv begins. T^^^ asy mptotic decay is determined by the largest Ruelle-Pollicott 

jg). i.e. the eigenvalue of the Perron- Frobenius operator, and is 



resonance ( Benenti et al. 



20031 



thus the same as the asymptotic decay of classical correlations. If there is a gap in the spectrum 
of the Peron-Frobenius operator, this decay will be exponential, otherwise it can be power law. 
Note that this asymptotic decay rate does not depend on the perturbation strength 6 but only 
on the phase space size L. For large times fidelity will decay towards the asymptotic value 
F = F{t — > oo) determined by the ratio of the initial packet size and the phase space size. In 
order to see the asymptotic regime of fidelity decay one has to look at F(t) — F. Furthermore, 
for the fidelity at the end of the Lyapunov decay ts to be larger than F we must have 5 > F. 

The decay rate for classical fidelity in chaotic systems, apart from the non-decaying "shoul- 
der" until tu, does not depend on the perturbation strength 6 (decay time borders do depend 
though) and is therefore independent of the perturbation itself. This must be contrasted with 
quantum fidelity, which as we will see does depend strongly on perturbation strength and type. 



Regular Systems 

The classical fidelity decay for regular systems has been explained ( Benenti et al. I . l2003al ) by 



studying changes in the action-angle variables caused by the perturbation. There is a compe- 
tition between two contributions to the fidelity decay. It can decay as a consequence of the 
ballistic separation of the perturbed and unperturbed packets caused by different frequencies of 
the perturbed and unperturbed tori on which the initial packet is placed, or it can decay due 
to different shapes of unperturbed and perturbed tori. In the latter case the decay is caused by 
transitions of the perturbed packet between unperturbed tori. Therefore, if the perturbation 
predominantly changes the frequency of tori, the classical fidelity will exhibit a ballistic decay. 



*In case of drift of packets one can get a super-exponential instead of exponential decay of fi delity, 
sufficiently small times quantum fidelity will show the same phenomena, see llSilvestroy et al. 1120031) . 



For 
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The shape of this balhstic decay is determined by the shape of the initial packet. If the packet is 
a coherent state having a Gaussian shape, the decay will be Gaussian. On the other hand, if the 
perturbation predominantly changes the shape of tori the decay will be algebraic F{t) ~ l/(6t)'^, 
in a system with d degrees of freedom. Which type of decay one gets, depends on the shape of 
the perturbation but not on its strength 6, provided it is small enough. The type of the decay 
can also depend on the position of the initial packet. 

Note that the above results were derived under the assumption that the angles across the 
packet are spread over 27r, i.e. t > 2TT/{Ajduj/dj), where Aj is the width of the initial packet in 
the action direction. Incidentally, this time is equal to the time ti (|4.50j) after which quantum 
fidelity freezes at a constant value if the time averaged perturbation is zero (see Section l4.2|l . 
Furthermore, the change of the actions across the packet caused by the perturbation must be 
small, i.e. Aj ^ 6. This last condition translates for coherent packets into 6 <^ h ' , also being 
the upper limit of validity of approximations used in the calculation of the quantum fidelity 
plateau (HSU- 

Comparing to quantum fidelity, the case of ballistic decay of classical fidelity corresponds to 
perturbations having a nonzero time average, discussed in Chapter |21 For such perturbations 
quantum fidelity agrees with the classical one under certain conditions. The case of the algebraic 
decay corresponds to the perturbation with a vanishing time average. In this case though, the 
quantum and the classical fidelity do not agree. What is more, quantum fidelity displays an 
intriguing new feature called freezing (see Section [4. 2|) and decays only on a much longer time 
scale ~ 1/5"^. Whereas only the functional dependence of the classical fidelity decay changes 
depending on the perturbation type, always decaying on a time scale ~ 1/(5, the decay of 
quantum fidelity drastically changes if one has a perturbation with a zero time average. 



1.1.3 Entanglement 

The literature on decoherence is exhaustive and we will here list only those more or less directly 
related to our work. 

Time independent perturbative expansion has been used bv lKiibler fc Zehl (|l97,'j ) to study 
the eigenva lues of th e reduced density matrix. A similar expansion was used much later for 
the purity (|Kim et al. . 1996). Note that these perturbative approaches are not equivalent to 
our linear response expressions presented in Chapter El as they correspond only to the short 
time regime in which the correlation function is constant. Among the early studi es of entropy 



growth and its relation t o the chaoticity of the underlying system is the one by lAlicki et al. 



( 1996| ). iMiller &: Sarkai (Il999l'l observed a linear entro py growth with the slope given by the 
Lyapunov exponent. ICorin fc SeligmanI (120021 . 1200,^ studied the purity decay for random 
matrix models. R andom matrix asymptotic value of the pur ity (or of linear entropy) was 
later re-derived in JBandvopadhvav & Laks hminaravanl . l2002lJ ) and numerically observed in 
chaotic systems. iProsen &: Seligman (2002 ) first defined the purity fidelity, generalising the 
purity to echo dynamics. The pur ity fidelity and its rel ation with the fidelity in chaotic and 
regular systems was stu died also in (IProsen et al. L l2003bh using a correlation function approach 
which was used also bv lTanaka et al. I (|200j) for studying chaotic systems. The connection be- 
tween th e purity and the fidelity was put on a firm gro und by a rigorous inequality between 



the two (jZnidaric fc Prosen . 20031: 



Prosen et al. 



2003a|). The classical analog of decoherence 



has been studied in ( Gong fc Brumeii l2003d |bllah . Using a perturbative approach the infiu- 



ence of the type of the pertur bation and of t he dy namics on the quantum-classical correspon- 

J, |20QJ) for a study of quantum-classical correspon- 



dence were explored, see also (jAngelo et al. 

de nce of entanglement . The entangleme nt in weakly coupled composite systems was studied 



in (jZnidaric fc Prosem l2003h as well as ( Fuiisaki et al. I . l2003l ). reaching the same conclusion, 
namely that the increase of chaos can inhibit the production of entanglement. Subsequently 
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the inequaUty between fideUty and purity was used in (ICucchietti et al. I. 2nf)3a[^ . Entanglement 



if 
production in a coupled baker's map has been studied in (IScott &: Caved . 120031) and in coupled 

kicked tops in ( Bandvopadhvav fc Lakshminaravanl . |2i]i}2^ . Similar semiclassical methods as 

for the evaluation of the fidelity have been used also for the purity (Jjacauod, 2004), predicting an 

exponential decay in chaotic systems, confirming prediction in (jZnidaric fc ProsenL l2003h . and 

algebraic decay in regular systems. The predicted power of the algebraic decay though does not 

agree with our theory and numerical results presented in Section f^.S.ll The entangleme nt under 



echo situation during a quantum computation has been studied i n (JRossini et al. ^ 2003h. Entan- 



glem ent in weakly coupled kicked tops has been studied also in ( Demkowicz-Dobrzahski fc Kua . 
I2OO4I ) . among other things also for random initial states. 



1.2 Outline 

In Chapter [21 we introduce the quantum fidelity and briefly review numerical models used for 
the illustration of the the ory. The material in sect ions d ealing with the avera ge fidelity has 
been mainly published in ( Prosen fc Znid arid. l2002h and (jProsen et al. L l2003d ). The core of 
Chapterl^Jhas been published in (|Prosen fc Znidaric, 200J), which is the main paper presenting 
the correlation function approach to fidelity. New, unpublished material consists of numeri- 
cal demonstration of a Gaussian distribution of diagonal matrix elements in chaotic systems, 
relevant for a Gaussian perturbative decay of fidelity. New is also the discussion about the 
average fidelity in regular systems as well as the last part of Section 13.31 with the illustration of 
the fidelity decay in terms of Wigner functions and some additional figures hopefully clarifying 
the relations between differen t dec av regimes. The ma,terial of Chapter [l] has been published 
in ( Prosen &: Znidarid . |2003b|) and ( Prosen fc Znidarid . |200J) . We unified the approach in mix- 
ing and regular situation making the exposition of a regular case more concise. New material is 
the explanation of the average fidelity decay in regular systems with zer o time average pertur- 
bation. About hal f of th e contents of Chapter 151 has been published in (|Prosen et al. 

i). 



Znidaric &: Prosenl . 120031 : JProsen et al. 



f )ertur- 
2003bl: 



2003s 



The new material consists of the purity fi- 
delity and the reduced fidelity calculation in regular systems beyond linear response, most of 
the discussion about the Jaynes-Cummings model, the section describing the freeze in a har- 
monic oscillator and the last section explaining the accelerated decoherence of cat states. The 
application of the fi delity theory to the improvement of the quantum Fourier transformation has 
been published in ((Prosen &: Znidaria 2001 ) an d an extension to the Ising model of quantum 
computer, not presented in the present work, in ( Celardo et a.l. I . l2003h . 



Chapter 2 

Fidelity 



Basic research is what I am doing when I 
don't know what I am doing. 

— Werner von Braun 



Quantum fidelity F between two general densi ty matrices p a. nd a, representing either a 
pure or a mixed state, has been used by Uhlmann ( UhlmannL Il976l ) . 



F{p,a) = [tr^ py^apy^y. (2.1) 



He called it a transition probability and the name fidelity was introduced by Jozsa ( Jozsal . |l99J) 



It is symmetric with respect to the exchange of p and a and in the case of one density matrix 
being a pure one, the general expression for fidelity simplifies into F = (iplalip) and if both are 
pure it is 

F{i;,ip) = Mip)\\ (2.2) 

where obviously p := \ip){ip\ and a := \{p){ip\. We will always deal with pure states and the 
latter definition (|2.2|) will be sufficient for us. Note that sometimes the name fidelity is used for 
a quantity without a square, i.e. for a fidelity amplitude. 

We will study stability of quantum dynamics with respect to the perturbation of evolution 
and the quantity studied will be the fidelity between states obtained by unperturbed and per- 
turbed evolutions, starting from the same initial state. Let us denote the initial state by |^(0)), 
and the states at time r as 

IV^(r)) = C/o(t)|V(0)), \Mr)) = Us{r)mO)), (2.3) 

where Uo{t) is a unitary propagator from time to time r and Us{t) is a perturbed propagator. 
The overlap between the perturbed and unperturbed states, denoted by F{t), will serve us as 
a criterion of stability, 

F{r) = \{Mr)mr))\' = |/(r)P, /(r) = (^(0)|M,(r)|V(0)), (2.4) 

where we introduced a complex fidelity amplitude /(r) and the unitary echo operator 

Ms{t) = uI{t)Us{t). (2.5) 

The fidelity as defined in (|2.4|) is a real quantity between for orthogonal states and 1 iff the 
two states are equal (up to a phase) and is a standard measure of stability. 
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Up to now we have not specified the propagators Uo{t) and Us{t) and to simphfy the matter 
we will limit ourselves to the case where a propagator for time r can be written as a power of 
some basic single-step propagator Uq. That is, we will use time index t measuring the number 
of basic units of duration tq, t = tTQ. We therefore write a single step propagator Uq = Uo{tq) 
and similarly for the perturbed evolution Us = Us{tq). The propagator for t steps is now simply 
the t-th power of a basic propagator Uo{Tot) = Uq. These discrete time formalism allows us to 
treat two interesting cases at once: the case of time-independent Hamiltonian and the case when 
the Hamiltonian is time periodic function H{t) = H{t + tq) with a period tq. The latter case 
includes the so-called kicked systems which will be used for the numerical demonstration of our 
theory. Note that the theory can be easily generalised also to time-dependent Hamiltonians, only 
notation becomes more cumbersome. For an example see Section [6.21 describing the application 
of fidelity theory to the quantum Fourier transformation algorithm. The general perturbed 
evolution can be written in terms of the perturbation generator V as 

Us = Uoexp{-[V5To/fi), (2.6) 

where 5 is a dimensionless perturbation strength. The above equation 1)2. 6() can be considered 
as a definition of a hermitian operator V, given the unperturbed and perturbed one time step 
propagators Uq and Us, respectively. In the case of Hamiltonian dynamics with Hq generating 
unperturbed evolution as Uq = exp (—iHotTo/h) and with the perturbed Hamiltonian of the 
form Hs = Hq + H'5 we have V = H' + 0{tq6). The difference between H' and V therefore 
goes to zero in the limit of either small perturbation strength 5 or in the limit of small time step 
tq. The latter limit, namely tq — > 0, corresponds to the case of time-independent Hamiltonians, 
where the time step tq can be chosen arbitrarily small. From now on we will drop the irrelevant 
parameter tq (i.e. take a unit of time to be tq) in all equations, so that the discrete time 
index t has units of time. We will use the same letter t for a discrete time index as well as 
for a continuous time (limit tq -^ 0) on few occasions. Whether t is a discrete time index or 
a continuous will be clear from the context. Our definition of fidelity (|2.4|) is in discrete time 
formulation 

F{t) = \{^mMs{t)\i,{m\ Msit) := Uq'uI (2.7) 

The fidelity is just the expectation value of the echo operator. It can be equivalently expressed 
in terms of the initial pure density matrix p{0) and the echo density matrix p (t) (sometimes 
refered to as the Loschmidt echo) obtained from p(0) by evolving it with the echo operator 

Ms{t) 

p{0) = mO))mO)l P'^it) ■■= Ms{t)p{0)Ml{t), (2.8) 

as 

F{t) = tr[p(0)p^(t)], fit) = tT[p{0)Ms{t)]. (2.9) 

On several occasions we will be interested in the fidelity averaged over some ensemble of initial 
states. In such a case the density matrix p(0) in the above definitions of fidelity and fidelity 
amplitude f{t) (|2.9|) must be replaced by appropriate mixed density matrix. We will frequently 
use a uniform average over whole Hilbert space of dimension M, p{0) = 1/M, see Section 12.2.21 
for details. 

Using our definition ()2.6() of perturbed dynamics in terms of the perturbation generator V, 
the echo operator can be rewritten in a more useful way by writing V in the interaction picture*, 
i.e. for an arbitrary operator A its interaction picture A{t) is 

A{t) = Uq'AUI (2.10) 



'It can also be considered as the Heisenberg picture of the unperturbed dynamics. 
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and so the echo operator is 

Ms{t) = exp {-iV{t - l)6/h) • • • exp {-iV{0)6/h) = T exp {-i^{t)6/h), (2.11) 

where T is a time-ordering operator and S(t) is the sum of operators V{j), 

S(t):=|]y(j). (2.12) 

In the case of the time-independent Hamiltonian one has S(f) = /q V{t')dt' . The form of 
the echo operator (|2.11j) is nothing but the interaction picture of the perturbed propagator, 
famihar in quantum field theor y. Methods of quantu m field theory have been actually used 
by Adamov, Gornyi and Mirlin ( Adamov et al. 1 . 1200,^ to calculate the fidelity in a disordered 



system. Remember that the time dependence of V{t) comes from the interaction picture 1)2.10(1 
and so V{t) is time dependent even if original perturbation V ()2.6|) is not, which is the case 
throughout our derivations. The same definition of fidelity, just by re placing density matrices 
with the densities in phase space, can be used for classical fidelity, see lVeble fc Prosen (2004 ). 
The echo operator Ms{t) can we written as an exponential function of a single operator by 
using the Baker-Campbell-HausdorfF (BCH) formula 

e^e''=exp^A + B + ^[A,B] + ^[A,[A,B]] + ^[B,[B,A]] + ...y (2.13) 

Applying the BCH formula on the product form of Ms{t) ((2.111) gives to the order (5^ in the 
argument of the exponential function 

Msit) = exp |-^ (^S(t)5 + ^r{t)5^ + ■■■)}' (2-^^) 

with the hermitian operator T(t) being 

rW = ^EE[nj),nfe)]- (2.15) 

j=0 k=j 

Both operators $](t) and T(t) standing in the expression for the echo operator have a well defined 
classical limit, provided V has a classical limit. The classical limit of T{t) can be obtained by 
replacing commutator with the Poisson bracket {•,•}, (— i/7i)[»,»] — > {•,•}. The BCH form 
((2.14() of the echo operator will be particularly useful in the case of regular dynamics. For now 
let us list three different possible fidelity decays depending on the behaviour of operator S(t). 
For mixing dynamics the fluctuations of $](t) give the dominant contribution, i.e. terms like 
(S'^(t)) grow linearly with time as ~ t. For regular dynamics second moment (E^(i)) ~ t^ 
grows quadratically with time, corresponding to the existence of a nontrivial time-averaged 
perturbation V(t). In certain cases we can have even (E^(i)) ~ t^ which results in the so-called 
freeze of fidelity. Only in this last case is the fidelity decay caused by the operator T{t). 

If we want the fidelity expressed as a power series to all orders in 6 the BCH formula becomes 
too cumbersome. An easier approach is to just expand the product form of the echo operator 
Ms{t) ((2.11(1 . giving the fidelity amplitude 

oo / ■\ms:m t—1 

/W = i+EMfi^ E 'r{v{n)v{j2)---v{j^)), (2.16) 

m=l JlvJm = 

where we introduced the notation (•) = (V'(0)| • |^(0)). One can see, that the fidelity is 
expressed in terms of ?7i-point quantum correlation functions of the perturbation generator V. 
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In the semiclassical hmit one is able to replace quantum correlation functions with the classical 
ones and therefore the quantum fidelity is expressed in terms of classical quantities. One should 
keep in mind that the classical fidelity behaves distinctively different from quantum fidelity (see 
Chapter ll.l.2() and therefore it is not obvious that quantum fidelity can be expressed in terms 
of classical quantities. Power series expansion of the classical fidelity for example does not yield 
a fruitful result (for chaotic systems) due to unboundedness of the classical operators and the 



classical fidelity can not be expressed in a similar way as the quantum fidelity (|2.16|) . 

Fidelity F{t) is now obtained by taking absolute value square of f(t) ()2.16|) . Only even 
orders in 6 survive as odd orders just rotate the phase of the fidelity amplitude f{t), the lowest 
order being quadratic in 6 

P(t) = 1 - S E C{j, k) + 0{6') = 1-^ {(S2(t)) - (S(t))n + 0{6% (2.17) 

"- j,fc=0 "" 

with C{j, k) being the quantum correlation function 

C{j, k) = {V{j)V{k)) - {V{j)) {V{k)) . (2.18) 

The second order expansion of the fidelity ()2.17|) is one of the central theoretical results. Al- 
though it is very simple and is just the lowest order expansion it contains most of the essential 
physics of fidelity decay. Furthermore, there is one very pragmatic reason why it is sufficient and 
higher orders are not needed. In all practical and experimental applications where the fidelity is 
the relevant quantity, one is mainly interested in a regime of high fidelity, i.e. of high stability. If 
the interesting range is say F(t) > 0.9, higher orders will give only corrections of order 0.01. To 
see this, let us denote the second order term with x = {6/h)'^ ^j,k C'(i> k). From the expansion 
(|2.16|) one can see that the term with 6'^"^ will be at most of the order [St/ti)^^ ~ x^"^. So if 
1 — F{t) ~ x^ <C 1 one can safely neglect higher orders in 5, i.e. terms of order x^™ with m > 1, 
irrespective of the values of individual parameters like 5, fi or time t. The range of validity of 
the lowest order expansion in (|2.17j) is limited only by the value of 1 — F{t). Furthermore, in 
certain cases such as mixing dynamics f Section 13. 1() . regular dynamics (Section 13. 2() or in the 
so-called freeze of fidelity (Chapter |1J) the series (|2.16|) can be resummed to all orders in 5 and 
one gets an expression for the fidelity valid in the whole range from 1 to its asymptotic value. 
Let us now discuss in more detail the lowest order term in the expansion of fidelity which 
can be in the case of continuous time-independent Hamiltonian written as 

1 - F{t) = ^ /7* C{t\ t") dt' dt". (2.19) 

fi Jo Jo 

This linear response expression is reminiscent of Green-Kubo like formulas. It says that the 
"dissipation of quantum information" 1 — F(t) equals the double integral of the correlation 
function of the perturbation V. An interesting and somehow counterintuitive conclusion can be 
drawn from it, namely, the smaller the integral of time correlation function the higher fidelity 
will be. In the semiclassical limit the quantum correlation function approaches the classical one, 
provided the initial state |V'(0)) also has a well defined classical limit, and one will see a fast 
decay of the correlation function for chaotic systems while for regular systems the correlation 
function typically will not decay. The double integral of the correlation function will therefore 
grow as oc t for chaotic systems and like tx t^ for regular systems. The fidelity will in turn 
decay more slowly for chaotic systems than for regular ones. Or in other words, the more 
chaotic the systems is, the slower decay of quantum fidelity it will have, i.e. the more stable 
it is to perturbations. This must be contrasted with the behaviour of the classical fidelity (see 
Section ll.l.2|) which is just the opposite. The more chaotic the system is, the faster classical 
fidelity will decay. 
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2.1 Models for numerics 

We will always use numerical simulations to compare with the theoretical derivations. Two 
classes of systems will be used. One group will be the so-called kicked top models, describing 
dynamics of a spin "kicked" by some external field. In the case of a kicked top the Hamiltonian 
is time periodic and the propagator Uq or Us represent a Floquet map over the period of one 
kick. The second example will be a Jaynes-Cummings system familiar in quantum optics. The 
Jaynes-Cummings model is a time-independent model with two degrees of freedom representing 
a harmonic oscillator (boson) coupled with a spin (fermion). In the following two subsections 
we will present both models and briefly describe their properties, while the actual numerical 
results will be presented in each chapter as needed. 

2.1.1 The Kicked Top 

The kicked top has been introduced by Haake, Kus and Sch arf (IHaake et al. . 1987i'l and has 



ine KicKea top nas been mtroaucea by iiaaKe, n.us ana ncn ari una 
served as a numerical model in numerous studies since then (IHaakd. 



1991 



Shack et al. 



Fox fc Elstoil . Il994l : \A]\ck\et al\ ll 994 JMiller fc SaxkH ll 999: 'Rrcslin fc M ilburTil . Il 9991 1 . In 



19941: 



addition the kicked top might also be experimentally realizable (, Haaka . i200Qil . 

In Chapter |3J discussing the fidelity decay for general perturbations, we will use a unitary 
one step propagator 

'2S 



Uq = C/(7, a) = exp (-i7Sy) exp -ia^ ) , (2.20) 



where 5'x,y,z are standard spin operators [Sk-,Si\ = iSkimSm and a, 7 are two parameters de- 
termining dynamical properties. The propagator for t steps is Uq. Half-integer (integer) spin 
S determines the size of the Hilbert space M = 2S + 1 and the value of the effective Planck 
constant h = 1/5. The perturbed propagator Us is obtained by perturbing the parameter a 

Us = U{j,a + 6), (2.21) 

so that the perturbation generator V is 

In the classical limit S ^ 00 the area preserving map corresponding to Uq can be obtained from 
the Heisenberg equations for spin operators, S(l) = UqSUq. The classical map is most easily 
written in terms of a unit vector on a sphere r = (x,y, z) = S/S as 

x' = C0S7 (xcos (az) — y sin(a2;)) + zsin7 

y' = y cos (az) + xsin (az) (2.23) 

z' = zcos'y — sinj {xcos{az) — ysm(az)). 

The classical perturbation generator p. 22)1 is simply ' V ^ v = z'^/2. The angle 7 in the 
propagator is usually set to tt/2 whereas we will use two different values, 7 = 7r/2 and 7 = 7r/6. 
For these two values of 7 the classical correlation function displays two different decays towards 
zero, a monotonic decay for 7 = tt/6 and an oscillatory decay for 7 = 7r/2. 

The phase sp ace of the classical map 1)2. 24() is regular for small values of a, at a ~ 3 (see 
e.g. ( Perea . Il99,^ ) most of tori disappear and for still larger a the system is fully chaotic. 



'We will use lowercase letters for classical quantities. 
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For 7 = 7r/2 the unperturbed propagator Uq (as weh as the perturbed one) commutes 
with the operator i?y of a vr rotation around y-axis, Ry = exp (— ivrSy). In addition there is 
an antiunitary symmetry so that t he Hilbert sp ace is reducible to three invariant subspaces. 
Fohowing notation in Peres's book ( Pered . Il995l ) we denote them with EE, 00 and OE with 



the basis states (here we assume S to be even) 

EE: |0),{|2m) + |-2m)}/V2 AAee = 5/2 + 1 

00: {\2m-l)-\-{2m-l))}/V2 Noo = S/2 (2.24) 

OE: {|2m) -I -2m)}/V2,{|2m-l) + | -(2m-l))}/V2 AAqe = 5, 

where m runs from 1 to S/2 and \m) are standard eigenstates of S^. For 7 7^ 7r/2 the subspaces 
EE and 00 coalesce as Ry is the only symmetry left and we have just two invariant subspaces. 
Except if stated otherwise the initial state will always be chosen from subspace OE (i.e. initial 
coherent state will be projected onto OE subspace) so that the size of the relevant Hilbert space 
win be AA = 5. 

In Chapter |1J dealing with the so-called freeze of fidelity, we will take a slightly different 
form of the propagator which will be presented in the mentioned chapter. Apart from the one 
dimensional kicked top we will also use a system of two coupled kicked tops to demonstrate the 
dependence of the fidelity decay on the number of degrees of freedom. The one step propagator 
for two coupled kicked tops is chosen to be 

Us = C/i(7, a)f/2(7, a) exp (-i(,5 + £)V/n). (2.25) 

Propagators Ui and U2 are the usual single kicked top propagators (|2.2()j) acting on the first and 
the second top, respectively, and the last term with operator V is responsible for the coupling 
of strength e for unperturbed evolution and e + 5 for perturbed one. The operator V will be 
left unspecified for now as we will use different l/'s. 

The initial state |^(0)) used for fidelity evaluation will be either a random state with the 
expansion coefficients Cm = {Tn\ip{Q)) being independent Gaussian complex numbers or a coher- 
ent state. A random state might be the most relevant for quantum computations for instance 
as it contains the most information and the states used in quantum computation are expected 
to be "random". Coherent state centred at the position r* = (sin t?* cos (^* , sin t?* sin (/?* , cos t?* ) 
is given by expansion 

^ / 2^? \'/' 
|r, (/?*)= y] L cos^+™(r/2)sin^-™(r/2)e-''"^*|m>- (2-26) 



Equivalently, it can be written in terms of a complex parameter r as 



l^*''/'*) = ,, , ,,„.g exp(Tg,)|5), r = e''^*tan(r/2), (2.27) 



with S± = Sx ± \Sy. In the semiclassical limit of large spin S the expansion coefficients of 
coherent state Cm go towards Cm x exp (— S'(?n,/5 — z*)^/2(l — z*^))e~""'^*. Coherent states 
have a well defined classical limit and this enables to compare the quantum fidelity for coherent 
initial states with the corresponding classical fideli t y. Th e initial classical phase space density 
corresponding to a coherent state is ( Fox &: ElstonL |l99J) 



2-5^ r -■■'" n*N2 , /,. ,.*^2„■ 2, 



PcxU^, ^) = \1 — exp {-S[{^ - ry + (99 - ipy sln^ ^]}. (2.28) 

The above density is normalised as Jp'^i^dQ = 1. 
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2.1.2 The Jaynes-Cummings model 



The Jaynes-Cummings model ( Javnes fc CummingsL Il96,4 iTavis fc Cumminga . 119681 ) , see also 



( Mevstre fc Sargent III Il990l . p. 336), is a system of a coupled harmonic oscillator and a spin. 



It can be realized experimentally in the cavity electrodynamics experiments (QED) by sending 
a beam of atoms trough a cavity. The electromagnetic field in the cavity is quantized with 
the Hamiltonian fiujo^a, and the spin degree of freedom fieS^z, of atoms interacts with the 
electromagnetic field. The dominant interaction is a dipolar d ■ E, with the monochromatic 
field E = eEQ{ae — a~^e~ ). For two level atoms only y-component of the dipole moment d 
is nonzero and is proportional to Sy = —^(5'+ — 5_). All this results in the Hamiltonian 

H = huja+a + heS^ + ^L {G {aS+ + a+5_) + G' (a+S+ + aS-)} , (2.29) 

\'2S 

with boson lowering/raising operators a,a'^, [a,a~^] = 1. Planck's constant is chosen as ?i = 1/S 
so that the semiclassical limit is reached for S ^ oo. If we are close to a resonance lo = e 
the rotating-wave approximation can be made in which the fast oscillating term G' can be 
neglected and only the G-term is retained (oscillating with small cu — e). In our discussion 
we will predominantly focus on this situation of G' = 0. In addition, frequencies uj and e are 
usually in a GHz regime, while the coupling G is of the order of kHz and is therefore small. If 
either G = or C = the model is integrable with an additional invariant being the difference 
or the sum of the spin and boson (harmonic oscillator) quanta. 

The initial state will be always chosen as a product state of spin and boson coherent states. 
The spin coherent state is the same as above (|2.26|) . while the coherent state for a boson is 

oo k 

\a) = e"'^^-*'^|0) = exp {-\a\'/2) J] ^\k), (2.30) 

with a being a complex parameter determining the position of the coherent state and |A;) is an 
eigenstate of operator a^a having k boson quanta. The classical density corresponding to the 
boson coherent initial state with a complex parameter a = a^ + iai is 

Pdass(g,p) = y^exp|-5 (q - .JyS a.f + {p - ^/yS aif |. (2.31) 

The normalisation is as usual Jp'^i^^dqdp = 1. 

The classical Hamiltonian is obtained by taking the limit S ^ oo and defining new canonical 
quantum operators, qi and pi for a boson and q2 and p2 for a spin 

(2.32) 

They satisfy commutation relations [91,^1] = ih and [(721^2] = ^^- In the limit 5 — > 00 they 
commute and can be replaced by the classical variables resulting in the classical Hamiltonian 




-f^class = -^(Pi+qi) + ep2 + G+sJ I - plqicosq2 - G-\l I - plpismq2 - — , (2.33) 

withG± = G±G'. 
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2.2 Average Fidelity 

Sometimes the average fidelity is of interest, i.e. the fidelity averaged over some ensemble of 
initial states. Such an average fidelity is also more amenable to theoretical treatment. Easier 
to calculate is the average fidelity amplitude f{t) which is of second order in the initial state 
1^) while the fidelity F{t) is of fourth order in j-i/;). We will show that the difference between 
the average fidelity amplitude and the average fidelity is semiclassically small. 

In a finite Hilbert space the fidelity will not decay to zero but will instead fluctuate around 
some small plateau value. The value of this plateau equals to a time averaged fidelity. For 
ergodic systems this time averaged value equals to the phase space averaged one. 

2.2.1 Time Averaged Fidelity 

We want to calculate the value of fidelity in the limit t ^ oo that is its asymptotic value for 
large time. For a finite Hilbert space size M the fidelity will start to fluctuate for long times 
due to a discreteness of the spectrum of the evolution operator. The size of this fluctuations 
can be calculated by evaluating a time average fidelity F 

1 m 

F= lim — VF(t). (2.34) 

m— »00 777 ^-^ 
t=0 

This is easily done if we expand the initial state in eigenbasis of the unperturbed propagator 
Uq and denote the matrix elements between unperturbed and perturbed eigenstates by Pm 

Uo\^i) = exp{-i4>iMi), Us\<pf)=exp{-i<t>f)\<pf), P„i = {^k\<ph- (2-35) 

We denoted the eigenphases of unperturbed and perturbed one-step propagator with (pi and (/>f , 
respectively. The matrix P is unitary and in the case when both eigenvectors can be chosen to 
be real it is orthogonal. This happens if Uq and Us commute with an antiunitary-f operator T 
whose square is identity §. The fidelity amplitude can now be written 

/(*) = E iP^P)imPmi exp {-[{cpf - cl)^)t), (2.36) 

Im 

with pim = {(f>i\p{0)\(f>m) being the matrix elements of the initial density matrix in the unper- 
turbed eigenbasis. To calculate the average fidelity F we have to take the absolute value square 
of f{t). Averaging over time t we will assume that the phases are nondegenerate 



exp (i((/)f - (pf, +(t)ra- (Pm')t) = ^mrn'^U'. (2.37) 



This results in the average fidelity 

F = Y.\ipP)ml?\Pml?- (2.38) 

ml 

The time averaged fidelity therefore understandably depends on the initial state p as well as on 
the "overlap" matrix P. 



*Antiunitary operator must satisfy (Tip\T(j)) = {iIj\(J))* . 

^Square of an arbitrary antiunitary operator is T^ = ±1. Time reversal operator for a system with spin is 
T — exp {—inSy)K, with a complex conjugation operator K. If the system has an integer spin (or half-integer 
spin and an additional rotational invariance symmetry) we have T'^ — 1. 
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For small perturbation strengths, say 5 smaller than some critical d^m, the unitary matrix P 
will be close to identity. Using P ^ 1 for 6 <^ 6rm (|^-38() gives us 

Ke.k = J2pl (2-39) 

I 

One should keep in mind that for the above result Fweak we needed eigenphases to be nonde- 
generate, (/>f ^ (pm l^'2.'A7^ . and at the same time P — > 1. This approximation is justified in the 
lowest order in 6, when off diagonal matrix elements are |-Pmip oc 6"^. 

On the other hand, for sufficiently large 6 ^ 5rm and complex perturbations V one might 
assume P to be close to a random matrix with independent real or complex matrix elements 
Pmi- Then we can average expression (|2.H8j) over a Gaussian distribution oc exp (— /3AA|Pm/p/2) 
of matrix elements Pmh where we have /3 = 1 for orthogonal P and (3 = 2 for unitary P. This 
averaging gives (|Pmz|^) = (4 — /?)/AA^ and (|-Pm,zP) = 1/AA for the variance oi Pmi (brackets (•) 
denote here averaging over the distribution of matrix elements and not over the initial state). 
The average fidelity for strong perturbation can therefore be expressed as 

I ly^m 

The point of crossover Jj-m from weak (|2.H9|) to strong (|2.4()j) perturbation regime is system 
dependent and can not be discussed in general apart from expecting it to scale with h similarly 
as a mean level spacing 6rm ~ ^ • We will discuss the value of F for three different initial states: 

(i) First, let us consider the simplest case when the initial state is an eigenstate of Uq say, 
P — 1 01 ) (01 1 with matrix elements pim = Si^i5m^i. For weak perturbations this gives 1)2. 39() 
-^weak = 1) therefore the fidelity does not decay at all. This result can be generalised to 
the case when p is a superposition of a number of eigenstates, say K of them, all with 
approximately the same weight, so that one has diagonal density matrix elements of order 
pii ~ l/K, resulting in -F„eak ~ l/K. On the other hand, for strong perturbation 5 ^ (5rm 
we get Fstrong = (4 — /3)/AA for an initial eigenstate. Summarising, for an initial eigenstate 
we have time averaged values of fidelity 



Aveak 



1, i^strong = (4 - (3)/M. (2.41) 



With this simple result we can easily explain the numerical result of Peres ( Pered . [l995|) 



where no-decay of fidelity was found for a coherent initial state sitting in the centre of 
an elliptic island, thus being a superposition of a very small number of eigenstates (it 
is almost an eigenstate). The behaviour in generic case may be drastically different as 
described in the present work. 

Second, consider the case of a random pure initial state \^) = J2m'^m\4'm), giving p^i = 
CmC*i . The coefficients Cm are independent random complex Gaussian variables with vari- 
ance 1/AA, resulting in averages (|/0/m.P) = 1/AA^ ioi m ^ I and {pfi) = 2/AA^ (average is 
over Gaussian distribution of Cm)- Using this in expressions for average fidelity ()2.39|) and 
(|Tin|) we get 

i^weak = 2/AA, Fstrong = l/M. (2.42) 

For random initial state there is therefore only a factor of 2 difference between finite size 
fiuctuating plateau for weak and for strong perturbation. The result for weak perturba- 
tion agrees with the case (i) where we had -Fweak ~ 1/K if there were K participating 
eigenvectors. 
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(iii) Third, for a uniform average over the whole Hilbert space, i.e. taking a non-pure initial 
density matrix p = 1/M, we have 



K 



weak 



l/AA, 



R 



strong 



(4-/3)/AA2 



(2.43) 



As expected, the fluctuating plateau is the smallest for an uniform average over the whole 
Hilbert space and strong perturbation. 
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Figure 2.1: Dependence of time averaged fidelity (multiplied by the Hilbert space size J\f = S) 
F on 5 is shown for a chaotic kicked top system and Hilbert space average p = 1/M, i.e. our 
case (iii). The transition from weak to strong perturbation regime is seen (|2.43() . Horizontal 
full lines are the theoretical predictions Kjtrong (|2.43j) . while the theoretical result for the weak 
regime corresponds to 1. 



Observe that the average fidelity F (|2.38)) is of fourth order in matrix elements of P, the same 
as the inverse participation ratio (IPR) of the perturbed eigenstates. Actually, in the case of 
initial eigenstate, our case (i), the average fidelity (|2.38|) can be rewritten as F = J^m l-fiml"^, 
exactly the IPR. The inverse of the IPR is a number between 1 and M which can be thought of 
as giving the approximate number of unperturbed eigenstates represented in the expansion of 
a given perturbed eigenstate. For an average over the whole space, case (iii), we have instead 



^ — ^im l^iml^/AA^) i-e. the average IPR divided by M. The time averaged fidelity is thus 
directly related to the localisation properties of eigenstates of Us in terms of eigenstates of 
Uq. However, except for the pathological case of the initial state being a small combination of 
eigenstates of C/q with weak perturbation, the fidelity fiuctuation is always between the limiting 
values 2/AA, and 3/AA^. Therefore, the fidelity will decay only until it reaches the value of 
finite size fiuctuations and will fluctuate around F thereafter. The time too when this happens, 
F{too) = F, depends on the decay of fldelity and will be discussed in subsequent chapters. 

To illustrate the above theory we have calculated the average fidelity (|2.38() for a kicked top 
with a propagator (|2.2flj) . As an initial state we used p = 1/M, i.e. the case (iii), where the 
Hilbert space size is determined by the spin value, M = S (OF subspace). We calculated the 
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dependence of -F on 5 for two cases: a chaotic one for kicked top parameters a = 30, 7 = 7r/2 
shown in Figure 12.11 and a regular one for a = 0.1, 7 = 7r/2 shown in Figure 12.21 In both 
cases one can see a transition from the weak perturbation regime Fweak = 1/AA to the strong 
regime -Fstrong = 3/AA^ for large 5. In the chaotic case the critical (5rm can be seen to scale 
as 5rm r^ h = 1/S. In the regular situation, the strong perturbation regime is reached only 
for a strong perturbation (5 ~ 4, where the propagator Us itself becomes chaotic. Namely, the 
transition from the regular to chaotic regime in the kicked top happens at around a = 3, see 
e.g. ( Perea . |l995). Still, if one defines drm as the points where the deviation from the weak 
regime starts (point of departure from 1 in Figure [22)) one has scaling (5rm ~ 1/S also in the 
regular regime. 
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Figure 2.2: The same as Figure l?!T] but for a regular kicked top. 



2.2.2 State Averaged Fidelity 

In previous section we used the initial density matrix p(0) = 1/M in expression (|2.9|) for the 
fidelity amplitude to calculate the average fidelity over the whole Hilbert space. As the fidelity is 
of fourth order in the initial state |V'(0)), whereas the fidelity amplitude is bilinear in |V'(0)) the 
average fidelity is not equal to the quantity obtained by first averaging fidelity amplitude and 
subsequently squaring it. In general the difference between the two depends on the ensemble of 
states over which we average. 

Let us look only at the simplest case of averaging over random initial states, denoted by 
((•)). In the asymptotic limit of large Hilbert space size M ^ 00 the averaging is simplified by 
the fact that the expansion coefficients Cm of a random initial state in an arbitrary basis become 
independent Gaussian variables with variance 1/AA. Quantities bilinear in the initial state, like 
the fidelity amplitude or the correlation function, result in the following expression 



UlAm =: {{A)) = {{cl^A^ici)) = ^t^A, 



(2.44) 



where A is an arbitrary operator. The averaging is done simply by means of a trace over the 
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whole Hilbert space. For the fidehty F{t) which is of fourth order in 1-0) we get 

{{F{t))) = {{cUMs{t)UciCp[Ms{t)];,.c:.)) = \{{f{t)))f + l. (2.45) 

The difference be tween the average fid ehty and the average fidelity amplitude is therefore semi- 



classically small (|Prosen et al. Ll2003a ) 



There are two reasons why averaging over random initial states is of interest. First, in the 
field of quantum information processing this are the most interesting states as they have the 
least structure, i.e. can accomodate the largest amount of information. Second, for ergodic 
dynamics and sufficiently long times one can replace expectation values in a specific generic 
state 10) by an ergodic average. If the system is ergodic on the whole space one can calculate 
the correlation function C{j,k) (|2.18j) by means of a simple trace (|2.44)) . so that it does not 
depend on the initial state, or even calculate it classically in the leading semiclassical order, 
which greatly simplifies theoretical derivations ". For regular systems, where the decay of 
fidelity depends on the initial state, the ergodic averaging differs from the average in a specific 
initial state, although one can still be interested in the behaviour of the average fidelity. Such 
an averaging will be discussed in the section describing the decay of fidelity in regular systems. 

We have seen that for sufficiently large Hilbert spaces there is no difference between averaging 
the fidelity amplitude or the fidelity or taking a single random initial state. For mixing dynamics 
the long time fidelity decay is independent of the initial state even if it is a non-random, 
wheras in the regular regime it is state dependent. For instance, the long time Gaussian decay 
f Section l3. 2. 2|) depends on the position of the initial coherent state. The fidelity averaged over 
this position of the initial coherent state might be of interest and will not be equal to the fidelity 
averaged over random initial states. We will discuss averaging over coherent states at the end 



of Section 1222] describing long time fidelity decay for coherent initial states. 



^In the case of autonomous systems the canonical or micro-canonical averaging should be used instead. 



Chapter 3 

General Perturbation 



An expert is a person who has made all the 
mistakes that can be made in a very narrow 
field. 

— Niels Bohr 



3.1 Mixing Dynamics 

Here we assume that the system is mixing such that the correlation function of the perturbation 
V decays sufficiently fast; this typically (but not necessarily) corresponds to globally chaotic 
classical motion. Due to ergodicity we will assume the initial density matrix to be p{0) = 1/M, 
so that all averages over a specific initial state can be replaced by a full Hilbert space average, 
(.) = tr(.)/AA (Section 12. 2. 2j) . For any other initial state (e.g. in the worst case for the minimal 
wave packet - coherent state) one obtains identical results on F(t) for sufficiently long times*, 
i.e. longer than the Ehrenfest time tE ~ ln(l//i)/A (for a classically chaotic system with maximal 
Lyapunov ex ponent A) needed for a minimal wave packet to spread effectively over the accessible 
phase space ( Berman fc Zaslavskvl . Il978l ). The state averaged quantum correlation function is 



homogeneous in time, i.e. C{j,k) = C{k — j), so we simplify the second order linear response 
formula for the fidelity (|2.17j) 

Fit) = 1-^1 tCiO) + 2 ^ (t - j)CU) I + 0(5'). (3.1) 

If the decay of the correlation function C{j) is sufficiently fast, namely if its integral con- 
verges on a certain characteristic mixing time scale tmixj then the above formula can be further 
simplified. For times t ^ tmix we can neglect the second term under the summation in (|3.1() 
and obtain a linear fidelity decay in time t (in the linear response) 

F{t) = 1 - 2{6/hfat, (3.2) 

with the transport coefficient a being 



*The exception might be systems with locahzed states. 

23 



24 



Chapter 3. General Perturbation 



In the continuous time case a is just the integral of the correlation function. Note that a 
has a well defined classical limit obtained from the classical correlation function and in the 
semiclassical limit this classical Uci will agree with the quantum one. 

We can make a stronger statement in a non-linear-response regime if we make an additional 
assumption on the factorisation of higher order time-correlations, n— point mixing. This implies 
that 2rn--point correlation iy{ji) ■ ■ ■ V{J2m)) is appreciably different from zero for j2m — Ji — ^ c>o 
only if all (ordered) time indices {jk ,k = 1 . . . 2m} are paired with the time differences within 
each pair, J2k — J2fc-i > being of the order or less than tmix- Then we can make a further reduction, 
namely if t S> mt^[^ the terms in the expansion of the fidelity amplitude f{t) (|2.16|) are 
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Figure 3.1: The classical correlation function 1)2. 18() of perturbation V 1)3. 6|) for chaotic kicked 
top and 7 = 7r/6 (top solid curve) and 7 = 7r/2 (bottom broken curve). The finite time 
integrated correlation function is shown in the inset. 
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The fidelity amplitude is therefore /(t) = exp {—6'^at/h ) and the fidelity is 



Fit) = exp (-t/r„ 



h' 



2(^Vd' 



(3.5) 



with a mixing decay time-scale t^ = 0{6~'^) and a classical limit of the transport coefficient 
cJci. We should stress again that the above result (|3.5j) has been derived under the assumption 
of true quantum mixing which can be justified only in the limit M — > 00, e.g. either in the 
semiclassical or the thermodynamic limit. Thus for the true quantum-mixing dynamics the 
fidelity will decay exponentially. The same result has been derived also by a quite different 
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approach, using a Fermi golden rule (|.Tacquod et al. 1 . 120011 : ICerruti fc Tomsovid 12002 ^ . That 
is why this regime of exponential fidelity decay is sometimes called a Fermi golden rule regime. 
To numerically check the above exponential decay, we will use the kicked top (|2.20|) with 
parameter a = 30, giving a totally chaotic classical dynamics. As argued before, one can 
calculate the transport coefficient a (|3.3|) by using the classical correlation function of the 
perturbation (|T^ . 



Vr 



cl 



(3.6) 



We consider two different values of kicked top parameter 7, namely 7 = 7r/2 and 7 = 7r/6. The 
classical correlation functions can be seen in Figure IXTl The correlation function (obtained by 
averaging over 10^ initial conditions on a sphere) is shown in the main frame. The correlation 
functions have qualitatively different decay towards zero for the two chosen 7's. In the inset the 
convergence of classical a (|3.3j) is shown, where one can see that the mixing time is imix ~ 5. 
The values of a^ are a^ = 0.00385 for 7 = 7r/2 and a^i = 0.0515 for 7 = vr/G. These values 
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Figure 3.2: Quantum fidelity decay in the chaotic regime for 7 = 7r/2 and three different 
perturbation strengths 6 = 5 x 10~^, 1 x 10"^ and 3 x 10"'^ (solid, dashed and dotted curves, 
respectively) is shown. The chain line gives theoretical decay 1)3. 5|) with the classically calculated 
a seen in Figure IXTl 

are used to calculate the theoretical decay of fidelity F{t) = exp {—6'^S'^2acit) (|3.5j) which is 
compared with the numerical simulation in Figures 13.21 and 13.31 We used averaging over the 
whole Hilbert space, p(0) = 1/5 and checked that that due to ergodicity there was no difference 
for large S if we choose a fixed initial state, say a coherent state. As fidelity will decay only 
until it reaches its finite size fluctuating value F (|2.43|) we choose a large S = 4000 in order to 
be able to check exponential decay over as many orders of magnitude as possible. In Figure 15?^ 
the decay of quantum fidelity is shown for 7 = tt/2. The agreement with the theory is excellent. 
Note that the largest 6 shown corresponds to Tm ~ 1 so the condition for n— point mixing 
t ^ tmix is no longer satisfied. The agreement with theory is still good due to the oscillatory 
nature of the correlation decay (see Figure 13. 1|) fullfiling the factorisation assumption 1)3.4(1 on 
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average. In Figure HOI for 7 = 7r/6 a similar decay can be seen. In both cases fidelity starts to 
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Figure 3.3: Similar figure as 13. 2( only for 7 = 7r/6 and perturbation strengths 5 = 1 x 10 ^, 
2 X 10~^ and 3 x 10~^ (solid, dashed and dotted curves, respectively). 

fluctuate around F calculated in the section r2.2. II for times larger than too- 

3.1.1 Long Time Behaviour 

We assumed that the quantum correlation function C{j) decays to zero and its integral converges 
to a. For a system having a finite Hilbert space of size M, the correlation function asymptotically 
does not decay but has a non-vanishing plateau C due to finite M, similarly as we have a 
finite asymptotic value of the fidelity F. This will cause the double correlation sum to grow 
quadratically with time. Because this plateau C is small, the quadratic growth will overtake 
linear growth 2a"cit only for large times. The time averaged correlation function C{j, k) 1)2. 18() 
can be calculated assuming a nondegenerate unperturbed spectrum cp^ as 



1 



t-i 



^ = }^^72 E CU,k)=Y.PkkiV, 



t2 



kk) 



j,k=0 



/ , Pkk^kk 
\ k J 



(3.7) 



where pky^ are diagonal matrix elements of the initial density matrix /9(0) and V^fc are diagonal 
matrix elements of the perturbation V in the eigenbasis of the unperturbed propagator Uq. 
One can see that C depends only on the diagonal matrix elements ' , in fact it is equal to the 
variance of the diagonal matrix elements. Since the classical system is ergod i c and mixing, we 
will use a version of t he quantum chaos conjecture (IF eingo ld Sz Pered . Il98q : IWilkinsonl . Il987 ; 



Feingold et al I. Il98fll: JProsen fc R.obnikL Il99,4 IProsenL IjJHJ) saying that Km are independent 
Gaussian random variables with a variance given by the Fourier transformation S{uj) (divided 
by M) of the corresponding classical correlation function Cci(j) at frequency to = 4>m — 4'n- On 



^The case when diagonal matrix elements are zero is the subject of Chapter |3 
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Figure 3.4: The finite time quantum correlation sum (7(t) = Y^- ^^q C(j, fc)/2i (solid curves) 
together with the corresponding classical sum a"ci(^) = J2i~k=oCci{j,k)/2t (dashed curves satu- 
rating at (Tci and ending at t ~ 1000) is shown for the chaotic kicked top. Quantum data are 
for a full trace p = l/M with S = 1500. Upper curves are for 7 = vr/G while lower curves 
are for 7 = 7r/2. Chain lines are best fits for asymptotic linear functions corresponding to 
Ct/2 = 0.0077t/5 for 7 = 7r/2 and 0.103t/S for 7 = vr/G. 



the diagonal we have a; = and an additional factor of 2 due to random matrix measure on the 
diagonal . Using 2acit = J2\ k C(i) k) = S{{))t we can write 



C 



25(0) 4aei 



M 



M ■ 



(3.8) 



Because of ergodicity, for large A^, C does not depend on the statistical operator p used in 
the definition of the correlation function, provided we do not take some non generic state like 
a single eigenstate \(j)k) for instance. Note that Equation (|3.8|) is valid on a single quantum 
invariant subspace. If Uq has symmetries, so that its Hilbert space is split into s components of 
sizes Mj, the average C will be different on different subspaces, Cj = Aad/Mj. Averaging over 
all invariant subspaces then gives 

SO that C is increased by a factor s compared to the situation with only a single subspace. 
The fidelity decay will start to be dominated by the average plateau 1)3. 8() at time in when the 
quadratic growth takes over, Ct\ « 2(TcitHi 



tn = -M^ h^'^, 



(3.10) 



which is nothing but the Heisenberg time associated with the inverse density of states. Again, 
if one has s invar iant subspaces, the Heisenb erg time is tu = M /2s. This crossove r time agrees 
with t he result of lCerruti fc Tomsovid ((2Qi]2) and for random matrix models with lGorin et al 
( 2004h . 
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In Figure l,S.4l we show numerical calculation of the correlation sum for the chaotic kicked 
top at a = 30. We compare the classical correlation sum (the same data as in Figure IXT]) and 
quantum correlation sum. One can nicely see the crossover from linear growth of quantum cor- 
relation sum 2(Tcit for small times t < tu, to the asymptotic quadratic growth due to correlation 
plateau C. In addition, numerically fitted asymptotic growth 0.103t/S and 0.0077^/5" nicely 
agree with formula for C, using AA = S and classical values of transport coefficients ad = 0.0515 
and 0.00385 for 7 = vr/G and 7 = vr/2, respectively. 
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Figure 3.5: Histogram of the normalised distribution of the diagonal matrix elements Vkk for 
the chaotic kicked top and S = 4000 on OE subspace (|2.24|) . The dotted line is the theoretical 
Gaussian distribution with the second moment C and the two chain lines are expected \/7^ 
statistical deviations if there are Ni elements in the i—th bin. Note the different x-ranges in 
two figures due to different ad for the two chosen 7. 

For times t > ^H; and provided 5 is sufficiently small, the correlation sum will grow quadrat- 
ically and the linear response fidelity reads 



Fit) 



^4aci 2 



(3.11) 



To derive the decay of fidelity beyond the linear response regime one needs higher order moments 
of diagonal elements of perturbation V. If we use the BCH form of the echo operator 1)2.14(1 and 
discard the term involving ^ T{t), we have the fidelity amplitude f{t) = X)fc ^xp {—iVkkSt/h)/M, 
where we choose an ergodic average p = 1/M. In the limit AA ^ 00 we can replace the sum 
with an integral over the probability distribution of diagonal matrix elements p{Vkk) = p(^)) 



f{t)= fdVp{V)exp{-iV6t/h). 



(3.12) 



For long times the fidelity amplitude is therefore a Fourier transformation of the distribution 
of diagonal matrix elements. For classically mixing systems the distribution is conjectured to 
be Gaussian with the second moment equal to C (|3.8|) . This is confirmed by numerical data in 
Figure 13.51 The mean value of diagonal matrix elements is perturbation specific and is for our 

{2S + 1){S + 1)/12S^. 



choice of the perturbation (j2:22jl J^k Vkk/{2S + 1) = (25 + 1)(S' + 1)/125^. From the figure we 
can see that the distribution is indeed Gaussian with the variance agreeing with the theoretically 
predicted C = 4(Tci/5'. The Fourier transformation of a Gaussian is readily calculated and we 



■'"in Chapter|l|we will see that the size of r(f) term grows at most linearly with time and so can be neglected 
because of S^ prefactor. 
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get a Gaussian fidelity decay 

F{t) 



exp 



(-(V^p) 



A^ h 



(3.13) 



^ V 4^ci 5 

In order to see a Gaussian fidelity decay for mixing systems the perturbation strength 6 must 
be sufficiently small. If it is not, the fidelity will decay exponentially (|3.5j) to its fluctuating 
plateau F before time tn when the Gaussian decay starts. Demanding that the mixing decay 
time Tm is smaller than tjj = AA/2 gives the critical perturbation strength <5p, 

For 5 < 5p we will have a Gaussian decay 1)3. 13() otherwise the decay will be exponential (|3.5I) . 
for details see Section rUH 
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Figure 3.6: Quantum fidelity decay for 6 < 6p in the chaotic regime. For 7 = 7r/2 data for 
6 = 1- 10^^ (solid curve) and 5 • 10^^ (dotted curve) are shown. For 7 = tt/6, 6 = 3- 10^^ 
(solid) and 1 • 10^^ (dotted) are shown. Note that for both 7 the curves for both 6 practically 
overlap. The chain curves are theoretical predictions (|3.13|) with classically computed ad- 

Again we numerically checked the predicted Gaussian decay for the chaotic kicked top with 
a = 30 and S = 1500 and a full trace average over Hilbert space. The results of numerical 
simulation, together with the theory are shown in Figure 15^ 



T he regime of Gaussian decay is sometimes referred to as the perturbative regime ()Cerruti &: Tomsovic , 

2002l : I.Tacquod et al. I . I2OOII ) because it can be derived using the lowest order perturbation theory. 
Writing first order corrections in phases (pf. = (pk+^kk^/^ ^-nd for overlap matrix Pki = 6ki+0{6) 
(|2.36|) one gets the Fourier transform formula (|3.12j) . 



3.2 Regular Dynamics 

As opposed to mixing dynamics which was characterised by a linear growth of the double sum 
of the correlation function, the regular regime will typically exhibit quadratic growth. By a 
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regular regime we again refer to the behaviour of the correlation function. Typically, its double 
sum will exhibit quadratic growth for integrable or mixed (KAM) systems. In such case, we 
can define a time average correlation function, denoted by C 

1 *-^ 

For nonzero C we will call this a "regular regime" . Of course, we have AA-dependent nonzero 

value of C also in mixing systems ()3.8() , as discussed in Section 13.1.11 Here we consider only 

systems where C exists in the limit of AA ^ oo, i.e. nonzero C is a consequence of the dynamics 

and not of the finite Hilbert space size. Due to non-ergodicity all expectation values (like C) 

depend on the initial state, and we can not use ergodic average over the whole Hilbert space as 

for mixing dynamics. The time in which time averaging of the correlation function converges 

will be denoted by iave- For times t ^ t^ve the linear response fidelity decays quadratically in 

time 

x2 _ 

F{t) = 1 - -^Ct^. (3.16) 

U 

In contrast to the mixing regime one does not need any further assumptions in order to go 

beyond the linear response formula. Namely, for times t ^ tave we can define a time average 

perturbation V 

y = lim - y V{k) = lim ^. (3.17) 

k=0 

Observe that V is by construction a constant of motion, [V, Uq] = 0. In a mixing regime V is 
trivial (in the limit M -^ oo), i.e. proportional to the identity ^, whereas for regular dynamics 
nontrivial V exists. The special case of y = will be considered in Chapter HI In the case 
of nondegenerate spectrum of the unperturbed propagator Uq, the time average is simply the 
diagonal part of V in the eigenbasis \(pk) of the unperturbed propagator, namely 



v = y2Vkkm{A\, (3.18) 



where as usual V^k = {4'k\^\4'k)- Note that the average correlation function is 

C = (y2) _ (y)2. (3.19) 

For times t ^ tave the operator S(t) is dominated by the linearly growing term Vt and one can 
neglect contributions not growing with time. The BCH form of the echo operator can therefore 
be simply written as 

Ms{t) = exp {-iV6t/n), t > tave, (3.20) 

where we neglected the term ~ T{t)6'^ as it becomes important only at large times ~ 1/5^, 
whereas the first term will cause the fidelity to decay already in time ~ 1/5. The above form 
of the echo operator will be the main ingredient of theoretical calculation of the fidelity decay 
in the regular regime. 

As the operator V commutes with the unperturbed propagator it is diagonal in the eigenba- 
sis I^A;)- In integrable systems the basis states {(pk) can be ordered in a very special way. Namely, 
there exist quantum numbers which are eigenvalues of canonical action operators having a very 
simple algebra. Using action-angle operators will make derivations easier and will furthermore 
enable us to use classical action-angle variables in the leading semiclassical order, thereby ap- 
proximating quantum fidelity in classical terms. So before proceeding with the evaluation of 
f{t) for various initial states, let us have a look at action-angle operators. 

^Actually it can be made zero by subtracting from V the identity operator which only rotates a phase of the 
fidehty amplitude and does not affect the fidelity itself. 
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3.2.1 Action-angle Operators 

Since we assume the classical system to be completely integrable (at least locally, by KAM 
theorem, in the phase space part of interest) we can employ action-angle variables, {jk, 6k,k = 
1 . . .d}, in d degrees of freedom system. In the present section, dealing with the regular regime 
as well as in Section [4.2.11 describing quantum freeze of fidelity in regular systems, we shall use 
lowercase letters to denote classical variables and capital letters to denote the corresponding 
quantum operators. For instance, the quantum Hamiltonian will be given as H{J, 0) whereas 
its classical limit will be written as h{j,0). 

As our unperturbed Hamiltonian is integrable, it is a function of actions only, i.e. ho = ho{j). 
The solution of classical equations of motion is very simple, 

Jit) = j, 

e{t) = e + u{3)t (mod 2tt) (3.21) 

with a dimensionless frequency vector 

-0) - ^. (3.22) 

The classical limit v{j,0) of our perturbation generator V can be written as a Fourier series 

Vij,e)= Y. ^rr.UV'^-', (3.23) 

where a multi-index m has d components. The classical limit of the time-averaged perturbation 
V \s V = voij), i.e. just the zeroth Fourier mode of the perturbation. 

In quantum me chanics, one quantises the action-angle variables using the famous EBK 
procedure (see e.g. (BerriJ, |l977|)) where one defines the action (momentum) operators J and 



angle operators exp(im • 0) satisfying the canonical commutation relations, 

[JA;,exp(im • 0)] = ?imfcexp(iTn • 0), /c = l,...,d. (3.24) 

As the action operators are mutually commuting they have a common eigenbasis \n) labelled 
by d-tuple of quantum numbers n = (ni, . . . , n^), 

J\n) = h{n + a.)\n) (3.25) 

where < a^ < 1 are the Maslov indices which are irrelevant for the leading order seiniclassical 
approximation we will use. It follows from l|3.24|) that the angle operators act as shifts 

exp(im, • 0)|n) = |n -|- m). (3.26) 

The Heisenberg equations of motion can be trivially solved in the leading semiclassical order by 
simply disregarding the operator ordering, 

J{t) = e^Hot/hj^-iHot/h ^ J 

in terms of the frequency operator ij^{J). Throughout this paper we use the symbol = for 
'semiclassically equal', i.e. asymptotically equal in the leading order in U. Similarly, time 
evolution of the perturbation observable is obtained in the leading order by substitution of 
classical with quantal action-angle variables in the expression (|3.23)) 

V{t) = e^iiot/hy^-iHot/h ^ J2 1;^( J)e^"^-"('^)*e'^-®. (3.28) 



m 
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Operator V{J) is diagonal in the eigenbasis \n) and therefore the expectation value of the 
echo operator ()3.2U|) in the initial density matrix p is 

f{t) = J2 exp {-i6tV{h{n + a})/h)Dp{hn), Dp{hn) = {n\p\n). (3.29) 

n 

For pure initial states Dp is just Dp = \{'>p\n)\'^. This is still the exact quantum mechanical 
expression of the fidelity amplitude. Now we make a leading order semiclassical approximation 
by replacing quantum V with its classical limit v and replacing the sum over quantum numbers 
n with the integral over classical actions j . The replacement of the sum with the action space 
integral (ASI) is valid up to classically long times ta_, such that the variation of the argument 
in the exponential across one Planck cell is small, 

ia = Tw^ - n'>/5. (3.30) 

\Ojv\d 

Subsequently we will see that the fidelity decays on shorter times and so the approximation 
is justifiable. By denoting with dp{j) the classical limit of Dp{hn) we arrive at the fidelity 
amplitude 

fit) ^ n-'^Jd^jexp (-i^t^(i)) dpij). (3.31) 

This expression will be used to evaluate the fidelity for different initial states. 

3.2.2 Coherent Initial States 

We proceed to evaluate fidelity amplitude (|3.31j) for coherent initial state centred at {j*,6*) 
in action-angle phase space. The expansion coefficients of the initial density matrix p = 
\j*, 0*){j* , 9*\ for a general coherent state in d degrees of freedom system can be written as 

{n\j*,e*)=(^^^ |detA|^/Sxp|-^(J„-j*).A(J„-r)-in.0*|, (3.32) 

where A is a positive symmetric dx d matrix of squeezing parameters and Jn = fi{n + a) is an 
eigenvalue of operator J in eigenstate |n). The classical limit of Dp (|3.29j) is therefore 

dpU) = {ni-nY'^ \detA\^/\xp{-{j - f) ■ Aij - f)/h). (3.33) 

The ASI for the fidelity amplitude can now be evaluated by the stationary phase method. To 
lowest order in 6t the stationary point js, i.e. the zero of [—i5tv — (jf — j*) ■ A{j — j*)], is at 

j,=r + ^A-V + 0(<52), (3.34) 

where the derivative of the average perturbation is 

C' = ^. (3.35) 

dj 

For 5t <^ 1 the stationary point is simply at j* and the fidelity is 



F{t) = exp {-(t/rr)^}, Tr = ]^^r^- (3.36) 

The fidelity amplitude on the other hand has an additional phase, 

f{t) = exp {-(t/Tr)V2} exp {-iv{j*)6t/h). (3.37) 
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Comparing the fidelity (|l-i.H(ij) with the linear response formula (|IH6|) we see that the average 
correlation function for a coherent initial state is (7 = 2^(^' ■ A~^i;'). 

We hav e already briefl y explained the behaviour of classical fldelity in Section 11.1.21 for 
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details see (JBenenti et al. L 12003a'). For regular dynamics one gets a ballistic decay of classical 



fidelity for perturbations predominantly changing frequencies of tori or power-law decay for 
perturbations predominantly changing the shape of tori. In the ballistic case, provided one is 
allowed to replace the perturbation v with its time average v{j), one gets the same Gaussian 
decay of the classical fidelity as for the quantum one ()3.36|1 . Therefore, in the regime where 
V determines the decay of classical fidelity, classical and quantum fidelity agree. From the ex- 
pression for the echo operator (|3.2()j) one can in this case interpret the fidelity as the overlap 
between the initial coherent state and the state obtained after the evolution with the Hamil- 
tonian 5v. Because v is a function of actions only, the classical equations of motion are very 
simple, namely only the frequency of classical motion is changed by the amount Acj = v' 5. 
This change in frequency causes the "echo" packet Ms{t)\'il:{Q)) to move ballistically away from 
its initial position and as a consequence fidelity decays. For coherent state the shape of the 
packet is Gaussian and therefore fidelity decay is also Gaussian. For other forms of localized 
initial packets the functional form of the fidelity will be different but with the same dependence 
on the ballistic separation "speed" v' . If y = we obviously have v' = 0, this situation will be 
discussed in Chapter 0] describing a freeze of fidelity. 

In one dimensional systems {d = 1) another phenomena should be observable. After long 
time the echo packet will make a whole revolution around the torus causing the fidelity to be 
large again. This will happen after the so-called beating time tb determined by 6v'{j*)t\^ = 2tt, 

t, = -. (3.38) 

This beating phenomena is particular to one dimensional systems as in general the incomensu- 
rability of frequencies will suppress the revivals of fidelity. 

Let us now illustrate the above theory with numerical experiment. We again take the kicked 
top ()2.20|) . but this time with 7 = tt/2 and a = 0.1 resulting in regular classical dynamics. The 
perturbation is the same as for mixing situation, i.e. a perturbation in a resulting in a classical 
perturbation v = z^ /2 1)2.22(1 . Let us denote by ??, (p the spherical angular coordinates measured 
with respect to the y-axis. For a = the action-angle variables are {j = cos"!? = y,9 = ip) and 
we can use the preceding theory. The squeezing parameter for spin coherent states is (|2.26|) 

A = l/sin2i? = 1/(1-2/2). (3.39) 

To calculate the decay time t^ we need the derivative of the average perturbation. For a = 
the evolution is very simple, just a rotation by 7r/2 around the y-axis, resulting in 

v = \{l-f), v' = -^-. (3.40) 

Theoretical decay of fidelity is therefore 

Fit) = exp [-6''StV{l - y')/8) , (3.41) 

where we used h = 1/5. In Figure 13.71 we show the results of numerical calculation for a co- 
herent packet placed at spherical coordinates (with respect to z-axis) {■&*, (p*) = 7r(l/\/3, ^/V2) 
resulting in action y = 0.77 (note that we projected the initial state onto the invariant OE 
subspace ()2.24j) ). Spin number is chosen S = 100. Excellent agreement with the theoretical 
Gaussian decay can be seen. In addition, we can observe a revival of the fidelity at tb- For our 
OE subspace the initial state is composed of two symmetrically positioned images, so that the 
beating time is tb = 2Tr/y5, which nicely agrees with the numerics. 
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Figure 3.7: Quantum fidelity for a coherent state of the regular kicked top, a 
7 = 7r/2 and 5 = 0.0025. The full curve is the theoretical decay 
indicates the theoretical beating time 1)3. 38|) . 



0.1, S = 100, 
and the vertical line is 



Average Fidelity 



■Jd''jF{t,j). 



(3.42) 



If F{t,j) is Gaussian fidelity decay (|3.36|1 for a coherent initial state centred at j, we can define 
the average fidelity 

Packing all j dependent terms in a non- negative scalar function g{j), the fidelity decay for a 
single coherent state can be written as F{t,j) := exp {—S^t'^g{j)/h). For large 6'^t'^/h the main 
contribution to the average will come from regions of small g{j) where the fidelity decay is slow. 
In general the function g{j) will have zeros in action space, for simplicity let us assume there is 
a single zero at j* of order rj. The asymptotic decay can then be calculated and scales as 






(3.43) 



where the sign x will denote "in the asymptotic limit" troughout the work. For infinite phase 
space the order of a zero r] can only be an even number, whereas for a finite space rj can also 
be odd. 

To illustrate the theory, we calculated the average decay for the kicked top and the per- 
turbation in a used before. The function g is g{j) = j'^{l — j^)/8. We have three zeros, two 
of order ry = 1 at the boundary of the phase space j = ±1 and one of order 77 = 2 at j = 0. 
Asymptotically the zero with r] = 2 will dominate, giving the decay 



(Fit)), 



27r 



(3.44) 



'' 6t^/S' 

Figure 13.81 shows that we indeed get asymptotic ~ t~^ decay. If one averages the fidelity 
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Figure 3.8: Average fidelity decay {F{t))- 1)3. 42() for a regular kicked top with a = 0.1, 7 = tt/2 
and perturbation of strength 6 = 10~^ in the parameter a. The averaging is done over 200 
coherent initial states randomly placed on a sphere and projected onto OE subspace (|2.24|) . 
The Dotted line is the theoretical power law decay ()3.44|) . Pluses show the absolute value 
square of the average fidelity amplitude. 

amplitude, the result is close to zero (pluses in the figure) because the phase present in the 
fidelity amplitude (|3.37|) averages to zero^. 

As a second example, we take the same regular kicked top system and parameters, only the 
perturbation is now in rotation angle 7, i.e. Us := U{'y + 6,a). The function g is in this case 
g{j) = (1 — j^)/2. As opposed to the previous case, now we have only two zeros of order rj = 1 
at j = ±1. The asymptotic decay is therefore 



int))^, 



1 



SH'^S' 



(3.45) 



Denominator is now quadratic in time, whereas for the perturbation S^ it was linear. The 
theoretical decay law is again confirmed by the numerical experiment shown in Figure IT^ The 
asymptotic power law decay of average fidelity is therefore system dependent. 

An entirely analogous method will be used in the next section to calculate the asymptotic 
decay for a random initial state. The only difference is that for a random initial state the zeros 
of V are the relevant quantity, whereas for an average fidelity decay of a coherent state the zeros 
of Tr, i.e. zeros of v' are important. Therefore, in a regular system the asymptotic decay for a 
random initial state is usually not the same as the average decay for a coherent state. 



3.2.3 Random Initial State 

In the previous section we have seen that for localized packets the fidelity decay is dictated 
by the separation of two packets due to the perturbation. The opposite possible choice of 

^One would have to have zeros of v and v' at the same position j* to get a nonzero average, which is generally 
not the case. 
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Figure 3.9: Average fidelity decay {F{t))- 1)3. 42() for a regular kicked top. Everything is the 
same as in Figure 13.81 apart from the perturbation being in the angle 7. Dotted line is the 
theoretical power law decay 1)3. 45() without any fitting parameter. Note that we have different 
power law decay as in Figure 13.81 The vertical dotted line shows the position of the beating 
time tb = T^/S (|3.38|) which is in this case independent of j. If tb depends on j (as in Figure IrTH]) 
such recovery is absent in the average fidelity. 



initial condition is to take completely delocalised state, a uniform mixture of all states p = 1/M 
being the extreme case. For uniform initial density matrix p = 1/M, the Dp takes the form 
Dp{J) = 1/M. The classical limit of Dp is obtained by calculating the number of states M 
using the Thomas-Fermi rule, resulting in the classical limit dp{j) = (27r?i) /V, where V is the 
size of the classical phase space. While for small times fidelity will decay quadratically ()3.16|1 . 
for large 6t/h the integral in f{t) (|3.31|) can again be calculated using the method of stationary 
phase. In contrast to a coherent state case, now one can have more than one stationary point. 
If we have p stationary points, Jn,r] = 1,. . . ,p, where the phase is stationary, dv{jrj)/dj = 0, 
the integral results in 



fit) 



{2^fd/2 



V 



'^/2 P exp{-iit;(j,)5//i 

r,=l 



\Ur,] 



detV^|i/2 



(3.46) 



where {V^jfc^ := d'^v{jn)/djkdji is a matrix of second derivatives at the stationary point 77, 
and i/fj = 7r(m+ — ?7i_)/4 where m± are the numbers of positive/negative eigenvalues of the 
matrix V^. We also assumed that the phase space is infinite. In a finite phase space we will 
have diffractive oscillatory corrections in the stationary phase formula, see end of Section 14.2.41 
or numerical results below. Note also that for a large Hilbert space dimension M the fidelity 
for a single initial random state will also decay according to the asymptotic formula ()3.46|) . 
The decay time of the quantum fidelity for random initial states scales therefore as t,- ~ h/6 
and is by a factor \/h shorter than for coherent initial states. We repeat that the stationary 
phase formula (|3.46|) is expected to be correct in the range const h/6 < t < const' /(5 so it will 
give asymptotic decay of fidelity. Most interesting to note is the power-law dependence on time 
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and perturbation strength, F{t) x [h/{t6)] . With increasing dimensionality (i of a system the 
decay gets faster. This allows for a possible crossover to a Gaussian decay wh en appr oaching 
the thermodynamic limit d — > oo, observed in a class of kicked spin chains bv IProsenI (|2fln2 ). 
Agreement beyond linear response with the Gaussian decay is frequently observ ed also for a 
finite d, e.g. in a spin model of quantum computation ( Prosen fc Znidar id . 120011 ) or even in a 
one dimensional kicked top seen in Figure IB. 101 
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Figure 3.10: Fidelity decay for the regular kicked top, a = 0.1, 7 = 7r/2, 5 = 0.01, 5 = 100. 
The solid curve gives the results of a numerical simulation for p = l/M on the OE subspace 
(|2.24|) . Isolated dots denote the difference between numerical calculation and analytic formula 
H3.48() for a ^ 0, |-Fnum. (0 ~ -^anaii.(*)l- The dotted line gives a predicted asymptotic decay 
(X t~^, and the dashed/chain curves are the predicted fidelity decays at small times, namely 
the second order expansion F{t) = 1 — (5^(5)^/180, and 'improved' one by an exponential. The 
wiggling dashed curve represents the numerics for a single random state. 



For numerical demonstration we take the same kicked top system as for the case of coherent 
initial states, a = 0.1, 7 = tt/2 and S = 100 and perturbation in parameter a, giving v = 
(1 — j^)/4. Rather than comparing the numerical results with the asymptotic formula (|3.46j) we 
will try to directly calculate the fidelity decay in order to demonstrate oscillatory corrections 
due to a finite Hilbert space. We start with the quantum expression for the fidelity amplitude 
1)3. 29(1 . using the eigenvalues {2m — 1)^ of operator S'^ on OE subspace. 



F{t) 






exp (-i(5t(2m - lf/4.S) 



For large S we can replace the sum with an integral and get 



vr 



''<'> = Wt 



erfi( 



ipW4, 



(3.47) 



(3.48) 



where erfi(2;) = ^-1= Jq^ e * dt is a complex error function with the limit lima;^oo |erfi(2e^'^'^y^)| = 
1 to which it approaches by oscillating around 1. We therefore have an analytic expression for 
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the fideUty (|l-i.48j) in the case of an uniform average over the whole Hilbert space or, equivalently, 
over one random initial state. Its asymptotic decay is proportional to t~^ which agrees with 
the general semiclassical asymptotic (j3.46|) but in addition we have an oscillatory erfi correction 
due to a finite space. For small times we can get initial quadratic linear response decay by 
expanding the full theoretical formula (|.S.48j) or more easily by simply calculating C according 
to formula (|3.19|) . One can even use the classical calculation. For a = we have v = {1 — y^)/4, 
so the classical average correlation functions is 

0^ = — (7-7^) = ^, (3.49) 

where the averages are performed over a uniform distribution on a sphere. In Figure 15. lUI we 
compare theoretical decay (|;-{.48j) with the numerical as well as with the linear response formula 
F{t) = 1 — {6S)'^Ccit'^- Also, beyond the linear response there is a significant agreement with 
the Gaussian approximation F{t) = exp (— J^S'^C'dt^) obtained by exponentiating the linear 
response expression. 

3.3 Time Scales 

Let us close the present chapter with an overview of time scales of different fidelity decays 
depending on parameters 6 and h and on the initial state. 

For regular dynamics we have only three relevant time scales, the classical averaging time 
tave in which the average perturbation operator converges to V, the quantum fidelity decay time 
Tr and the time too when the fidelity reaches a fluctuating plateau due to a finite dimension of 
Hilbert space. For times smaller than tave decay is system and state specific and can not be 
discussed in general. After tave the fidelity decay is first quadratic in time as dictated by the 
linear response formula (|3.16j) . The decay time t^ scales as ~ yh/6 (|3.36|1 for coherent initial 
states and as ~ h/5 ()3.46|) for random initial states. Beyond linear response the functional 
dependence of the decay is Gaussian for coherent initial states and power law for random initial 
states with the before mentioned decay time t^. This decay persists until the finite size plateau 
F (|2.38|) is reached at time too. The time too again depends on the initial state as well as on 
the Hilbert space dimension M. For random initial states the power law decay gets faster with 
increasing dimensionality d of the system, and is conjectured to approach a Gaussian decay in 
the thermodynamic limit. 

For mixing dynamics we have a more complicated situation. There are five relevant time 
scales (even six for coherent initial states): The classical mixing time t^ix on which correlation 
functions decay; the quantum decay time of the fidelity Tm; the Heisenberg time tn after which 
the system starts for "feel" finiteness of Hilbert space; the decay time Tp of perturbative Gaussian 
decay present after tn; the time too when the fidelity reaches finite size plateau; for coherent 
initial states we have in addition the Ehrenfest time tE up to which we have quantum-classical 
correspondence. Depending on the interrelation of these time scales, i.e. depending on the 
perturbation strength 6, Planck's constant h and the dimensionality d, we will also have different 
decays of fidelity. All different regimes can be reached by fixing h and increasing 6. Let us follow 
different decay regimes as we increase 6 (shown in Figure 15. 11|) : 

(a) For 6 < 6p we will have tn < t^- This means that at the Heisenberg time, the fidelity 
due to exponential decay (|3.5|1 will still be close to 1, F{tu) ~ 1, and we will see only a 
Gaussian decay due to finite Hilbert space 1)3. 13() . The critical 6p below which we will see 
this regime has already been calculated (|3.14|) and is 
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For 6 < 5p the fidelity will have Gaussian decay with the decay time Tp (|l-i.l,'-{j) 



As we increase 5, we will eventually reach a regime in which we will see initial exponential 
decrease and then at in crossover into a Gaussian decay until the plateau is reached at 

(b) For 5p < 6 < 6s we will have a crossover from the initial exponential decay ()3.5() to the 
asymptotic Gaussian decay (|3.13() at time tu . This regime will take place if Tm < in < ^oo • 
With increasing perturbation, t^o will decrease and the upper border 6s is determined by 
the condition ioo = ^H- Denoting a finite size plateau by F ~ l/M'^, with jj, lying 
between 1 and 2, depending on the initial state (see Section 122^ j we have the condition 
exp (— (tn/'^p)^) = F which gives 

h 



\/o-ciA/ 



x/filnM = 6px/filnAf. (3.52) 



Further increasing the perturbation, we reach perhaps the most interesting regime, in 
which quantum fidelity can decay faster the more chaotic the systems is. In this regime 
the exponential decay persists until the plateau is reached. 

(c) For 6s < 6 < 6mm {6e) we will have an exponential decay until too- The upper border 6raix 
is determined by the condition t^ = imix which is a point where the argument leading to 
the factorisation of n— point correlation function breaks down. For random initial states 
6mix does not depend on h and we get 



^ ' ' ^ (3.53) 



V2(Tcltmix V ^^ 



mix 



Note that the relative size of this window 6mix/6s = ^AA/2/ztmix In AA increases both in 
the semiclassical ^ — > and in the thermodynamic d ^ oo limit. 

For coherent initial states the quantum correlation function relaxes on a slightly longer 
time scale, namely on the Ehrenfest time ^e ~ —lnh/\. Until t-^ quantum packet follows 
the classical trajectory and afterwards interferences start to build leading to the break- 
down of quantum-classical correspondence. Equating t^ = ^e gives the upper border for 
coherent states 

(d) For 6 > 5raix the perturbation is so strong that the quantum fidelity decays before imixj i-e. 
perturbed and unperturbed dynamics are essentially unrelated and fidelity decays almost 
instantly. 

For coherent initial states the upper border of regime (c) is at 6-e which is smaller than the lower 
border of regime (d) Jmix which opens up the possibility of another regime between (c) and (d), 
namely for 6-e < 6 < ^mix the fidelity will decay within the Ehrenfest time. In this regime the 
decay of quantum fidelity is the same as the decav of classica l fidelity and can be explained in 
terms of classical Lyapunov exponents ( Veble Sz Prosenl . |200J) . Note that the relative width of 



this regime (^mix/'^E = V^nTV^OT^^^mix grows only logarithmically in 1/h, i.e. much slower than 
the width of regime (c). 
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In {l/h 



Figure 3.11: Schematic view of different fidehty decay regimes for mixing dynamics. For details 
see text. 



The fidelity decay time scales as ~ 1/5 for regular systems, while it is ~ 1/(5^ for mixing 
dynamics. This opens up an interesting possibility: is it possible that the fidelity would decay 
faster for regular systems than for chaotic ones? The answer is yes. Demanding t^ < t^ we find 
that for sufficiently small 5 one will indeed have faster fidelity decay in regular systems. This 
will happen for 



6< 



6r = nC^/^/2ad0^fi 

6c = n^^' 



' -A-^v', 



iaf, oc ?i^/^ 



4 



random init. state 
coherent init. state 



(3.55) 



We explicitly wrote the result for regular initial states 5^ and coherent initial states 6c as 
the two have different scaling with h. We can see that for random initial states 6^- scales in 
the same way as 5mix and so one has faster decay of fidelity in regular systems just provided 

5 < (5r ~ 5mix- For a coherent initial state this can be satisfied above the perturbative border 

6 > 6p only in more than one dimension d > 1. In one dimensional systems 6p and 6c have 
the same scaling with h and whether we can observe faster decay of fidelity in regular systems 
than in chaotic ones depends on the values of ad and v' . We stress that our result does not 
contradict any of the existing findings on quantum-cla ssical cor respondence. For example, a 
growth of quantum dynamical entropies ( Alicki et al. I . 11996 : Miller Sz Sarkaii Il999il persists 
only up to logarithmically short Ehrenf est time ^f,, which is a l so the upper bound for the so - 
called "Lyapunov" decay of the fidelity (j.Talabert fc Pastawskil . bODll : ICucchietti et al. I . l2nr)2ah 
and within which one would always find F^^^{t) > F'^^^(t) (for coherent states) above the 
perturbative border 6 > 6p, whereas our theory reveals new nontrivial quantum phenomena 
with a semiclassical prediction (but not correspondence!) much beyond that time. If we let 
?i — > first, and then 6^0, i.e. we keep 6 ^ (^r,c(?i), then we recover the result supported 
by classical intuition, namely that the regular (non-ergodic) dynamics is more stable than the 
chaotic (ergodic and mixing) dynamics. On the other hand, if we let (5 — > first, and only 
after that h ^ 0, i.e. satisfying ()3.55|) . we find somewhat counterintuitive result saying that 
chaotic (mixing) dynamics is more stable than the regular one. We can conclude the section 
by saying that we have three non- commuting lim^its, namely the semiclassical limit h ^> 0, the 
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perturbation strength limit 5^0, and the thermodynamic limit d ^ oo, such that no pair of 
these limits commutes. 

Similar regimes as for the fidelity decay were also obtained bv lCoherJ (|l999l . l200Cl ): ICohen &: Kottod 



200(1 ) when studying parametrically driven systems and energy dissipation. 



3.3.1 Illustration with Wigner Functions 

In order to demonstrate faster decay of fidelity for a regular than for a chaotic system we did 
a numerical simulation on the kicked top. For chaotic situation we choose a = 30, while for 
regular one we take a = 0.1 as before. Other parameters are the same for the chaotic and 
the regular case, 7 = 7r/2, S = 100 and a coherent initial state at {'d*,ip*) = 7r(l/\/3, l/\/2)- 
Perturbation of strength 6 = 1.5 • 10~^ is in parameter a. In Figure 13.121 we show the fidelity 
decay for both cases and an illustration of states in terms of Wigner functions. Spin Wigner 
function is a distribution on a sphere with a nice property that if we have two states pi and p2 
and their corresponding Wigner functions Wp-^ and Wp^ the following equality holds, 

tiipiP2) = JWp,Wp,dn. (3.56) 

For an exact definition of a spin Wigner function see Appendix Wigner functions enable 
us to represent fidelity, which is an overlap of the initial state p{0) with the echo state p (t) 
(|2.8|1 . as a phase space overlap integral of two Wigner functions. Also, the classical quantity 
analogous to Wigner function is just the classical density in phase space. But note that Wigner 
function is not necessarily positive whereas the classical density is, so the positivity of a Wigner 



function is a necessary condition (but not sufficient) for the classicality of a state. In Figure I3TT21 
we show for the chaotic (triangles in the fidelity plot and pictures above) and the regular case 
(pluses in the fidelity plot and pictures below) two series of Wigner functions at different times 
(t = 0,60,120): the Wigner function after the unperturbed forward evolution (row labeled 
"forward") and the Wigner function after the echo (row labeled "echo"). We also show the 
structure of the classical phase space. In the inset the same data for the fidelity decay is shown 
on a longer time scale and the vertical line shows the theoretical position of tb- The fidelity is 
the overlap between the echo Wigner function and the initial Wigner function. For our choice of 
coherent initial states, the initial Wigner function is a Gaussian in the semiclassical limit ()2.28|) . 
For chaotic dynamics the forward Wigner function develops negative values around the Ehrenfest 
time after which the quantum-classical correspondence is lost. For the regular dynamics this 

1/9 

correspondence persist much longer, namely until the "integrable" Ehrenfest time ~ h ' after 
which the initial wave packet of size ~ ^ ' spreads over the phase space. For a detailed study of 

J — — 1 r~ ^ 1 1 j \ 

Wigner functions in chaotic systems see (jHorvat fc ProsenL 12003 : Lomba rdi fc Seligmanl . 119931 ) 
and references therein. The echo Wigner function for regular dynamics moves ballistically from 
the initial position, causing the Gaussian decay of fidelity. But note that for regular dynamics 
the echo Wigner function does not have negative values even if they occur in the forward Wigner 
function. In our case the quantum fidelity agrees with the classical one for regular dynamics. 
In a chaotic case on the other hand, the echo image stays at the initial position and diffusively 
decays in amplitude, causing the fidelity to decay slower than in regular case. Classical fidelity 
follows quantum fidelity in the chaotic regime only up to the chaotic Ehrenfest time. 

Finally, we would like to illustrate the dimensional dependence of the fidelity decay. We 
have seen that in d > 1 border b^ is above the perturbative border and so one can get faster 
regular decay in a wider range of b. In addition, strong revivals of fidelity in one dimensional 
regular systems for coherent initial states should be absent for d> \. 
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Figure 3.12: Fidelity decay for chaotic (top curve and pictures) and regular (bottom curve and 
pictures) kicked top. Initial conditions and the perturbation are the same in both cases (see 
text for details). Wigner functions after forward and echo evolution are shown. 



3.3.2 Coupled Kicked Tops 

As already remarked, for a one dimensional system [d = 1), the 'surprising' behaviour of the 
regular decay time being smaller than the mixing one, t^ < t^, is for coherent initial states 
possible only around the border (|3.5U|) 5p (unless a^ is very small) where the exponential decay 
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in the mixing regime goes over to a Gaussian decay due to a finite N . However, for more than one 
degree of freedom, such behaviour is generally possible well above the finite size perturbative 
border (5p. In order to illustrate this phenomenon we will now briefly consider a numerical 



example of a pair of coupled kicked tops where d = 2. 
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Figure 3.13: Fidelity decay for two coupled kicked tops, 5 = S ■ 10~^ and S = 200. The 
upper curves are for e = 20 (mixing regime), the solid curve for a coherent initial state and 
the dashed curve for a random one. The lower dotted curve is for e = 1 (regular regime) 
and a coherent initial state. The exponential and Gaussian chain curves give, respectively, the 
expected theoretical decays (|3.5() and (|3.36|) . with the Gaussian decay time determined by the 
best fit. 

We take a simplified version of coupled kicked tops ()2.25|) with a unitary propagator 

[/(e) = e-if ^iye"'5^2yg-i€5i.52z/5^ (3_57^) 

where S\ and ■S'2 are two independent quantum angular momentum vectors. The perturbed 
propagator is obtained by perturbing the parameter e, so that Us = U{e + 5). The perturbation 
generator is therefore 

V = ^Si,S2., (3.58) 

with h = 1/5. We have used the propagator (|3.57j) over the full (25" + 1)^ dimensional Hilbert 
space, without taking into account the symmetry classes of the double kicked top. 

The classical limit is obtained by 5 ^ 00 and writing the classical angular momentum 
vectors in terms of two unit vectors on the sphere ri^2 = 81^2/ S. In component notation we get 
the classical map 



2^1,2 

2/1,2 

I 

^1,2 



Zl,2 



(3.59) 



yi,2 cos(e2:2,i) + 2;i,2 sin(e2;2,i) 
-xi,2 cos(ex2,i) + yi,2 sin(e2;2,i) 
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We have chosen two regimes, non-ergodic (KAM) regime for e = 1 where a vast majority of 
classical orbits are stable, and mixing regime for e = 20 where no significant traces of stable 
classical orbits were found and where the correlation sum was to a very good accuracy given by 
the first term only 



''»'i''(°> = 2^"'*-*i = liO + l 



(3.60) 



The value of a if the whole correlation sum is calculated is o" = 0.058, only slightly larger than 
1/18 given by C(0). Our motivation here was to compare the regular and mixing fidelity decays 
for the coherent initial state which is here the product of spin coherent states (|2.2f)|) for each 
top 

\'d,^) = \'d2,ip2)®\'&l,Vl). (3.61) 

In Figure 13.131 we show the fidelity decay at 5" = 200 and 5 = 8- 10~^ in regular and mixing 
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Figure 3.14: Numerically calculated decay time of the quantum and classical fidelity for the 
double kicked top. The two full curves show quantum fidelity, the one with a higher slope for the 
chaotic regime and the other for the regular regime. Symbols are decay times of classical fidelity, 
diamonds for the regular and squares for the chaotic regime. The two dotted straight lines are 
theoretical decay time predictions for the regular 1)3. 63(1 and chaotic (|3.62|) regimes. The dotted 
curve (overlapping with squares) is the theoretical classical decay time Tdas = In (0.25/(5) /A. 
Vertical lines show the theoretical position of perturbation borders (|3.64)) . The shading and 
the letters (a),(c) and (d) correspond to the regimes described in Figure 13. Ill The Zeno regime 
corresponds to very short times r < 1 (i.e. strong perturbations 5 > 0.4). 

cases starting from the coherent state (|3.6H) with (-!?!, c/?i) = ('!92,(/:?2) = 7r(l/\/3, l/\/2). We 
find excellent agreement between the theoretical predictions (|3.5() and ()3.36|) and the numerics. 
Note that we are here already in the regime 5 < 5c where the fidelity decay in the mixing 
regime is slower than in regular regime. In mixing regime (e = 20) we show for comparison also 
the fidelity decay for a random initial state for which the decay is (due to ergodicity) almost 
identical to the one for coherent initial state. Overall the fidelity decay here is similar as in a 
one-dimensional case, however, the scaling of various time and perturbation scales on ^ = 1/ S 
is different as discussed in Section [3.31 Observe also that in the regular regime the revivals of 
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m 



fidelity (quantum recurrences at t^) are much less pronounced in d = 2 than in d = 1 (e.g 
Figure 13 .Tf) . 

To furthermore illustrate theoretical regimes of the fidelity decay as explained in the previous 
section (Figure rs.ll|) we numerically calculated the dependence of the time r at which the fidelity 
falls to value 0.37 on the perturbation strength 5 for the double kicked top model (|3.57|) . The 
initial state is the same product coherent state as before. We also computed the decay time for 
the classical fidelity in order to compare the quantum and the classical fidelity. The results for 
S = 100 are on Figure 13.141 For both the classical and the quantum fidelity we calculated two 
sets of data, one for regular regime at e = 1 and one for chaotic regime at e = 5. Numerical data 
is then compared with the theoretical predictions. For the chaotic decay time we use previously 
calculated a = 0.058 ^UTil to get r^ dSISJ 



8.6 



(3.62) 



In the regular case we used fitting of the decay in Figure 13.131 to get the theoretical prediction 
(jHiMl), i.e. to obtain Vv'A-^v', 



4.5 



6VS 



Vv'A~^v' = 0.31. 



(3.63) 



Note that the coefficient Vv'A~^v' has been obtained by numerical fitting only for convenience. 
In principle it could be obtained from classical dynamics, but we would again have to resort to 
numerical calculations as the system at e = 1 is in a mixed KAM-like regime. The values of a 
and Vv'A~^v' can then be used to calculate various perturbation borders as discussed in the 
beginning of Section 13.31 For our choice of S* = 100 we get ()3. 5013. 5513.5413.5^ 



0.0005, dc = 0.0019, 6e = 0.013, (5„ 



0.029. 



(3.64) 
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Figure 3.15: Numerically calculated dependence of quantum (solid lines) and classical (symbols) 
fidelity decay times on the parameter e for the double kicked top. Different curves are for 
different perturbation strengths 6. By increasing e the classical dynamics goes from regular 
to mixing regime, see also right figure showing the dependence on e of numerically calculated 
Lyapunov exponents. By decreasing 6 on the other hand, we go from the regime of quantum- 
classical correspondence for 6 > 6e towards a genuinely quantum regime, where in a chaotic 
regime we can increase the decay time by increasing chaoticity - shaded region for the three 
smallest 5. 
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These theoretical borders are denoted with vertical lines in Figure K-i.141 We can see, that in 
regular regime the quantum and the classical fidelity agree in the whole range of 6 and further- 
more agree with the theoretical t^- In the chaotic regime things are a bit more complicated. By 
decreasing perturbation from 6 = 1 we are at first in regime of very strong perturbation where 
the fidelity decay happens faster than any dynamical scale and it does not depend on whether 
we look at chaotic or regular system or quantum or classical fidelity. This is the so-called Zeno 
regime. For smaller 5 the regular and chaotic regimes start to differ but in the chaotic regime 
we still have quantum-classical correspondence. This correspondence breaks down around (5mix 
where the quantum fidelity starts to follow the theoretical Tm, while the classical fidelity decay is 
Tcias = log (0.25/5) /A, with 0.25 being a fitting parameter (depending on the width of the initial 
packet) and A = 0.89 is the Lyapunov exponent, read from Figure 15. 151 For a brief explanation 
of this classical decay see Section 11.1.21 Incidentally, in our system at e = 5 the classical mixing 
time is very short", tmix ~ Ij and we see that the correspondence breaks already slightly before 
(5e. The quantum fidelity decay time Tm is valid until a perturbative border Sp is reached, when 
a finite Hilbert space dimension effects become important and the decay time becomes equal to 
Tp. Note that for 6 < 5c we indeed have a faster fidelity decay for a chaotic dynamics than for 
a regular one. 

Another interesting aspect of our correlation function formalism is that the decay rate of 
the fidelity in a mixing situation is proportional to the integral of the correlation function a. 
As stronger chaoticity will usually result in a faster decay of C{t) and therefore in smaller a, 
this means that increasing chaoticity (of the classical system) will increase quantum fidelity, 
i.e. stabilise quantum dynamics. Of course, for this to be observable we have to be out of the 
regime of quantum-classical correspondence. All this is illustrated in Figure 15. 15| where we show 
similar decay times as in Figure 13.141 i.e. the same system, initial condition and S = 100, but 
depending on the parameter e for six different perturbation strengths 6. Parameter e controls 
the chaoticity of the classical dynamics. At e = 1 we are in the regular regime and for larger e 
we get into the chaotic regime, also seen from the dependence of the Lyapunov exponent. We 
can see that in the regular regime (e < 2) the classical fidelity agrees with the quantum one 
regardless of 6. In the chaotic regime though, the agreement is present only for the two largest 
5 shown, where we have 6 > 5mix (|3.fc)4|) . For 6 < 6c and for chaotic dynamics (three smallest 6) 
we get into the non-intuitive regime (shaded region in Figure l3.15|) where the quantum fidelity 
will increase if we increase chaoticity. Note that this growth of the decay time stops at around 
e ~ 4 because the classical mixing time tmix gets so small that the transport coefficient is given 
by its time independent term a = C(0)/2 alone. 



"This is the reason why a 13.601 1 is almost entirely given by C(0). 



Chapter 4 

Special Case: Zero Time Averaged 
Perturbation 



One should always keep an open mind, but 
not so open that one's brains fall out. 

— Bertrand Russell 



In the previous chapter we considered general perturbations, for which the double correlation 
sum was growing with time either linearly for mixing dynamics or quadratically for regular 
systems. What about the third possibility, namely if it does not grow with time. This situation 
will be the subject of the present chapter. We will demand that ($]^(i)) ~ t^ for any initial 
state, i.e. all matrix elements of S^ must be 0{t^), 

E2;,(i)oct°, j,k = l,...,J\f. (4.1) 

We will write the perturbation V as the sum of its time average V and the rest, called the 
residual part Vres 

V=:V + V,es- (4.2) 

For a nondegenerate* spectrum of the unperturbed propagator Uq we have seen that the time 
averaged perturbation equals to the diagonal part of V ()3.18|) therefore, the residual V^es is 
just the off-diagonal part of V and has zeros on the diagonal, (Vres)fcfc = 0. The operator S(t) 
can then be written as T,{V,t) = Vt + 0{t^), where the second part depends just on T4es and 
does not grow with time. One can conclude that to satisfy condition (|4.H) the time averaged 
perturbation must be zero, V" = 0. The subject of the present chapter are perturbations with 
V = 0, also called residual perturbations because V = T^es- For non-residual perturbations all 
essential physics was contained in the operator S(i). For residual perturbations on the other 
hand, the second term involving T(t) in the BCH form 1)2.14(1 of the echo operator will also be 
important. Matrix elements of r(i) in the eigenbasis of Uq are 

^^^^1^ = lT.\y^k\'cot[l{cPk-4>,)]+0{t-'), (4.3) 

t n 2sm[^[(j)k- (pj)\ 



*For degenerate spectra we have V = J2k.i ^^'t''' ~ 4>i)Vki\'l)k){4>i\- 
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Note that the matrix elements of T{t) can not grow faster than linearly with t, despite the 
double sum over time in the definition of T. We see that, provided the perturbation is residual, 
the limit of time-averaged T(t) defined as 



t-i t-i 



hm I^ = A lim \Y.Y. [^(*'), V{i!')] (4.4) 



i'=Ot"=t' 

exists and is diagonal in the eigenbasis of Uq: 

r = Y.^n\cl^,)m, fn = \Y. \V,k? coi[\{cku - </>,)]• (4-5) 

The operator T is again a constant of motion, [?7o,r] = 0. 

Any residual perturbation can be defined in terms of another operator W by the following 
prescription 

V = W(l) - W{Q), W{t) := U^^WUl (4.6) 

In the continuous time case we would have definition V =: {d/dt)W = ^[Hq, W]. Indeed, given 
a residual perturbation one easily determines the matrix elements of W as 

^^^ '■= /v^^'\^^ 1 - (^•'^) 

exp {i{(j)j - (pk)) - 1 

Note that instead of the unperturbed evolution operator Uq one could use any other unitary 
operator that has the same degeneracy and eigenvectors as Uq in the definition of V in terms 
of W (|4.6|) . With the newly defined operator W, the expression for S(t) is extremely simple, 

S(t) = W{t) - W{0). (4.8) 

Similarly, the expression for T(t) is also considerably simplified, 

r(t) = SR(t) - ^[H^(0), W{t)], R := ^[W{0),W{1)], (4.9) 

and 

t-i 

SR(t) := J2 «(*'), Rit) ■■= Uo'RK (4.10) 

t'=0 

The operator T(t) is, apart from the term ^[W, W{1)], similarly as S(t) the sum of R{t). In the 
continuous time case we have R := ^[W, {d/dt)W] = n~^[W, [W,Ho]] and r(t) = /g dt'ii(t') - 
■^[W, W{t)]. The fidelity decay will be given by the echo operator 



Ms{t) = exp |-^ (^S(t)<5 + ^Tit)5^ + •••)}' (^"^^^ 



with S(t) (|4.8|) and T{t) (|4.1()j) expressed in terms of the operator W. We see that for small 
times t < t2, with ^2 ~ 1/^, the second term involving T{t) can be neglected. Therefore, for 
t < t2 the fidelity amplitude is simply 

fit) = (exp i-i6iWit) - Wm/h)) . (4.12) 

Expanding f{t) to the second order in 6, we find F{t) = I — t7{kq + k^ — C{t) — C{t)*) where 
nl ■= {W'^{k)) - {W{k)f, C{t) := {W{t)W{<d)) - {W{t)) (VF(0)). Using the Cauchy-Schwartz 
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inequality |C(t)| < K-oKt and the fact that for a bounded operator W the sequence Kt is bounded, 
say by r, we find a freeze of fidelity 

l-F(t) <4-2r2, t<t2 (4.13) 

n 

for arbitrary quantum dynamics, irrespective of the existence and the nature of the classical 
limit. For residual perturbation the fidelity therefore stays high up to a classically long time 
t2 and only then starts to decay again. After t2 the second term T{t) will become important 
and this will cause the fidelity to decay. We have to stress that the freeze of fidelity is of purely 
quantum origin. The classical fidelity does not exhibit such behaviour. 

Although residual perturbations might seem artificial at first sight, there are two possibilities 
how they can arise. First, the average perturbation V commutes with the Hamiltonian Hq 
generating Uq and can sometimes be put together with the original Hq into the unperturbed 
Hamiltonian, thereby resulting in a residual perturbation. This moving of V into Hq just changes 
the eigenenergies and is often done in various mean field approaches. Another possibility to 
have y = is due to symmetries. For instance, having a unitary symmetry P commuting with 
the unperturbed evolution, [Uq, P] = 0, and perturbation V anticommuting with the symmetry, 
VP = —P V, will result in a residual perturbation. But note that in order to have a well defined 
operator W (|4.7j) with no singularities near the diagonal, so that our theory can be applied, the 
matrix elements of V must increase smoothly away from the diagonal. 

In the next two sections we will discuss two examples for which one can use the echo 
operator (|4.11|) to calculate fidelity decay to all orders in 5. This can be done for completely 
mixing dynamics and for the opposite extreme of regular dynamics. In both cases we will also 
assume the operators V, W, and therefore also R, to have well defined classical limits, so we 
will be able to use semiclassical arguments. For small times t < t2 the fidelity will freeze to a 
constant value - the plateau. The value of the plateau will be determined by the non-increasing 
S(t), 

^exp(-i-{Wit)-W{0)i 



i"plat(i) 



h 



(4.14) 



For large times t > t2 instead, the fidelity decay will be determined by the operator R in r(t). 

4.1 Mixing Dynamics 

4.1.1 The Plateau 

We begin with the linear response evaluation of the fidelity plateau (|4.14() . For classically 
mixing dynamics we can assume that for times larger than some mixing time ti the time 
correlations vanish in leading semiclassical order, C{t) -^ 0{fi). Also, quantum expectation 
values become time independent and equal to their classical values in leading order. We will 
denote by {A)^^ •= /d/i^d the classical average value of observable A. Therefore, for times 
ti < t < t2 the linear response expression for fidelity plateau, i.e. its time-independent value is 

i^plat = 1 - ^(k2 + 4), ^2 .^ ^^2(^)^ _ ^^(^)^2 ^ (^^^^^ 

where k^j := {w^)^i — (^)ci i^ time-independent classical limit of k^ with the classical limit 
w{q,p) of operator W. We will consider two kinds of initial states, namely coherent states and 
random states. For coherent initial states (CIS) the mixing time ti is equal to the Ehrenfest 
time, ii ~ — In fi/X. As the initial state variance of W is proportional to the spread of the 
CIS, Kq cx h, it can be neglected with respect to k^j. The fidelity plateau for CIS is therefore 
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FS^ = 1 — (5/?i)^K^[. For random initial states (RIS), the mixing time is h independent, ti (x h . 
Also initial state average for RIS is equal to ergodic average and the fidelity plateau is in this 
case -Fi}at — 1 — 2(5/7i)^K^j. Summarising, the linear response value of the fidelity plateau is for 
CIS and RIS 

1-^pM-^'^ci, l-Fp^i?-2^4, (4.16) 

i.e. it is twice as large for RIS than for CIS. 

One can go beyond linear response in approximating (|4.14j) using a simple fact that in the 
leading order in h quantum observables commute, and as before, that for t > ti the time 
correlations vanish, namely (exp(-if (W(t) - W))) = (exp(-if W"(t)))(exp(if VF)). For t > h 
expectation values become time independent and so 

i^piat = Kexp {-iw6/h))^, (exp {iW6/h))f . (4.17) 

Note that the the right average is an average over initial state whereas the left average is the 
classical average over the invariant measure. The fidelity plateau can be compactly expressed 
in terms of a generating function G{5/Ti) 

G{z) := (exp(-izu;))^i = ^ j d'^qd'^p ei,v{-izw{q,p)}. (4.18) 



For CIS one can neglect the initial state average over a localized packet, i.e. the second term, 

and gets 

i-, 

plat 

For RIS on the other hand, initial state average is equal to the classical average in the leading 
order and we have 

F^i^\G{6/h)\\ (4.20) 

We have a universal relation between the plateau for CIS and RIS, namely F^^ = {F^^^)'^. If 
the argument z = 5/Ti of the generating function is large, the analytic function G{z) can be 
calculated generally by the method of stationary phase. In the simplest case of a single isolated 
stationary point x* in N dimensions we obtain 



F^^ = \G{6/n)\\ (4.19) 



\G{z) 



n 

Yz 



^/2 -1/2 

deta,,a,„t/;(a;*) . (4.21) 



This expression gives an asymptotic power law decay of the plateau height independent of the 
perturbation details. Note that for a finite phase space we will have oscillatory diffraction 
corrections to Eq. (|4.21|) due to a finite range of integration /d/u which in turn causes an 
interesting situation for specific values of z, namely that by increasing the perturbation strength 
5 we can actually increase the value of the plateau. 

We tested the above theory by numerical experiments. For the system we choose the kicked 



top, with slightly different order of factors as before (|2.20|) . One step propagator is 

Uo := exp f-i«^ j exp (-if^y) . (4.22) 

We take a = 30 ensuring fully chaotic classical dynamics whereas spin size is taken S = 1000 
giving effective Planck's constant fi = 1/ S = 1 • 10^'^. Again we choose two initial states, a 
random one and a coherent one centred at (-(?*,(/?*) = (1,1) and both are projected onto OF 
subspace (jTM . To get a residual perturbation we take W := 5^/25^ so that 1^(1) = S'^/2S^ 

giving the perturbation generator 

q2 _ ^2 
^:=^^^- (4.23) 
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Figure 4.1: F{t) for the kicked top with perturbations 5 = 10~^ (a), and 6 = 10~^ (b). In each 
plot the upper curve is for CIS and the lower for RIS. Horizontal chain lines are theoretical 
plateau values, linear response H4.16() in (a) and full H4.2UI4.19|1 in (b). Vertical chain lines 
are theoretical values of ^2 ()4.29|) . The full circles represent calculation of the corresponding 
classical fidelity for the CIS which follows quantum fidelity up to the Ehrenfest (log h) barrier 
and exhibits no freezing. 



The perturbed propagator is as always Us = Uoexp{—i5V/Tl). We choose two perturbation 
strengths. A weak 6 = 10"'^ to check the linear response expressions (|4.16j) . The results are 
seen in Figure l^lTk . The classical value of k^j used in theoretical formulas is easily calculated for 
w = z^/2 and one gets k^j = 1/45. In Figure l^?Tb we show numerical results for larger 6 = 10~^ 
where the plateau is very low and full formulas must be used. Generating function ()4.18|) can 
be calculated exactly without stationary phase approximation, resulting in 



G{5S) 



^ erfi f e-/^ 



26S 



6S/2 



(4.24) 



We have an asymptotic power law decay as predicted by the general stationary phase formula^^ 

(ion 



In addition, due to a finite phase space we have a diffractive oscillatory erfi correction. 
Theoretical prediction for the plateau using G{5S) (|4.24|) agrees well with the numerical results 
shown in the figure. Small quantum fluctuations around the theoretical plateau values lie beyond 



^Note that w = z^/2 has one stationary point in A'' = 1 dimension, despite the phase space being two- 



dimensional. 
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the leading semiclassical approximation used in our theoretical derivations. Note also that the 
quantum fidelity and its plateau values have been expressed (in the leading order in h) in terms 
of classical quantities only. Yet, the freezing of fidelity is a purely quantum phenomenon as one 
can also see in Figure WTH where the classical fidelity does not exhibit freezing. 

4.1.2 Beyond the Plateau 

Next we shall consider the regime of long times t > ^2 • Then the second term in the exponential 
of the echo operator (|4.11j) dominates the first one, however even the first term may not be 
negligible for large 6. Up to terms of order 0{t6^) we can factorize Eq. (|4.11jl as Ms{t) ~ 
ex.p{—i-^{W (t) — W)) ex.p{—i^r{t)) . When computing the expectation value we again use the 
fact that in leading semiclassical order operator ordering is irrelevant and that, since t ^ ti, 
any time-correlation can be factorized thus the second term ^[W^, VF(t)] of T{t) ()4.9() vanishes 
and we have 

exp -i;^SR(t) 



F{t) ^ Fpiat 



2h "' 7 



t > t2. (4.25) 



This result is quite intriguing. It tells us that apart from a pre-factor Fpiat, the decay of 
the fidelity due to a residual perturbation is formally the same (in the leading semiclassical 
order when time ordering is not important) as the fidelity decay with a generic non-residual 
perturbation, eq. (|2.11|) . when one substitutes the operator V with R and the perturbation 
strength 6 with 6ji = (5^/2. Thus we can directly apply the semiclassical theory of fidelity 
decay for general perturbations explained in Section 13.11 using a renormalised perturbation R 
of renormalised strength 5r. Here we simply rewrite the key results in the 'non-Lyapunov' 
perturbation-dependent regime, 5fi < fi. Using the classical transport rate ur, 

an := lim ^^r(^))c1 " (^i^)ci ^ (4.26) 

t^oo 2t 

we have either an exponential decay 



F{t) ^ Fpiat exp I --^cTRt 1 , t < in, (4.27) 

or a (perturbative) Gaussian decay 

F{t) ^ Fpiat exp f -— ^aR— j , t > tn- (4.28) 

The crossover between the two decay regimes happens at the Heisenberg time tn = M jl (|3.1fl|) . 
Of course, the same consideration regarding the asymptotic saturation value of the fidelity F 
due to a finite H applies here as well (see discussion of time scales in Section IX^ . The prefactor 
-^piat can be calculated as described in the previous section, Eq. (|4.17j) . and depends on the initial 
state. The exponential terms of 1)4. 2714. 28|) on the other hand do not depend on the initial state. 
Prom the two possible regimes of long-time fidelity decay we can also more precisely specify 
time ^2 when the plateau ends. Comparing the two factors in 1)4. 2714. 28|) with the plateau value 
(Eq. 14.161 for 5 < /i), we obtain a semiclassical estimate of f2 



,.^.„^,,^|,_^J. ,4,29) 

In Figure 14.11 we can see nice agreement of theoretical t2 with the results of the simulation 
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Figure 4.2: Long time Gaussian decay ()4.28|) for the same CIS and parameters as in Figure l^lTk . 
The fuU curve is a direct numerical simulation, empty circles are numerical calculation using a 
renormalised perturbation R of strength 5r (|4.3U|1 . while the chain curve gives the theoretical 
decay. 



Interestingly though, as we have another time scale 12-, the duration of the plateau, not 
present for a general case of non-residual perturbation, the exponential regime ()4.27|) can only 
take place if ^2 < ^H; otherwise we immediatelly get a Gaussian decay after the plateau. If 
one wants to keep -Fpiat ~ 1 (i-e- high) and have an exponential decay in the full range until 
the asymptotic F ~ 1/AA, the condition on dimensionality is imposed. Namely, demanding 
high plateau S/fi < 1 for t < t2 and low fidelity in the limit ?i — > at the Heisenberg time, 
F{tu) = ex.p{—{6/h)^aji/4:h ~ ), gives condition on dimensionality d > 2. In one dimensional 
systems and in the semiclassical limit the exponential decay (|4.27)) therefore can not be seen. 

We again compared theory with numerics. The system is the one dimensional kicked top 
already used in the previous section describing the plateau. The perturbation generator R for 
our choice of W is 

R = -^{S.SyS, + S^SyS^), (4.30) 

having a classical limit R^i = r = —xyz. We numerically calculated the integral of the classical 
correlation function (|4.26|) giving the transport coefficient ur = 5.1-10"^. In Figure lX2l we show 
three different sets of data. Direct numerical calculation, "a renormalised calculation" obtained 
by taking perturbed dynamics as Us = Uoex.p{—i5^R/fl), i.e. perturbation generator R with 
the strength (Jr = 5'^ /2 and theoretical Gaussian decay 1)4. 28(1 where all parameters have been 
calculated classically. Apart from the plateau prefactor Fpiat (which is close to 1) the decay 
does not depend on the initial state. 



Double Kicked Top 

To demonstrate that for d > 1 we can also get an exponential decay of the fidelity, we also 
consider a two dimensional system. We look at two {d = 2) coupled tops Si and 5^2 described 
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Figure 4.3: Long time fidelity decay for two coupled kicked tops. For strong perturbation 
6 = 7.5 • 10~^ in (a) we obtain an exponential decay 1)4. 27() . and for smaller 6 = 2- 10~^ in (b) 
we get a Gaussian decay (|4.28|) . The three curves have the same meaning as in Figure 14.21 



by a propagator 



Uq = exp (-ie5'ziS'z2) exp (-i7r5yi/2) exp (-i7rS'y2/2), 



(4.31) 



with a perturbation generated by 



2 /r,c2 



W = Ai^l + l^A2, A = St/2S 



(4.32) 



We set S = 1/Tl = 100, and e = 20 in order to be in a fully chaotic regime. The initial 
state is always a direct product of spin coherent states centred at ((/?i,2,^i,2) = (1)1) which 
is subsequently projected on an invariant subspace of dimension M = S{S + 1) spanned by 
{Woe ® "Wrlsym, where TC^ = Ti. \ Ti-oE and {-jsym is a subspace symmetric with respect to 
the exchange of the two tops. The results of numerical simulation are shown in Figure 14.31 
We show only the long-time decay, as the situation in the plateau is qualitatively the same 
as for d = 1. For sufficiently large perturbation one obtains an exponential decay shown in 
Figure 14.3b . while for smaller perturbation we have a Gaussian decay shown in Figure 14.3b . 
Numerical data have been successfully compared with the theory 1)4. 2714. 28|) using a classically 
calculated ur = 9.2- 10~^ together with theoretical -Fpiatj and with the "renormalised" numerics 
using the operator R, similarly as in Figure W^ for a one dimensional system. 
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4.2 Regular Dynamics 

The procedure of calculating the fidelity decay for the case of regular dynamics will be similar 
as for general perturbations in Section f3. 21 We will use the classical action-angle variables (for 
definition see Section [3. 2. 1() and semiclassical methods to calculate the fidelity plateau as well 
as its long time decay. The only difference will be that in contrast to the case of a general 
perturbation, where only $](t) was important in the BCH form of the echo operator, here the 
second term involving T{t) will also be relevant for long times. The same form of the echo 
operator (|4.11|) has already been used for the case of a residual perturbation in mixing systems 
as described in the previous section. The existence of action-angle variables enables us to expand 
the classical limit v of the quantum perturbation generator V into Fourier series 

V{j,e):= Y.^rrrUy'^-'. (4.33) 

The fact that the perturbation is residual is reflected in the zeroth Fourier coefficient which is 
zero, vq = 0, so we explicitly excluded this term from the summation. Classically this means 
that the perturbation only changes the shapes of tori as the average change of the frequency 
along the unperturbed tori is z ero. For the ex planation of the decay of classical fidelity in case 



of a residual perturbation see (|Benenti et al. . . ■2003a '). Similarly, we can expand the classical 



limit w{j^ 0) of a quantum observable W , 

w{j,e):= y^Wrr^Uy^-r (4.34) 






The zeroth Fourier mode wq can be set to zero as it cancels in the definition of v in terms of 
w. Using these expansions we can easily calculate the leading semiclassical forms of S(t) and 
r(t). For T.{t) (glHI) we have 

S(i) ^ w{J, + u;{J)t) - w{J, &)= Yl Wm.iJW'^"^^''^' - 1)6*"^-®. (4.35) 

Note that this is still an operator (capital J and 0) and is correct in the leading semiclassical 
order (the sign =). Coefficients Wm can be also expressed in terms of Vm as 

wmij) = -ie-i— 0-)^/^ . ;-^-^") (4.36) 

2sm(im • u){j)) 

In the continuous time case we have Wm = —Wm/i^ ' ^• 

The operator S(t) will determine the plateau. Long time decay on the other hand will be 
dictated by the operator T{t). For long times when this decay will take place we can define 
an average F ()4.4|) and approximate r(t) ~ Ft. This approximation is justifiable as we will 
see, because the fidelity decay will happen on a long time scale ~ l/(^^, whereas the average 
of F(i) converges in a much shorter classical averaging time tave- As one can see from the 
definition of F(i) in Eq. (|4.9jl . the average F is nothing but the time averaged operator R. For 
R the semiclassical limit r can be calculated using the Poisson brackets instead of commutators, 
r = — {w, w{l)}. When we average r over time, only the zeroth Fourier mode survives resulting 
in 

f (J) ^ f{J), r{j) = -J2m-dj{ \wm{j)\^ sin [m • lj{j)]} . (4.37) 

In the continuous time case we would have f = — J2m^o "^ " ^j {l^mP m ■ u}. 

From the equation for Wm (|4.36j) in terms of Vm we see that the denominator sin (m • uj/2) 
could cause problems at points of vanishing m ■ u:{j) = 0. This divergence carries over to all 
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operators like f or S(i). In classical perturbation theory this so-called "small denominator" 
problem is well known and cannot be avoided. For our quantum mechanical calculation though, 
there is an easy remedy. It is sufficient to remember that we are dealing with a residual 
perturbation, i.e. the one for which all matrix elements Vjk are zero if the corresponding 
eigenphases are equal, (pj = (pk- But the term m • a; is nothing else than the semiclassical 
approximation for (pj — cpk. Therefore, one can see that actually all the diverging terms have to 
be excluded when evaluating expectation values. This is furthermore confirmed if we compare 
the semiclassical expression for f (|4.37|) with the matrix elements of operator T{t) ()4.3() . An 
example of such a calculation will be presented when describing the decay for random initial 
states. 

Using semiclassical expressions for S(t) and T we can write the echo operator as 

M5(i) = exp|-^ (^^{t)6 + ht6^^y (4.38) 

with r(t) and T given in Eqs. (|4.35|4.37j) . The third order term in the BCH form of the echo 
operator ()4.38|) can be shown to grow no faster than 6'^t and can therefore be neglected. For 
times smaller than t2 (specified later) the term involving T can be neglected and the fidelity 
will exhibit freeze. The fidelity amplitude of the plateau is /(t)piat — (exp (—i6T,{t)/h)) with 
semiclassical S(t). Furthermore, for t > ti the fidelity plateau is time independent and can be 
calculated by averaging the fidelity amplitude over time 

1 * 

/piat = /hn-^/(0. (4.39) 

"^°° t'=0 

Time averaging the classical version of S(t) = Y,{j,6{t)) (|4.35|) is equivalent to averaging over 
the angle 0, resulting in the fidelity plateau 

/plat = /exp (i^wiJ, 0)) J ^ exp (-i|^( J, ^)) ) , (4.40) 

where we assumed ergodicity oi m ■ u:t so we could replace a sum over t with an integral over 
angles. 

For long times t > ^2 the second term in the echo operator becomes dominant. If the plateau 
is small {6 < h), the first term with S(t) can be neglected and the fidelity is 



F{t) 



'"'A-'U^^-^'^ 



(4.41) 



For strong perturbations, when the plateau is not negligible, both terms should be retained. But 
as opposed to the mixing situation where we could factorize the two contributions, here one can 
not make any general statements. In the following theoretical calculations we will focus on the 
case when the plateau can be neglected in comparison with the long time decay. The formula 
1)4.41(1 is completely analogous to the long time fidelity decay for a non-residual perturbation 
and we can use the same ASI representation in terms of an integral over classical actions. The 
role of V is played by F and the perturbation is here 5^/2. The fidelity amplitude in the ASI 
formulation is therefore 

fit) - n-^d^jexp (-i^^j)) dpij), (4.42) 
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with dp{j) being the classical limit of Dp{fin) = {n\p\n) and recall that r is the classical limit 
of r H4.37() . The AST representation is valid in a time range ti < t < tg,, where the upper limit 
ta is determined by the variation of the phase over one Planck's cell, 

Before going on with the evaluation of the plateau (|4.4()|) or of the long time decay (|4.42|) 
for random and coherent initial states, let us justify why we were allowed to time average the 
plateau. 

Justification of Time Average Plateau 

By expanding the echo operator into power series over 5, the calculation of the plateau is turned 
into calculation of expectation values of moments S'^(t), 

k 

S'=(t)^ Y, l[wrr^,{J){e''^'^'^^■^^' -l)e''^'-®. (4.44) 

TTll,...,Tn.fe^O 1 = 1 

We can average over the initial density matrix p with matrix elements Pn,n' in the action 
eigenbasis. Taking into account that the exponential of operator acts like a shift operator 
p. 26(1 and that eigenvalues of J are h{n + a), we see that the expectation value will be a sum 
of terms of the form 

^ Y.g{hn)e'^-'^^^^^'pr,,r^+m', (4.45) 

T7ai,...,m.fe7^0 'IT- 

where the function g depends on indices mi, index m' := J2i "i^l ^^'^ ^^ •= J2i simi with s; = 
or 1 depending on which terms we take from the product in Eq. (|4.44() . We also neglected Maslov 
indices as they are negligible in the leading semiclassical order. For our derivation only the inner 
sum over n is important. Let us consider two cases of initial states p, random initial state and 
coherent initial state. 

For a random initial state we average over a random ensemble, resulting in Pn^m -^ ^n,m/-^- 
From this we immediately conclude that index m' in Eq. ()4.45|) must be zero, and as a con- 
sequence also m, = is zero. Therefore, for random initial states only terms with all mi = 
survive the averaging over random ensemble (which is the same as if we would average over time 
instead of over random ensemble). The time scale ti after which approximation with a time 
average is permissible is determined just by frequencies, ti ~ l/\m ■ u:\ and does not depend on 
h or perturbation strength 5, i.e. for random initial states we have 

h ~ 5°n°. (4.46) 

For localized initial states a little more work is needed to show that the expectation value 
equals the time average. For definiteness we will consider coherent initial states, for which the 
matrix elements have a Gaussian distribution (JSI32I)- Furthermore, we will assume that the 
number of relevant Fourier modes is smaller than the width (in n) of the initial packet. As 

1/9 

the width of the packet grows as h ' this is justifiable in the semiclassical limit provided 
w{j,6) is sufficiently smooth, so that its Fourier coefficients decrease sufficiently fast. In this 
approximation we have 

Pr^,r^+rr^' = t^p(?^n)e'^'■^^ (4.47) 

with dp a classical limit of the quantum structure function (|,3.33|) . The inner summation over 
quantum actions hn in Eq. ()4.45|) can now be replaced with an integral over the classical action 



= 9[ 
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j and the method of stationary phase can be used to calculate the resulting integral. Expanding 
frequencies around the centre of the packet j* , ijj{j* + a;) = i^{j*) + i^x + • • •, where 

is a matrix of frequency derivatives, we can calculate the sum 

^^■*ymMr)t (!!l\ |detA|^/2 f d'^xexp (-^x-Ax + itm-nx 

= 9(i*)e*"^'^(^'*)*exp \-^rn ■ nA'^n^rnj . (4.49) 

We see that all terms with m ^ decay to zero. The longest decay time scale of Gaussian 
envelopes is estimated as 

1/9 

- min (m ■ nA-^n^m)] oc h~^/'^. (4.50) 

4 m^o V J J "■ ' 

This means that for times longer than ti only term with zero m will survive which is equivalent 
to time averaging the classical expression. Note that the Gaussian decay (|4.49|1 is absent if f2 = 0, 
e.g. in the case of a d-dimensional harmonic oscillator. Fidelity decay for a residual perturbation 
in a harmonic oscillator will be discussed for a Jaynes-Cummings model in Section [5.51 There 
may also be a general problem with the formal existence of the scale ti 1)4. 50() if the derivative 
matrix fi is singular, but this may not actually affect the fidelity for sufficiently fast converging 
or finite Fourier series ((2221) • Time ti can be interpreted as the integrable Ehrenfest time, up 
to which quantum-classical correspondence will hold. After ti a quantum wave packet of size 
~ \/n will spread over a classically large region of size ~ h and interference phenomena will 
become important. This will also be reflected in the fidelity. As the perturbation Wm will couple 
different tori, the echo packet will also start to exhibit interferences after ti and therefore we 
can expect agreement between quantum and classical fidelity only up to time ti. This must 
be contrasted with the case of a general non-residual perturbation, where there were no upper 
limits on the correspondence between quantum and classical fidelity for localized initial packets 
under certain conditions. 

4.2.1 The Numerical Model 

For the sake of numerical demonstration in the present section of residual perturbations in a 
regular dynamics we take an integrable kicked top model with a slightly different unperturbed 
one step propagator 

C/o = expf-i|5||-/3}'V (4.51) 

where we introduced a second parameter /?. The classical inap corresponding to Uq is simply a 
twist around z-axis 

xt+i = xt cos {a{zt - l3)) -yt sin {a{zt- P)) 

yt+i = yt cos {a{zt - P))+xt sin {a{zt- (3)) (4.52) 

Zt+l = Zt. 
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Unperturbed evolution represents a continuous time system with the classical Hamiltonian /iq 
generating a frequency field uj{j) 

hoij) = l^ij - P?, u;{j) = a{j - 13). (4.53) 

Here we used a canonical transformation from a unit-sphere to an action-angle pair (j, 9) G 

[—1,1] X [0,27r), namely 



x = Jl—pcos6, y = Jl — p sinO, z = j. (4.54) 

We perturb the Hamiltonian by periodic kicking with a transverse pulsed field in x direction, 

1 / °° 

hsij,e,r) = -a{j-l3f + 6^1- P cose ^ 5(r - fc). (4.55) 

fc=— oo 

Perturbed quantum evolution is given by a product of two unitary propagators 

Us = Uoexp{-[6S,), (4.56) 

so the perturbation generator is 

V = SJS. (4.57) 

The classical perturbation has only one Fourier component, namely 



viJ,e) = ^l-pcose, v±iij) = -^^l-f, (4.58) 

whereas zeroth Fourier component is zero, vq = 0, indicating that we have a residual perturba- 
tion, V = and V = 0. The classical limit w of W is also readily calculated, giving 

/ ■ /.^ 1 r, :; sin (6* - w/2) ,^ ^^, 

2 ^ sm (uj/2) 

with u> = a{j — 0). The reason we introduced parameter (3 is to be able to control a possible 
singularity in w{j,9) at points where uj = 0. 

4.2.2 Coherent Initial State 

The Plateau 

For times t < ti quantum fidelity will follow the classical fidelity and will be system specific. 
For times ti < t < t2 quantum fidelity will exhibit the plateau while the classical fidelity will 
continue to decay. To calculate the quantum plateau we have to evaluation the time-average 
formula for the plateau (14.40(1 . We shall make use of the fact that for coherent states we have 
the expectation value 

{expi-{i6/h)g{J,&))) ^ exp{-{i6/n)g{r,e*)). (4.60) 

for some smooth function g, provided that the size of the wave-packet ~ vh is smaller than the 
oscillation scale of the exponential ~ h/6, i.e. provided 5 <C h}''^. Then the squared modulus 
of /plat 11301) reads 



pCIS 
r. 



plat /'9„\2d 



(27r)2 



Jd''eexp(^-'-^w{r,e) 



(4.61) 
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It is interesting to note that the angle 6* does not affect the probability -Fpiat as it only rotates 
the phase of the amplitude /piat • Note the similarity of this result with the plateau for coherent 
initial states and mixing dynamics, Eq. (|4.19|) . In the mixing case we had an equilibrium average 
of w in 2d dimensional phase space, whereas for regular dynamics we have an average over only 
d dimensional angle-space of 0. For weak perturbation 6 < h only the linear response expression 
for the plateau is needed. Expanding general -Fpiat we obtain 

1 - ^PM = ^^cis, i^cis = E \^mif)\'. (4.62) 

Results of the numerical simulation for (3 = 0, a = 1.1 are shown in Figure ^31 For our 
choice of the perturbation 1)4. 59(1 the integral in the semiclassical expression of the plateau 1)4. 61(1 
is readily calculated. Actually, the calculation of -Fpiat can be analytically carried out for any 
perturbation with a single nonzero Fourier mode ±mo, with the result 



^pU? = ^o(2-|«;^„(r)|j, (4.63) 

were Jq is the zero order Bessel function. For a more general multi-mode perturbations we have 
to evaluate the integral (|4.61l) numerically. For our single mode perturbation w 1)4. 59() we have 



plat -^o y"^2sm{aj*/2) 
The linear response expansion of this general result reads 



^P^i? = 4 ^S ^frT^J ■ (4-64) 



1 - r' 



F-^ = i5Sf ^c.s, ^c.s = ^-^^y (4.65) 



For numerical illustration in Figure 14.41 the initial coherent packet has been placed at 
{'&*,ip*) = (1,1). We can see that until time ti (|4.5fl|) quantum fidelity follows the classical 
one (circles in Figure ^71]) . After ti quantum fidelity freezes at the plateau, whose value nicely 
agrees with the linear response formula (|4.65|) or with the full expression (|4.64j) for strong per- 
turbation. Vertical chain lines show theoretical values of t2, which is the time when the plateau 
ends. Also, at certain times the quantum fidelity exhibits resonances, i.e. strong revivals of 
fidelity. These "spikes" occurring at regular intervals are prominent also in a long time de- 
cay of fidelity in Figure 14.51 These will be called the echo resonances and are particular to 
one-dimensional systems. 

Long Time Decay 

To obtain long time decay of quantum fidelity one has to evaluate the ASI in Eq. (|4.42|) . Ev- 
erything is analogous to the case of general non-residual perturbation described in Section [3.21 
In the formula for F{t) = exp (— (i/Tr)^) 1)3. 36() we only have to replace the perturbation 5 with 
(5^/2 and v with f, so that the fidelity decay is 



Fit) ^ exp {-(t/r„)^}, r„ = ^^j^r^^ (4-66) 

with the semiclassical f given in Eq. ()4.37|) and its derivative f' := djf. The decay time 
scales as r^r ~ ^ ' S~'^ and is thus smaller than the upper limit t^^ r^ fi (5~^ of the validity of 
the stationary phase approximation used in deriving the Gaussian decay. Remember that the 
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Figure 4.4: Short time decay of the fidehty for a quantized top is shown for the coherent initial 
state, for S = 200 (a), and S = 1600 (b), with a fixed product 5S = 0.32, where the plateau is 
well described by linear response. In (c) we show S = 1600 and a stronger perturbation with 
6S = 3.2. Note that the time axis is rescaled as t/ti. Symbols connected with dashed lines 
denote the corresponding classical fidelity. The horizontal chain line denotes the theoretical 
value of the plateau (|4.64j) , while the vertical chain line denotes the estimated theoretical value 
for ^2 ()4.67|) . In (b,c) we also indicate fractional 2iTk/p resonances with k/p marked on the 
figure (see text for details). 



above formula is valid only if the plateau is small and the $](t) term can be neglected, i.e. for 

— 1/2 

6 < i^Qjg h. For such small 5 the crossover time t2 from the plateau to the long time decay can 
be estimated by comparing the linear response plateau (|4.65j) with the long time decay (|4.66|1 . 
{S/hyucis ~ (*2/Trr)^, resulting in ^2 ~ Tj-j-d^/ucis/h ~ h~^''^6~^. For stronger perturbations, 
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Figure 4.5: Long time ballistic decay of the fidelity for the kicked top with coherent initial state 
is shown for cases S = 200 (a), and S = 1600 (b), of weak perturbation 6S = 0.064, and for 
strong perturbation S = 1600, 6S = 3.2 (c). Chain curves indicate a theoretical Gaussian 1)4. 66(1 
with analytically computed coefficients, except in case (c), where we multiply the theoretical 
Gaussian decay by a prefactor 0.088 being equal to the theoretical value of the plateau ()4.64|) . 
and re-scale the exponent of the Gaussian by a factor 0.8 taking into account the effect of a 
non-small first term in the exponent of ()4.38() . Note that in the limit S* — > oo the agreement 
with semiclassical theory improves and that the size of the resonant spikes is of the same order 
as the drop in the linear response plateau. The insets show the data and the theory on the 
normal scale. 



namely up to 5 ~ v^ time ^2 can be simply estimated by Trr- We therefore have 

t2 = min{l, :rfcig}Trr = min{ const ft ' 5~ , const fi~ ' 5~ }. (4.67) 



Time scale t2 can be seen in Figure 14.41 as the point of departure of fidelity from the plateau 
value. Using our model with a = 1.1 and /3 = and the position of the initial coherent state 
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at {'d*,ip*) = (1,1) this can be calculated to be t2 = inm{0.57\'S/6,0.57/{5'^vS)} (similarity 
of numerical prefactors is just a coincidence). These theoretical positions of ^2 are shown with 
vertical chain lines in Figure 1^^ and are given by TrrJ^Qi^5/h in Figures l4.5b ,.b while it is t„ for 
a strong perturbation 5S = 3.2 in Figure H^ . Long time decay of fidelity is shown in figure H31 
Theoretical Gaussian decay (|4.66|) . shown with a chain curve, is again confirmed by numerical 
results. The average f given by classical f is ()4.37|1 

8sm^ [aj/2) ^ J 

The derivative f' is 

r =- (4 cot {aj/2) - a Uj + a(l - f) cot {aj/2)\ /sin^ (aj/2)| , (4.69) 

8 L L J J 

which gives using A = 1/(1— j^) for spin coherent states the decay time (|4.6()j) t„ = 0.575~^?i^'^. 
Note that we do not have any fitting parameters, except in the case of a strong perturbation 
{5S ^ 1, Figure ESt) where the prefactor and the exponent of a Gaussian had to be slightly 
adjusted due to the non-negligible effect of the first term in 1)4. 38|) [see caption for details]. 

Average fidelity 

We should remark that, although we obtain asymptotically Gaussian decay of fidelity for a 
single coherent initial state, one may be interested in an effective fidelity averaged with respect 
to phase-space positions of the initial coherent state, similarly as for the case of general per- 
turbations 1)3. 42() . For times t < t2 when we have a plateau, the average fidelity {F{t))- will 
also have a plateau of the same height as the plateau for a random initial state calculated in 
Section [4.2.41 Equation 1)4.82(1 . In the linear response regime this plateau is just twice as large 
as for a single coherent state. Long time decay of the average fidelity {F(t))- will be given by 
the phase space average of the regular decay time t„- The calculation is very similar to the 
one for a general perturbation described at the end of Section 13.2.21 so we wont repeat is here. 
Asymptotically we will have a power law decay with the power determined by the order rj of 
zeros of f' • A ~^f' and dimen s iona lity d (see Eq. I3.43() which is not universal as claimed by 



some authors ([.Tacquod et al. 1 120031^ . To demonstrate asymptotic power law decay (|3.43|) we 

take the same regular kicked top model (|4.51|) as used throughout this section with parameter 

a = 1.1, P = and spin size S = 1000. We take perturbation ()4.57|) of strength 6 = 5- 10"'^ and 

average fidelity over 200 coherent initial states placed randomly over sphere. The results are in 

Figure Uni The t^-^, i.e. f' • A~^r' has for our r' ()4.69() two zeros of the first order, rj = 1, at 

j = ±1. The derivative (i.e. the first nonzero term in the expansion around zero) of the term 

f' ■ A~^f' evaluated at j = — 1 is 

2a — sin a , , ^ . 

c = T— — -. (4.70) 

2sin2(a/2) ^ ' 

The asymptotic theoretical decay ()3.43() should then be 

(m), X ^. (4.71) 

where for a = 1.1 the coefficient is c = 2.4. This theoretical decay agrees with numerics in 
Figure 14.61 The theoretical value of the plateau according to the Equation ()4.84() for a random 
initial state should be 0.58 which also agrees with numerics. 
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Figure 4.6: Average fidelity decay {F{t))- ()4.71|) for the regular kicked top with freezing (no 
fitting parameters!). The dotted line is the theoretical power law decay (|4.71|) . The horizontal 
line gives the theoretical value of the plateau (|4.84|1 . 



Illustration in Terms of Wigner Function 

All the phenomena described theoretically in the preceding subsections can be nicely illustrated 
in terms of the echoed Wigner function — the Wigner function WpU (cos ??, (p) of the echo- 
dynamics. For details regarding the spin Wigner functions see Appendix \K\ and the discussion 
at the end of Section EHl The fidelity F{t) is given simply by the overlap of the echoed Wigner 
function and the Wigner function of the initial state. Therefore, the phase-space chart of the 
echoed Wigner function contains the most detailed information on echo-dynamics and illustrates 
the essential differences between different reg imes of fidelity decay. Thi s is shown in Figure E771 
for the quantized top, see also online movies ( Prosen fc Znidario . l2^^3a^ . In the initial classical 
regime, t < ti, the echoed Wigner function has not yet developed negative values and is in 
point-wise agreement with the Liouville density of the classical echo-dynamics. In the plateau 
regime, ti < t < t2, the echoed Wigner function decomposes into several pieces, one of which 
freezes at the position of the initial packet providing significant and constant overlap — the 
plateau. At very particular values of time, namely at the echo resonances, different pieces of 
the echoed Wigner function somehow magically recombine back into the initial state. In the 
asymptotic, ballistic regime, t > t2, even the frozen piece starts to drift ballisticly away from 
the position of the initial packet, thus explaining a fast Gaussian decay of fidelity. 



4.2.3 Echo Resonances for Coherent Initial States 

Let us now discuss the behaviour of the fidelity for initial wave-packets in the regime of linear 
response approximation in some more detail. We shall consider possible deviations from the time 
average plateau ()4.4U|) . Specifically we will explain the resonances observed e.g. in Figure ETH 
For such a resonance to occur the phases of terms (S (i)) 



(|4.44|1 have to build up in a 
constructive way and this is clearly impossible in a generic case, unless: 
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Figure 4.7: [Movies online ((Prosen fc Znidaric 
echo-dynamics for a quantized top a = 1.1, /3 = 



2003 



Snapshots of the Wigner function of 
0,7 = with S = 200 and for 6 = 1.6 ■ 10"^ 
(same as in Figures I4.4bl4.8b ). The upper hemisphere is shown with j = cosi} £ [0, 1] on the 
vertical axis and 9 = ip € [0, 27r] on the horizontal axis. From top to bottom we show: the 
initial state at t = 0, the state at i = 14 ~ ii when we are around the regular Ehrenfest time, 
at t = 300 in the middle of the plateau, and at t = 100000 in the ballistic regime. 



(i) We have one dimension d = 1, so we sum up over a one-dimensional array of integers n 
in action space'''. 

(ii) The wave-packet is localized over a classically small region of the action space/lattice such 
that a variation of the frequency derivative duj{j)/dj over this region is sufficiently small. 

In this subsection we thus consider a one-dimensional case, d = 1. For simplicity we will 
focus on a linear response regime, so we will study time-dependent terms (S^(t)) and (S(t)) 
with 5](t) given in Eq. ()4.35|) . In the justification of the time average plateau. Section [4.21 we 
replaced a sum over n with an integral over action space. Here the time is not small enough 



*In more than one dimension we would clearly need a strong condition on commensurability of frequency 
derivatives over the entire region of the action lattice where the initial state is supported. 
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Figure 4.8: Structure of echo resonances for coherent initial states of the modified quantum top 
Uo = exp{-iS[a{SJSf/2 + -f{S^/S-j*f/6]), S = 200,6 = 1.6 • lO'^, for increasing values 
of a;" = 7 = (a), uj" = 1 (b), and oj" = 4 (c), which weakens and broadens the resonances. 
Note that in (a), uj" = we have the same data as in Figure H^lb . Vertical chain lines show 
theoretical times tr/2 for vr, and tj- for 27r resonances. 

and furthermore the very graininess of quantum actions is the source of resonances and we will 
have to explicitly consider sums over n. We seek a condition, such that consecutive phases in 
the exponential build up an interference pattern. 



27r-resonance: 

Let us expand the frequency around the centre of the packet 



1 



a;(i) = u;* + (i - ry + -(j - r) v + 



(4.72) 
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where co* = uj{j*), to' = duj{j*)/dj and to" = d'^io{j*)/dj'^. The phases in the sums of the form 
H4.49() with factors e""'^'-^") come into resonance, for m = 1 and hence for any higher m > 1, 
when they change by 27r per quantum number, which happens at time t^, 

2tt 
huj't, = 27r, tr = ^-7, (4.73) 

nuj' 

and its integer multiples. It is interesting to note that these resonant times correspond precisely 
to the condition for revivals of the wave-packet in the forw ard evolution (apart from a phase- 
space translation) studied in ( Braun fc Savichevl . ll99q ^ and ( Leichtle et al. Lll996l l. In addition, 



we need that the coherence of linearly increasing phases is not lost along the size of the wave- 
packet, i.e. 

.. ^ 9 hmuj"t ^ / , „ . X 

C := mu;"f A? = — < 27r, (4.74) 

where we denoted by A^ := ((j — j*)'^) = y^/v the action-width of the initial wave-packet. 
Therefore, if C ^ 27r we will observe echo resonances at integer multiples of 27r-resonance time 

The shape of the echo resonance can also be derived. Let time t be close to fctr, A; G Z, and 
write t = ktr + t' where t' <C tr, so that hco't' <^ 2tt. We can estimate the general time dependent 
sum over n in entirely analogous fashion to Eq. (|4.49|) by: (i) shifting the time variable to t' , (ii) 
incorporating the resonance condition, and (iii) approximating the resulting sum by an integral, 
due to the smallness of t' , 

J2 g{hn)e'"''^'^^''^'clp{hn) ^ ^ 5(7in)e™("*+(^"-^*)"'+5(^«-i*)'-")(*^+*')rf^(;j„) 

n n 

= e'™"^** Y. 5(?in)e^™((^"-^*)'^'*'+^(^"-^'*)'"^"*)(ip(nn) 

n 

= e-*,(/)^_=exp(^ Al^TTe)- ^'-''^ 

We see that the resonance has a Gaussian envelope, modulated with an oscillation frequency 
~ w*. At the resonance centre t' = 0, and assuming ( <^ 1, we can easily calculate the linear 
response terms in the fidelity, getting (S^(t)) = (S(t)) . Therefore, for small ^ the fidelity at a 
27r echo resonance is 1. We get a perfect revival of fidelity. For non-negligible (", the resonance 
height dies as (1 -|-C^)^^ . As C depends on time t, the resonances will decrease in height with 
the increasing order k of the resonance. That is, for the k-th. resonance occurring at t = kt^ we 
have 

( = k--, (4.76) 

and the magnitude of the resonances therefore falls as ~ 1/k. We will get strong and numerous 
resonances, i.e. small C) provided that either the second derivative co" is small, or the initial 
state is squeezed such that A ^ 1. For example, if the second derivative vanishes everywhere, 
uj" = 0, then the resonances may appear even for extended states. This is the case for our 
numerical model, where resonances can be seen also for a random state in Figure I4.1U1 The 
Gaussian envelope of the resonance has an effective width 



A,= '-^^l + er.k/Vh. (4.77) 
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In the semiclassical hmit, the resonance positions scale as tj- (x h^ , while their widths grow 
only as Af oc ?i~ ' , so they are well separated. On the other hand, with increasing order k the 
resonance width grows, and they eventually start to overlap at /c ~ 1/yh. This overlapping 
takes place at time ~ fL~ ' and is smaller than t2 provided 5 < i'qi^ ^• 

The structure of 27r— resonance is nicely illustrated in a numerical example in Figure 14.81 
where we consider a slightly modified model with 



Uo = exp {-■iS[a{S,/SY/2 + ^{S,/S - fY/Q]) 



(4.78) 



corresponding to /io(j) = ctj^/2 + 7(j — 3*Y /Q and a coherent initial state at (??, 93) = (1, 1). 
For such a system uj" = 7 may not be identically vanishing and uj' = a. 
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Figure 4.9: Two-time correlation function C{t', t") 1)2. 18(1 for the quantized top with a = 1.1, 7 = 
0,S = 16, and coherent initial state at {■&, ip) = (1, 1). The structures giving rise to a vr and 27r 
resonance can be seen. 



We note that one may obtain a resonance condition for a single time-dependent term (|4.75j) 
with fixed m even at a shorter time, namely for t = tj./m. This is trivially the case for 
perturbations with many, or at least more than one Fourier components with |?Ti| > 1. However, 
in such cases only selected time-dependent terms of the moment will be affected, so the fidelity 
will generically not come hack to 1, even in the linear response regime and in the strongly 
resonant case C <C 1. Such (incomplete) resonances at fractional times {k/m)t^ will be called 
27r/?TT.-resonances. 

However, we may obtain a resonant condition at t = tr/2 even for the first Fourier com- 
ponent m = 1 of the perturbation, as taking the square of the operator S(t) produces Fourier 
components m + m! = ±2. Such a resonant behaviour at times {k + 2)tj- will be called a 
vr— resonance. 

So for perturbations with a single Fourier mode m = ±1, or more generally with only 
odd-numbered Fourier modes m = 2/ + 1, the tt— resonance can affect only the term having 
gi(m+m. )uj(nn)t ^^ ^^iq expression for the second moment. All other terms (having a single m in 
the phase) result in their time-averaged values. To see this, we observe that the time dependent 
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parts of the form g{hn)e^'^^^ are proportional to X^n '^p(^'^)5'(^"')(~l)™"- ^^ "^ i^ ^'^ ^^^ 
number and dp{hn)g{hn) is a smooth function of n, this sum averages to zero. This allows us to 
again compute explicitly the fidelity in linear response close to the peak in a strongly resonant 
case. We find 

(S(t)) = - E Wm{f)e'^'\ for |t-(A: + ^)t,|«Ai, and C « 1, 

m=2l+l 

(S2(i)) - ( ^ ^„(/)ei-e*)'+( ^ «;„(r)e^™(^*+'^**))', (4.79) 

m=2l+l m=2l+l 

6' ( 

1-F{t) = 72 E \wm{f)\cos{mio*t + f3„ 

"' \m=2l+l 

where (3m are phases of complex numbers Wmij*)^^"^^* ■ So we have learned that the fidelity 
at the peak of a vr-resonance displays an oscillatory pattern, oscillating precisely around the 
plateau value -FSf ()4.62p with an amplitude of oscillations equal to 1 — F9^ so that the fidelity 
comes back to 1 close to the peak of the resonance. 

Again, our numerical example illustrates such an oscillatory structure of tt— resonance in 
Figure 14.81 The resonances can also be nicely seen in 'short-time' Figure 14. 4[ and because 
C = also in the 'long-time' Figure 14.51 In Figure 14.91 we depict the structure of vr— and 
27r— resonance as reflected in the two-time correlation function C{t',t"). Note that the first 
intersection of the soliton-like-trains for t' — t" = const and t + t' = const happens at tj-/2 and 
produces a vr— resonance, while the second intersection at tr produces a 27r-resonance. 

In analogy to the emergence of a tt— resonance as a consequence of the contribution from the 
second moment of S(t), even for the first Fourier mode m = 1, we shall eventually obtain also 
fractional 27r/p-resonances at times {k/p)tj- in the non-linear-response regimes where higher 
moments (SP(t)) contribute to F{t) ~ {exp{—[T,{t)5/h)). This is illustrated numerically in 
FigureEH: showing the case of strong perturbation 5S = 3.2 so that higher orders are important. 
One indeed obtains fractional resonances, some of which have been marked on the figure. 



4.2.4 Random Initial States 

The second specific case of interest is that of random initial states. Here we shall assume that 
our Hilbert space has a finite dimension J\f, like e.g. in the case of the kicked top or a general 
quantum map with a finite classical phase space, or it is determined by some large classically 
invariant region of phase space, e.g. we may consider all states \n) between two energy surfaces 
El < hQ^hn) < E2 of an autonomous system. In any case we have the scaling 

where V is the volume of the relevant classical phase space. Throughout this section we will 
assume that the Hilbert space size M is sufficiently large, so that the difference between the 
expectation value in a single random state and an average over the whole space can be neglected. 
Also, as discussed in Section [2.2.21 the difference between the average fidelity amplitude squared 
and the average fidelity is semiclassically small. 



70 



Chapter 4. Special Case: Zero Time Averaged Perturbation 



1 
0.98 



8S=0.32 




S=1600 



S=200 



10 



100 
t 



1000 



10000 



Figure 4.10: Short time fidelity for a quantized top with [3 = 1.4 (having no singularities) and a 
random initial state. In order to reduce statistical fiuctuations, averaging over 20 random initial 
states is performed for S = 1600, and over 100 initial states for S = 200. The horizontal chain 
line shows the semiclassical theory (|4.82|) . Echo resonances (two spikes for S = 200) are here 
present due to the special property u)"{j) = and will be absent for a more generic unperturbed 
system. The main figure shows the case of the weak perturbation 5S = 0.32, whereas the inset 
is for the strong perturbation 5S = 3.2. 



The Plateau 

To calculate the plateau we replace the quantum expectation value (•) with a classical phase 
space average, resulting in a fidelity amplitude 

2 
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RIS ~ 
plat 



(2vr)" 
V 



/d^, 



d'^6' 

(27r)^ 



■ exp 



w{3,0] 



(4.81) 



Interestingly, the plateau for a random initial state is just the average plateau for a coherent 
state squared, where averaging is done over the position of the initial coherent state. If we 
denote the plateau for a coherent initial state centred at j* by F^\^{3*) (Eq- 14.611) . then the 
plateau for a random initial state F^^ is simply 

{2^Y 



pRIS ~ 
-^plat — 



V 



■ jd^JF, 



plat \3 ) 



(4.82) 



Similarly as in mixing dynamics there is a square relationship between the plateau for RIS and 
CIS. In the mixing case the plateau for CIS does not depend on the position of the initial state 
due to ergodicity, here though the plateau for CIS does depend on j* and one has to take a 
square of the average plateau for CIS. In the linear response approximation one has a simple 
formula for the plateau 
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(4.83) 



TTIT^O 



4.2. Regular Dynamics 



71 



1 



0.8 



0.6 



0.4 



0.2 




10 



100 
t 



1000 



10000 



Figure 4.11: Short time fidelity decay for /3 = (we have singularities) and a random initial 
state. The chain line shows the theoretical value of the plateau as computed from Eq. (|4.84)) . 



In Section IT^ we already discussed possible divergence problems in Wm- If we express Wm in 
terms of Vm, Eq- (|4.36|1 . we have denominators sin (m • u>{j)/2). We therefore have divergences 
at points in phase space where m ■ u{j) = 2'Kk with k an integer. For a coherent initial state 
this was not a real problem as it would occur only if we placed initial packet at such a point. 
For random initial state though, there is an average over the whole action space in the plateau 
formula and if there is a single diverging point somewhere in the phase space it will cause 
divergence. The solution is very simple as explained in Section [4.21 The integral is actually an 
approximation for a sum over fin and so we have to retain the original sum over the eigenvalues 
of the action operator and exclude possible diverging terms. The formula for the plateau in the 
case of such divergences is 






1 

Jf 



E 



pCIS 
^plat 



{fin) 



(4.84) 



where in the summation over n we exclude all terms for which for any constituent Fourier mode 
m we would have m • u{fin) = Ink. 

Again we find an excellent confirmation of our theoretical predictions in the numerical 
experiment for the same system as for a coherent initial state (|4.51|) . In the first calculation 
we choose the shift (3 = 1.4 so that we have no singular frequency in the action space. In 
Figure H. 101 we demonstrate the plateau, which in the case of random states starts earlier than 
for coherent states, namely at ti oc fi 6^ (|4.46|) . The value of the plateau can be immediately 
written in terms of the result for coherent states, 
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(4.85) 



The integral has to be calculated numerically. Horizontal chain lines in Figure H.IOI correspond 
to this theoretical values and agree with numerics, both for the weak perturbation 5S = 0.32 and 
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the strong perturbation 6S = 3.2 (inset). Small echo resonances visible in the figures are due to 
the fact that the Hamiltonian is a quadratic function of the action and therefore uj" = 0, so that 
the resonance condition (|4.7.S)) is satisfied also for extended states (|4.74|) . For a more generic 
Hamiltonian these narrow resonant spikes will be absent. In Figure H. Ill we also demonstrate 
the plateau in the fidelity for the zero-shift case /? = with a singular- frequency, uj{j = 0) = 0, 
where we again find an excellent agreement with the theoretical prediction ()4.84|) . In this case 
the theoretical value has been obtained by replacing an integral in (|4.85j) with a sum as implied 
by Eq. 1)4. 84() and summing over all quantum numbers except n = 0. Observe that the value of 
the plateau is much lower than in the case of a non-zero shift /? = 1.4 in Figure H. 101 as quantum 
numbers around n = will be close-to resonant. 



Long Time Decay 
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Figure 4.12: Long time power law fidelity decay for random states in the kicked top with /? = 
and a = 1.1, for S = 200 (a), and S = 1600 (b). Here averaging for S = 200 is performed 
over 1000 initial random states, otherwise all is the same as for Figure 11.101 The heavy chain 
line shows the theoretical asymptotic decay (|4.87|) with an analytically computed prefactor (no 
fitting parameters). The inset in the bottom figure shows the diffractive quotient between the 
numerical fidelity and the asymptotic formula (|4.86|) (chain line in the main figure). 



After a sufficiently long time t > t2 fidelity will start to decay. To calculate this decay we 
have to evaluate the ASl formula ()4.42|) . For a uniform average over Hilbert space we have 
Dp = 1/J\f with the classical limit dp{j) = {2Trh) /V. The fidelity amplitude is therefore 
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The stationary phase procedure is completely analogous to the situation for the case of non- 
residual perturbation, described in Section [3.2.31 Eq. (|3.46|) . We will only write the asymptotic 
result, 

'f \'^ ft 

'ran \ . . "- 



Fit) X (^^j , t,an = const X ^. (4.87) 

Here we should remember that the asymptotic formula (|4.87|) has been obtained as a stationary 
phase approximation of an integral in the limit of an infinite action space. If we have a finite 
region of the action space, as is the case for the kicked top, the stationary phase approximation 
gives an additional oscillating prefactor, whose amplitude is damped as {h/t)^'"^ for 7i — > 
and/or t ^ oo, and which can be interpreted as a diffraction. This oscillating prefactor can be 
seen in numerical data for the fidelity in the inset of Figure 14.121 For random states the time 
t2 when the plateau ends is independent of h and is determined by the ratio of two competing 
terms in the BCH form or the echo operator, t2 ~ \5T,{t) /h\/\6'^t /2h\, therefore we have 

t2 ~ 1/6. (4.88) 

This agrees with the numerical results shown in Figures l4. 101 and 14. Ill 



Chapter 5 

Coupling with the Environment 



When a man tells you that he got rich through hard 
work, ask him: 'Whose?' 

— Don Marquis 



The fidelity might not always be the relevant measure of stability. Coupling with the envi- 
ronment is usually unavoidable so that the evolution of our system is no longer Hamiltonian. 
To preserve the Hamiltonian formulation we have to include the environment in our description. 
We therefore have a "central system", denoted by a subscript "c", and an environment, denoted 
by subscript "e" . The names central system and environment will be used just to denote two 
pieces of a composite system, without any connotation on their properties, dimensionality etc. 
The central system will be that part which is of interest and the environment the rest. The 
Hilbert space is a tensor product 7i = Tic ® Wo and the evolution of a whole system is deter- 
mined by a Hamiltonian or a propagator on the whole Hilbert space 7i of dimension N = A/'cA/'g. 
The unperturbed state \4'{t)) and the perturbed one |V'(5(0) are obtained with propagators U{t) 
and Us{t) (|2.3|) . Fidelity would in this case be the overlap of two wave functions on the whole 
space. But if we are not interested in the environment, this is clearly not the relevant quantity. 
Namely, the fidelity will be low even if the two wave functions are the same on the subspace of 
the central system and differ only on the environmental part. 

5.1 Reduced Fidelity, Purity Fidelity 

We can define a quantity analogous to the fidelity, but which will measure the overlap just on 
the subspace of interest, i.e. on the subspace of the central system. Let us define a reduced 
density matrix of a central subsystem 

Pc{t) := irMt)], pf{t) ■■= tre[/.^(t)], (5.1) 

where tre[»] denotes a trace over the environment and p^{t) = Ms{t)p{0)Ms{t)'^ is the so-called 
echo density matrix. Throughout this chapter we will assume that the initial state is a pure 
product state, i.e. a direct product, 

1^(0)) = iVc(O)) iVe(O)) =: iVc(O); Vc(0)), (5.2) 

where we also introduced a short notation \il>c;ipe) for pure product states. The resulting initial 
density matrix p(0) = |^(O))('0(O)| is of course also pure. The fidelity was defined as F{t) = 
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ti\p( 0)p (t)] (|2.9j) and in a similar fashion we will define a reduced fidelity (jZnidaric fc ProsenL 



4c 

i) 



2003) denoted by FR(t), 

Fnit) := tv,[pMpfit)]. (5.3) 

The reduced fidelity measures the distance between the initial reduced density matrix and the 
reduced density matrix after the echo. Note that our definition of the reduced fidelity agrees 
with the information-theoretic fidelity (|2.1|) on a central subspace Tic only if the initial state is 
a pure product state, so that /Oc(0) is also a pure state. 

One of the weirdest features of quantum mechanics is entanglement. Of some interest is 
therefore also whether the coupling with the environment will produce entanglement and how 
fast. Due to the coupling between the central system and the environment the initial product 
state will evolve after an echo into the pure entangled state Ms(t)\'ijj{0)) and therefore the 
reduced density matrix p^{t) will be a mixed one. For a pure state |V'(^)) the criterion for 
entanglement is very simple. It is quantified by a purity I{t), defined as 

lit) := tr,[p2(t)], p,(t) := tre[|V(t))(V'(t)|]. (5-4) 

Purity, or eq uivaleri t lv von Neumann entropy tr(plnp), is a standard quantity used in decoher- 



ence studies ()ZurekLll99lh . Iff the purity is less than one, / < 1, then the state {tp) is entangled 
(between the environment and the central system), other wise it is a product state. Similarly, 
one can define a purity after an echo, called purity fidelity i Prosen fc Seligmanl . l2002h Fp{t), 



Fp(t) := trc[{p^it)n (5.5) 

All three quantities, the fidelity F(t), the reduced fidelity F^{t) and the purity fidelity Fp{t) 
measure stability to perturbations. If the perturbed evolution is the same as the unperturbed 
one, they are all equal to one, otherwise they are less than one. The fidelity F(t) measures the 
stability of a whole state, the reduced fidelity gives the stability on "He and the purity fidelity 
measures separability of p^{t). One expects that fidelity is the most restrictive quantity of the 
three - p{0) and p (t) must be similar for F{t) to be high. For Fp^{t) to be high, only the 
reduced density matrices Pc(0) and p^{t) must be similar, and finally, the purity fidelity Fp{t) 
is high if only p (t) factorizes. It looks though as fidelity is the strongest criterion for stability. 

5.1.1 Inequality Between Fidelity, Reduced Fidelity and Purity Fidelity 

Actually, one can prove the following inequality for an arbitrarv vure state IV') and an arbitrary 
pure product state |(/)c;0e) ( Znidaric Sz ProsenL 120031 : JProsen et al. Ll2003al ^. 



K<Ac;<Ae|V')l' < \{MPcm\^ < tvM, (5.6) 

where pc := tre[|V')(V'|]- 

Proof. Uhlmann's theorem ( UhlmannL Il97fih . i.e. noncontractivity of the fidelity, states 



that tracing over an arbitrary subsystem can not decrease the fidelity, 

tr[|0e;0e)(0c;0e||V')(V'|] < tr[|0e) (</'c|Pc]. (5.7) 

Then, squaring and applying the Cauchy-Schwartz inequality |tr[^T^]p < tr[^^' ] tr[i?i?T] we 
immediately obtain the wanted inequality 1)5. 6() . 

The rightmost quantity in the inequality / = ti[p'^] is nothing but the purity of state {ip) 
and so does not depend on \(j)c',(pe)- One can think of inequality 1)5. 6|) as giving us a lower 
bound on purity. An interesting question for instance is, which state \(j)c',<Pc) optimises this 
bound for a given \ip), i.e. what is the maximal attainable overlap |(0c; </'o|V')l (fidelity) for a 
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given purity. The rightmost inequahty is optimised if we choose |0c) to be the eigenstate of 
the reduced density matrix pc corresponding to its largest eigenvalue Ai, Pc\<i>c) = ^i\4>c)- To 
optimise the left part of the inequality, we have to choose |(/)e) to be the eigenstate of pe := trc[p] 
corresponding to the same largest eigenvalue \^,Pp\(i)p) = Xi\4)e)- The two reduced matrices pe 



and pc have the same eigenvalues ( Araki fc LiebL 1197(1 ^. Ai > A2 > • • • > A^v;- Foi' such choice 



of \(j)c',4>c) the left inequality is actually an equality, [(c/ic; i;^e|'0)l = l('?^c|Pc|</'c)|^ = Af and the 
right inequality is 

\l<tT[pl]=Y,\l (5.8) 

i=i 
with equality iff Ai = 1. In the case when the largest eigenvalue is close to one, Ai = 1 — e, 
the purity will be / = (1 — e)^ + O(e^) ~ 1 — 2e and the difference between the purity and the 
overlap will be of the second order in e, / — |((/>c; i?!)e|V')l ~ ^^ ■ Therefore, for high purity the 
optimal choice of \(j)c] V^s) gives a sharp lower bound, i.e. its deviation from / is of second order 
in the deviation of I from unity. 

For our purpose of studying stability to perturbations, a special case of the general inequal- 
ity (|5.6j) is especially interesting. Namely, taking for 1-0) the state after the echo evolution 
^^5(^)1-0(0)) and for a product state |0c;</'e) the initial state |V'(0)) (|5.2|) . we obtain 

F^{t)<Fl{t)<F^{t). (5.9) 

Immediate consequence of this inequality is that if the fidelity is high, the reduced fidelity and 
the purity fidelity will also be high. In the case of perturbations with a zero time average in 
Chapter |1J the fidelity freezes at the plateau and from the inequality we immediately know that 
the same phenomenon will be present for the reduced fidelity and the purity fidelity. 

5.1.2 Uncoupled Unperturbed Dynamics 

Special, but very important case is when the unperturbed dynamics Uq represents two uncoupled 
systems, so that we have 

Uo = Uc® Ue. (5.10) 

This is a frequent situation if the coupling with the environment is "unwanted", so that our 
ideal evolution Uq would be an uncoupled one. The reduced fidelity -FR(t) and the purity fidelity 
Fp{t) have especially nice forms in such case. 
The reduced fidelity (|5..S|) can be rewritten as 

FK{t) = tr,[pMP^it)] = trc[Pc{t)pM, (5.11) 

where Pc{t) is the unperturbed state of the central system and /jf(t) := tVe[Us{t)p{0)UJ{t)] 
the corresponding state obtained by perturbed evolution. Whereas for a general unperturbed 
evolution the reduced fidelity was an overlap of the initial state with an echo state, for a 
factorized unperturbed evolution it can also be interpreted as the overlap of the (reduced) 
unperturbed state at time t with a perturbed state at time t, similarly as the fidelity. 

The purity fidelity can also be simplified for uncoupled unperturbed evolution. As the Uq 
is in factorized form, we can bring it out of the innermost trace in the definition of the purity 
fidelity and use the cyclic property of the trace, finally arriving at 

Fp{t) = tr,[{p^{t)}'] = tiMm'] = m. (5.12) 

The purity fidelity is therefore equal to the purity of the forward evolution. The general in- 
equality gives in this case 

FHt) < Flit) < I{t), (5.13) 
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and so the fidehty and the reduced fideUty give a lower bound on the decay of purity. Because 
the purity is frequently used in studies of decoherence this connection is especially appealing. 

In all our theoretical derivations regarding the purity fidelity we will assume a general 
unperturbed evolution, but one should keep in mind that the results immediately carry over 
to purity in the case of uncoupled unperturbed dynamics. Also a large part of our numerical 
demonstration in next two sections will be done on systems with an uncoupled unperturbed 
dynamics as this is usually the more interesting case. 

5.1.3 Linear Response Expansion 

We proceed with the linear response expansion of the reduced fidelity (|5.3j) and the purity 
fidelity 1)5. 5() . We will use notation p^ := |V'c(0))(V'c(0)| for initial pure density matrix on a 
central system and pe := lV'e(0))(V'e(0)| for the environment. The perturbed propagator is 
defined in exactly the same way as for the fidelity ()2.6|) in terms of the perturbation generator 
V . To order 0{5'^) we get, 

l-F(t) = (^ (S(t)(l®l-pe®/>e)S(t)) 

5-^ 



1-Fn{t) = ^-j (S(t)(l-pe)®lS(t)) 

l-Fp(i) = 2('0 (S(t)(l-pc)®(l-/Oc)S(t)), (5.14) 

where (•) = tr[{pc Pe)*] denotes the quantum expectation value in the initial product state 
and S(t) is the sum of V{t) (|2.12|) . If the expectation values are written explicitly in terms of 
expectations in the base states \j; u), j = 1, . . . ,Mc, u = 1, . . . , A^, with the convention that the 
first base state |1; 1) := [■i/'c! V^e) is the initial state, we have 

l-F{t) = (^0'{(l;l|S2(t)|l;l)-(l;l|S(t)|l;l)2} 
l-FR(t) = (0'|(l;l|S2(t)|l;l)-^|(l;^|S(i)|l;l)|4 (5.15) 

1 - Mt) = ^{i)\ (1; i|s'(t)|i; 1) - E 1(1; HS(t)|i; i>|' - E 10; M^m; i)\' 

As one can see, the linear response expansion of course also satisfies the general inequality (|5.9|) . 
The difference between -FR(t) and F{t) as well as between Fp{t) and F(t) is in off-diagonal ma- 
trix elements of operator S(t). Somehow reminiscent perturbative expansion, although without 
time dependence, has be en obtained in studying the eigenvalues of the reduced density ma- 
trix ( Kiibler Sz Zelj . ll973l ). Depending on the growth of linear response terms with time we will 



again have two general categories, that of mixing dynamics and that of regular dynamics. 

Mixing Dynamics 

For mixing dynamics the correlations decay and the linear response term will grow linearly with 
time. For large times one can argue that S(i) should look like a random matrix and the terms 
giving the difference between the fidelity and the purity fidelity and the reduced fidelity can be 
estimated as 

Ef=2l(i;i|s(0|i;i)l' E,IF(t)]o;i),(i;i)P j_ 

(i;i|s2(t)|i;i) '^E,,JP(i)](i;i),a,.)P^A4' ^- ' 
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because there are more terms in the sum for fidelity. Therefore we can estimate the difference 
Fp (t) - F2(t) ~ 1/Mc + l/AAe and F-R,{t) - F{t) ~ l/Mc- Provided both dimensions A^ce are large 
and for sufficiently long times, so the "memory" of the initial state is lost and the correlations 
decay, we can expect the decay of all three quantities to be the same. 

Regular Dynamics 

For regular dynamics on the other hand, S(t) will not approach a random matrix but will grow 
with time with a well defined long time limit E(t) — > Vt. This will happen for times larger 
than the averaging time tave (|3.17|) . Expectation value of 5]^(t) will then grow quadratically 
with time. In a similar way as we defined the average correlation function C (|3.15j) . we can also 
define a time average of the correlation functions occurring in the linear response expressions 
for F^{t) and Fp(t), 

Ch := U.n <^""^-g>^^^"» =(V-(.-Pc)0lF) 

In the linear response regime purity fidelity and reduced fidelity will therefore decay quadrati- 
cally with time, if Cr and C-p are nonzero. 

For coherent initial states the average correlation function C is proportional to fi 1)3. 36(1 
provided v{j) is sufficiently smooth. In the semiclassical limit we can make an expansion 
around the centre of the wave packet v{j) ~ v{j*) + v' [j*){j — j*} + • • •. The second moment 
(^'^)ci — (^')ci is then proportional to the dispersion in j of the initial packet (|3.33j) . i.e. to fi. As 
we take a product initial coherent packet we have a dc + de dimensional squeezing matrix A 

and we can write derivatives of classical v as 

^=:iX.K). ^;:=^. f^^-^. (^-l^) 

where we split actions j into two components j =: (jcjc)- The formula for C (|3.3H|) can then 
be written as 

c = h [KK'K + KK'K) , (5.20) 

By a similar method one can calculate also Cr, (|5.17() and gets 

CK = ln[v',A-X)- (5.21) 

For Cp on the other hand, the lowest order expansion of v' used for C and Cr gives zero. To 
get Cp we must expand v to the second order in {j — j*) 

vij) ^ vif) + v'{f){j - f} + i{j - f} . v"{f){j - j*} + • • • , (5.22) 

where v"{j*) is a symmetric matrix of second derivatives evaluated at the position of initial 
packet j* , 
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.;;.(.-) :^ ^. (5.23) 



80 Chapter 5. Coupling with the Environment 

The only nonzero contribution to Cp comes from off-diagonal terms t;^^ and v'l^ = (t'ce) ; with 
the final result being 

C^ = (\n\ trM, n := A- VeAc~M'c- (5-24) 



2 

The result for Cp is very interesting as it means that the purity fidelity in regular systems will 
decay as Fp{t) = 2{5/TCf'C-pt^ on % -independent time scale because of C-p oc % . Note that to 
reach this conclusion of 7l-independent decay, we need only the existence of a smooth classical 
limit of 5](t), therefore F-p[t) will decay on ?i-independent time scale even before V converges, 
i.e. for t < tavg- 

A simple special case of time averaged perturbation V is the tensor product form V = Vc^dVe- 
Then we have v'^^, = v'^®v'^/v (note that in order to comply with our previous notation we have 
v'^ = (dcVc)ve)- If in addition Ac,e are diagonal, then the average Cp is 

Cp = [^^hj {<KX) {v'eK'<) ■ (5.25) 

As we can see, for a tensor product forms of the time averaged perturbation there is a relation 
between C, Cr and Cp, 

v^Cp = Ck{C-Ck). (5.26) 

The purity fidelity will be high, i.e. decoherence will be slow, if we either make Cr, small or 
make Cr close to C. Note that if Cr, is small, then the inequality (|5.9|) already tells that Cp 
will be also small, but here we have a stronger result, that Cp is of second order in h whereas 
Cr is only of first order. 

5.1.4 Numerical Illustration 

Here we would like to briefly demonstrate quadratic decay of the purity fidelity for regular 
systems and linear decay for chaotic systems. Furthermore, we will show that one can also have 
a situation where the purity fidelity Fp{t) decays slower for chaotic systems than for regular. 
We take a Jaynes-Cummings model (Section I2.1.2p which is a two degrees of freedom model of 
a harmonic oscillator coupled to a spin. As initial state we always take a product state of two 
coherent states, for the spin the coherent state 1)2. 26(1 is centred at (i9*, (/?*) = (1, 1) and for the 



boson coherent state (|2.3fl|) we take a = 1.15. We take two different sets of parameters. For the 
chaotic regime we choose to = e = 0.3 and G = C = 1, for which the energy of the initial state 
is E = 1.0 and the classical Poincare section shows a single ergodic component. The second set 
of parameters is lo = e = 0.3 and G = 1, while C = 0, for which the systems is integrable and 
the energy is E = 0.63. The perturbation is chosen to be in the spin energy e, i.e. a detuning, 
corresponding to the perturbation generator V 

V = hS^. (5.27) 

Because the simulation is very time consuming for the chaotic case, we choose a small spin size 
S = 4, resulting in an effective Planck's constant H = 1/4. 

Note that, as opposed to previous kicked top model, here the system is conservative, and 
the time index t is continuous. First we want to check the growth of linear response terms 
with time. In Figure W/H we show the correlation integral Sp{t) for fidelity and the difference 
S'F(i) — Sp{t) between the correlation integrals for fidelity Spit) and for purity fidelity Sp{t), 

Spit) := ^(S(t)(l®l-/>c®/?e)S(t)) 

Sp{t) := ^(S(t)(l-pe)»(l-/>e)S(t)). (5.28) 
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Note that for chaotic systems 3^(1) should converge to 2a ^'^.'A\f . where a is the transport 
coefficient, i.e. the integral of the correlation function. The difference Sp^t) — Sp{t) should 
on the other hand be of the order ~ 2cr(l/A4 + 1/A4) ~ cr/2 (due to 1/A4 + 1/A4 ^ l/S). 
These two predictions are nicely confirmed by the lower two curves in Figure 15.11 We can see 




Figure 5.1: The correlation integral SF{t) and Sp{t) — Sp{t) (|5.28j) in the Jaynes-Cummings 
model for regular dynamics (upper curves) and chaotic dynamics (lower curves). The horizontal 
dashed lines (for the chaotic case) are at 2a = 0.2 and 0.2/4. The two linearly increasing dashed 
lines for regular dynamics have slopes C = 0.046 and (7(1 — 0.98/4). 

that the correction to the purity fidelity is really of order l/S in the chaotic case and would 
therefore vanish in the semiclassical limit. For regular dynamics we expect Sp{t) to grow linearly 
with time, with the slope given by C, SF{t) -^ Ct, with C tx 1/5 1)5. 20() . The purity fidelity 
integral should grow as Sp{t) — > Cpt with Cp oc 1/S^ (|5.24j) . The difference should therefore 
be SF{t) — Sp{t) ^ 5f(1 — const/5). This is again confirmed in Figure I^TTl We checked that 
5p oc 1/5^ also for larger 5, up to 5 = 24. 

Next, we would like to demonstrate a faster decay of the purity fidelity in the regular 
regime than in the chaotic one. We set the perturbation strength to 5 = 0.005 for which we 
expect the purity fidelity to become lower in the regular regime at the time determined by 
1 - Fp{t) = 26'^S'^Cpf ^ 1 - F{t) = S'^S'^2at, i.e. at t ^ 9, where we used a = 0.1 and 
Cp = 0.011 as given in Figure E^ In Figure \^7^ we show the purity-fidelity for chaotic and 
regular regimes. We can see the crossover at the predicted time t ~ 9. We also illustrate 
the evolution of the purity fidelity with the square of the Wigner function corresponding to the 
reduced density matrix p^{t) on the central spin subspace. If W^m {-d, (p) is the Wigner function, 
purity fidelity can be written as a phase space integral 



Fp{t)= fdnw^Mi^,^^). 



(5.29) 



For details about the spin Wigner function see Appendix Near the top and bottom of 
Figure 15.21 we see this evolution for the chaotic and the integrable Hamiltonian respectively. 
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Figure 5.2: Echo dynamics in the Jaynes-Cummings model. The square of the Wigner function 
for chaotic dynamics (top diagrams) and integrable dynamics (bottom diagrams) is shown as a 
function of time at times corresponding to the axis. The purity fidelity is shown in the main 
frame on the same time scale and for short times in the inset. Red curves give the integrable and 
blue curves the chaotic evolution. In the inset full curves show the complete numerics, symbols 
the evaluation starting from the numerical correlation integrals of Figure E3 and dashed curves 
the linear or quadratic approximation using a and Cp. 



In the centre of the figure we plot the purity fidelity on the same time scale as the Wigner 
functions in the main frame and an amplification of short times in the inset. We observe 
detailed agreement of numerics with the results obtained from the numerical values of the 
correlation integrals 3^(1) and Sp{t), reproducing the oscillatory structure of the decay. From 
the same correlation integrals we obtained the coefficients for the linear and quadratic decay, 
which agree well if we discard the oscillations. It is important to remember that the integral 
over the square of the Wigner function gives the purity fidelity and therefore the fading of the 
picture will be indicative of the purity fidelity decay. On the other hand the movement of the 
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centre is an indication of the rapid decay of fidelity (not shown in the figure) . 



5.2 Mixing Dynamics 



Discussing linear response results (Section l5.1.3|) in the case of mixing dynamics we have shown 
that the linear decay is the same for all three quantities. Similar random matrix arguments as 
for the linear response can be used also for higher order terms and therefore one expects that 
in the semiclassical limit of small l/AAc + l/AAe we will have the same exponential decay ()3.5|) 



Fp{t) ^ Flit) ^ F\t) = exp (-2t/T„ 



(5.30) 



with the decay time Tm = h'^ /26'^ac\ (|3.5() independent of the initial state. This result is expected 
to hold when S(t) can be approximated with a random matrix for large times (|5.16|) and V 
does not contain terms acting on only one subspace. Such terms could cause fidelity to decay 
while having no influence on purity. 
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Figure 5.3: Decay of F^(t),Fp^(t) and I{t) (dotted curves) in the mixing regime of the double 
kicked top. The solid line gives the theoretical exponential decay ()5.3U|) with Tm calculated from 
the classical Uci = 0.056. Horizontal chain lines give the saturation values of the purity and the 
reduced fidelity, 1/200 and 1/400^, respectively. 

For numerical verification of this result we chose a double kicked top system ()2.25() with 
e = 0, so that the unperturbed evolution is uncoupled. For the double kicked top model we will 
always have an uncoupled unperturbed evolution in all our numerical demonstrations in this 
chapter. Other parameters are 7c, e = '7r/2.1 and ac,e = 30, ensuring chaotic classical dynamics. 
The spin size is chosen to be S* = 200, so that we have Mc.e = 25 + 1. The perturbation is 



V 



^Oz yy Ozj 



(5.31) 



with the strength 5 = 8 • 10 ^. The coherent product initial state is placed at t?*g = vr/-v/3, 
Vce = vr/-v/2. We show in Figure 1^31 the decay of the fidelity F{t), the reduced fidelity Fp{t) 
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and the purity I{t). Clean exponential decay is observed in all three cases, on a time scale t^ 
H5.3U() given by the classical transport coefficient a^- We numerically calculated the classical 
correlation function 

Cel(t) = [(Zc(t)^c(O)) J2, (5.32) 

where we took into consideration that the unperturbed dynamics is uncoupled and is the same 
for both subsystems and that (z)^[ = 0. Taking only the first term Cci(O) = 1/9 would give 
o"ci = 1/18 (|3..S|) while the full sum of Cci(i) gives a slightly larger value cJci = 0.056. Exponential 
decay, of course, persists only up to the saturation value determined by a finite Hilbert space 



size (see Section l2.2.1jl 



5.3 Regular Dynamics 

5.3.1 Beyond Linear Response 

Purity Fidelity 

For regular dynamics and coherent initial states one can calculate purity fidelity to all orders 
in the semiclassical limit and not just the linear response expansion 1)5.24(1 . For times larger 
than the averaging time iave the echo operator goes towards Ms{t) — > exp{—i6Vt/h). In the 
semiclassical limit we can use the classical limit v instead of V. Classical v{J) is a function of 
action operators only and so similarly to the evaluation of fidelity one can use the ASI (|3.31|) 
for the evaluation of purity fidelity. A partial trace over the environment gives an integral 
over je and due to a square in the definition of purity fidelity we end up with an integral over 
2{dc + de) = 2d dimensions, if dc is the dimension of jc and de of je, 



Fp{t)^n-^'' /dFexp 



.6 ,,. . 



vU'cJe) + vij'^Ji) - Vijcj'e)} 



dp{j)dp{j'), (5.33) 



where dF = d'^jd'^j', j = (jc, Je), j' = (jc> j'e) ^^"^ dp{j) is the classical limit of {n\p\n) ()3.31|) . 
Next we expand v{j) around the position of the initial packet (|5.22j) . The constant term v{j*) 
and the linear terms cancel exactly as well as the diagonal quadratic terms, regardless of the 
position of the initial packet. The argument of the exponential function which remains is then 



.6t 



^T[{Jc-j'c)-v'L{f)Ue-j'e] 



h 



(5.34) 



with the matrix v" of second derivatives given in Equation ()5.23() . The resulting expansion 
can be used in -Fp(t) (|5.33|) to calculate the purity fidelity for initial states having well defined 
classical dp{j). For coherent initial states dp{j) are Gaussian ()3.33|) and so the whole integral 
is also a Gaussian, giving 

Mt) = ^=, (5.35) 

det {1 - -i5tL-W') 



where L and V" are 2d dimensional matrices 
/Ac 



L:-- 






Vo 





Ac 










Ac 







Ac/ 



V" :-- 













V 



■^ec 
;r." 



-Vc 











(5.36) 



/ 



with matrices of squeezing parameters Ac and Ag ()5.18|1 of dimensions dc and dg, respectively, and 
dc X de dimensional matrix u"g (|5.23|) . Note that the determinant is real despite the imaginary 
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unit. Tlie determinant of tlie matrix 1 — i6tL^^V" can be simplified using the following identity 
for block matrices 

det (^ ^ ^ ) = ^^* *^^) ^"^^ (^ ~ CA-^B) = det (D) det {A - BD-^C), (5.37) 

with m X n dimensional matrices B and C^. Noting that the matrix 1 — iStL^^V" has 1 on 
the diagonal, we obtain 

det (1 - [6tL-^V") = det (1 + {6t)^Z), (5.38) 

with Z being the 2dc x 2dc dimensional matrix 

"^ •" 2 I -A- V' A- V. A-V.A-1^:' )• ^^■'^^> 



"ce e ^ec c "ce e ^ec 



To simplify the determinant of 1 + {5t) Z we use the following identity 

detf'^ ^^ = det {DA- BC), (5.40) 

for square matrices A,B,C and D and commuting B and D, [B,D] = 0. Using this finally 
gives 

Note that u is the same dc x dc matrix we had in the expression for Cp 1)5. 24() . This very simple, 
yet important result deserves a little discussion. It gives the purity fidelity (and thereby as a 
special case also the purity) decay for regular systems and coherent initial states to all orders in 
6. Its validity is limited to times larger than the averaging time tave in which V converges and 
for sufficiently small h. Planck's constant must be small to allow the replacement of V with its 
classical limit v and v must be smooth on the scale of the wave packet size \/h. Furthermore, 
we replaced a sum over the quantum numbers with an integral over the actions, the validity 
of this being given by the condition Stlv'^Jfi^ /h <C 1, i.e. (5tu < l/h, if tu is the upper time 
limit of validity of the ASI. As we will see, the asymptotic decay of purity fidelity as given by 
Equation 1)5. 41() is between l/{6t) and l/{6t)'^'', depending on the matrix u. The Fp(t) at the 
upper border t^ is therefore between h and H ". The purity fidelity will for long times saturate 
at the plateau given by the finite Hilbert space size, Fp ~ l/Mc ~ ^ ''. Comparing this, we see 
that the upper border tu coincides with the point where the asymptotic saturation ^p is reached 
for one degree of freedom systems. For dc > 1 the ASI break time tu is only by a constant factor 
smaller than the time when we reach Fp. Summarising, the purity fidelity decay ()5.41() is valid 
from tave all the way to the asymptotic plateau -Fp provided h is sufficiently small, without any 
bound on 5. 

Let us now explore the asymptotic decay of Fp(t). To simplify theoretical arguments, we 
assume that the squeezing parameters are all equal to one*, L = 1 so that we have u = "Uce^cJ- 
In 1 + dg degrees of freedom systems (dc = 1), the matrix u is just a number, and the purity 
fidelity decays as 

^'« = TttW' * = '■ "■''' 

A single parameter u is already fixed by the linear response, i.e. by the value of Cp (|5.24|) . 
Asymptotically we get Fp{t) x l/{St) decay regardless of the second dimension d^. For general 

* Usually this can be achieved by the right choice of actions. 
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systems with dc > 1 one can see that the determinant det{l + {6t)^ u} is a polynomial of 
order at most dc in (5t)^. Furthermore, as the matrix u is symmetric and positive definite, its 
eigenvalues are positive, meaning that all the coefficients of the polynomial are positive and the 
Fp{t) is always less than one. In a special case, when the matrix u can be written as a tensor 



the dimensions involved and we get 



product of two vectors, like u = I'ce^ce ='■ x ® y, the determinant is again simple regardless of 



^p(^) = -7T=7 ux2 u = x®y. (5.43) 

In such case we again get the asymptotic -Fp(t) ^ ^/{^i) decay, but here regardless of both 
dimensions. Matrix u has such a form for instance if all the matrix elements of v"g are the same 
(e.g. for the perturbation v = J2k i iJc (^ jc)ki)- Our result of course a pplies als o to th e purity 
decay in weakly coupled systems and does not agree with the result of l.Iacquodl (|2nr)4l ) . 

Reduced Fidelity 

For the reduced fidelity a similar procedure based on the ASI can be used. In the expansion 
of v{j) around j* (|5.22() the first non vanishing term is a linear one. The resulting Gaussian 
integrals are analogous to the ones in the calculation of the fidelity decay for coherent initial 
states (but there are twice as many), the final result being a Gaussian decay 



FR(t) = exp(-^CRt^l, (5.44) 

with a single parameter Cr = 2^ (v'^A~^v'^) given by the linear response alone (|5.2H1 . 

Numerical Illustration 

To numerically confirm the above formulas for the decay of purity fidelity and reduced fidelity 
we take a double kicked top model H2.25() . To be in a regular regime we take ac,e = and 
7c = 7r/2.1 and 7e = vr/v7. Different unperturbed frequencies 7c 7^ 7e are chosen in order to 
have a general situation, i.e. that the subsystems are not in resonance. Unperturbed dynamics 
is again uncoupled, e = 0, so that the purity fidelity equals purity. The spin size is chosen 
S = 100 and the initial product coherent state is placed at {{}*,ip*)c = 7r(l/-v/3, l/-v/2) and 
(??*, (p*)c = 7r(l/\/3; 3/\/7)- The perturbation with strength 6 = 0.01 is of the form 

V = j^S', Si (5.45) 

Unperturbed classical evolution is very simple, namely rotation around y-axes for angles 7c and 
7e for the central system and the environment, respectively. The classical limit of V is readily 
calculated and expressed in terms of the actions jc = Vc and jc = ye, 

^ = i(l-j2)(i_j-2). (5.46) 

The derivatives of v are v'f. = — jc(l — ie)/^ and similarly for v'^. The squeezing parameter for 
the spin coherent states is A = 1/(1 — j^) and the average correlation function C (|5.20)) is 

C = ^[i-£)(i-f^){fji-£)+f^(i-f^)}, (5.47) 
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Figure 5.4: Decay of F{t),F^{t) and I{t) in the regular regime of the double kicked top. For 
parameters see text. The theoretical Gaussian decay for the fidelity (|3.36j) and the reduced 
fidelity ()5.44|) , with the theoretical value of C 1)5. 47() and Cr 1)5. 48() , overlap with the numerics 
(symbols) within the line width. Theoretical purity decay (|5.42|) (dotted curve) as determined 
by the theoretical u (|5.49j) also nicely agrees with the numerics (full curve) up to the plateau. 



which for our choice of the initial coherent packet gives C = 0.024/5". For the reduced fidelity 
we get 

CK = lj!{l-f,){l-j!), (5.48) 



which evaluates to Cr = 0.021/S'. The above two parameters C and Cr completely determine 
the Gaussian fidelity (|3.36() and reduced fidelity 1)5. 44() decays. As we have 1 + 1 degrees of 
freedom system, the purity is determined by a single parameter u (|5.42j) . For our perturbation 
we get 

^ = JcJe(l-Jc)(l-Jc), (5.49) 

giving u = 0.032 for our initial product coherent state. As our perturbation V is of the product 
form, we could as well use a universal relation between Cp, Cr and C 1)5. 26() to calculate Cp 
(or equivalently u) . The numerical results are shown in Figure 15.41 together with the theory. 
Agreement is excellent. 



5.3.2 The Jaynes-Cummings Model 

In this section we will consider in detail the stability of a Jaynes-Cummings model, described in 
Section r2.1.2| under various perturbations. The Jaynes-Cummings system can be experimentally 
realized and so it is a possible model on which one could experimentally study the quantum 
stability. In experiments one usually has only a co-rotating term, i.e. G' is zero, therefore 
we will focus on a situation when we have C = in the unperturbed dynamics and so the 
classical limit of the unperturbed system is integrable. For the initial state we will always 
choose a product of coherent states for a harmonic oscillator, given by a real parameter a 
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(|2.3()|) . and a spin, given by the initial position {"d,^) (|2.2H|) . In studies of the reduced fidelity 
and purity fidelity we will consider the spin as a central system, and the harmonic oscillator 
as an environment. For regular dynamics with coherent initial states the decay of fidelity and 
reduced fidelity is Gaussian with the decay time determined by a single parameter C and Cr, 
while the decay of the purity fidelity is Fp{t) = l/-\/l + u {5t)^ (|5.42j) and is again determined 
by a single parameter u. With this in view, it is sufficient to determine only the linear response 
parameters C (|5.20|) . Cr (|5.2H) and Cp (|5.24|) . to get the decay to all orders in 6. Therefore, 
we will focus only on the calculation of these coefficients. The calculation of the expectation 
values in the coherent initial state is described in Appendix IbI 
Our unperturbed system will be 

Ho = hiva+a + heS^ + G^={aS+ + a+5„). (5.50) 

v25' 



For the perturbation we will look at four different situations: 

• variation of to, corresponding to y = ha~^a =: [5uj] 

• variation of e, corresponding to V = hS^ ='■ [de] 

• variation of G, corresponding to ^ = ^=^(a5+ + a'^S-) =: [dO] 



variation of C, corresponding to ^ = -^{aS- +0+5+) =: [SG'], 



where we introduced short notation V = [6»] symbolically denoting the perturbation in the 
parameter •. In the last case of ^ = [SG'] and if G = all correlations average to zero, 
corresponding to the residual perturbation (i.e. freezing), and this case will be considered 
separately in Section 1^31 

As an additional simplifying assumption, we will assume that we have w = e, i.e. spin 
and an oscillator are in resonance. For G 7^ we have numerically checked that the average 
correlation functions are quite insensitive to the resonance condition. The only exception being 
the V = [5G] perturbation, for which C, (7r and Cp vanish if we are out of resonance and 
have G = 0. Because this freezing case will be studied separately, we do not loose generality 
by taking uj = e but gain in the simpler theoretical calculations. Namely, in resonance the 
G term commutes with the unperturbed Hamiltonian, [aS+ + a^S-,HQ] = 0. In addition, in 
experimental implementations we are usually close to the resonance condition. 

Perturbation of the frequency lo 

We will first analyse the case when the perturbation consists of the variation of uj, so that the 
perturbation generator is 

V = ha+a. (5.51) 

Due to the resonance condition uj = e the unperturbed propagator Uq factorises into a part 
proportional to the G term and the rest. The perturbation in the Heisenberg picture is therefore 



To calculate the average correlation function one needs an integral of the perturbation S(t) = 
Jo y{'T)<^T in the limit t — > 00. As the constant G just re-scales time in the perturbation V{t) 
1)5.52(1 . the G's which are an infinite time limit property of S(t) will not depend on G, except 
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Figure 5.5: Numerically calculated dependence of C's on the position of the spin coherent state 
(/? and cos ■(9 for uj = e = 0.3, G = \ (everything is independent of G provided it is nonzero), 
G' = 0, a = 1.15 and J = 4. Left column shows the results for e perturbation and the right for 
the perturbation in uj. 



possibly at G = where the symmetry of K^ changes and we might have an effect due to 
degeneracies. Direct calculation in this simple case of G = gives 



G 



1^ 

5 S ' 



Gr = Gp = 0, (ifG = 0), 



(5.53) 



for the initial product coherent state (|2.26|2.3()j) . Note that o? j S has a well defined classical 
limit, namely the energy of the oscillator. The last two results are expected. The Gp = because 
the perturbed dynamics is also uncoupled and Gr = because the perturbation is only in the 
oscillator ("environment") part of Hq. For G > we get an additional angle dependent term 
of order fi = 1/ S"^ which can be seen in numerically calculated G's in Figure 15.51 Therefore, 
there is a discontinuous jump in all three G's at G = 0, the discontinuity being proportional to 
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h and thus of higher order for C while it is the leading order term for Cr and Cp . Otherwise 
the values of the plateaus C are independent of the coupling G, but of course the time scale on 
which we get the convergence of C"s scales as 1/G. For small G the discontinuity will happen 
at large times. 

Perturbation of the spin energy e 

In this case the perturbation generator is 

V = hS^. (5.54) 

Everything is analogous to the previous case and in the simple uncoupled case of G = (and 
cj = e) we get 

- — 1 sin ^9 - 

^ = ^R = ;5^-' C'p = 0, (ifG = 0). (5.55) 

The C now does not depend on the oscillator coherent state parameter a but instead depends on 
the position of the spin coherent state. Also, the Cr is now nonzero as we make the perturbation 
in our central system. At G = there is again a discontinuity of order h in all three C"s, 
otherwise they are independent of G. For G and Cr this discontinuity is of higher order in h 
and can be neglected in the semiclassical limit. Dependence of numerically calculated C's on 
the initial position of the spin coherent packet is shown in Figure 1^31 One can see, that Cp is 
equal for [5uj] and [6e] perturbations. Therefore the decoherence (purity fidelity) is insensitive 
to which frequency we detune, whereas fidelity and reduced fidelity are not. 

Perturbation in the coupling G 

The perturbation generator is 

V = -^{a+S^+aS+). (5.56) 



We again take the resonant condition to = e. Now the perturbation is constant in time V{t) = 
V{0) regardless of the value of G and we can calculate all three C's for an arbitrary G. These 
are written in Figure 15.61 next to each plot showing its dependence on the initial spin packet 
position. There is no dependence on G and also no discontinuity at G = 0, contrary to previous 
two examples. Also the C and Cp are strictly larger than in the case of u; or e perturbation 
whereas the Gr is of order h and is smaller than for the e perturbation. Now Gp oc ^ is of 
course nonzero and gives us the decay time scale of purity fidelity (or purity in the case of 
G = 0). 

5.4 Separation of Time Scales 

Until now we have discussed two broad categories, mixing systems in which the whole correlation 
function decays and regular systems where we had a plateau in the correlation function. In this 
section we will consider situation where the time scale of the environment is much smaller than 
that of the central system and the correlation function can be considerably simplified. We will 
furthermore consider perturbations of the product form 

V = Vc0 Ve. (5.57) 

For numerical illustrations we will use a coupled double kicked top with e = 0, i.e. an uncoupled 
unperturbed system. To simplify notation we will use {Ac) to denote the average of operator Ac 



5.4. Separation of Time Scales 



91 



V = [SG] = ^{a+S. + aS+) 



C = Ck + ^ sin^ ^ 



0.15 

0.1 

0.05 







- - 2 / \ ^ 







^P — "02 cos 77 




Figure 5.6: Dependence of C"s on the position of tlie spin colierent initial state {(p and cost?) 
for uj = e = 0.3, G = 1 (everything is independent of G), G' = 0, a = 1.15 and J = 4. The 
left column shows exact theoretical results whereas the right one shows numerical plots of this 
dependence. All is for the perturbation in the coupling constant G. 



acting only on the central system in the initial state of the central system, {Ac) = trc[AcPc{0)], 
and similarly for the environment. If the operator acts on the whole system, then (A) denotes 
the expectation in the whole initial state, as before. When the time scale of the environmental 
correlations {Vc{t)Vc{t')) is much smaller than the time scale of the central systems' correlations 
{Vc{t)Vc{t')) , time averaging over the fast environmental part of the perturbation can be per- 
formed. Regarding the environmental correlation function two extreme situations are possible. 
If the correlations of the environment decay, we will call such a case "fast mixing environment" , 
and we have a finite integral of the environmental correlation function. If the correlations of the 
environment do not decay, we will call it a "fast regular environment" , and we have generically 
a non vanishing average correlation function of the environment. 
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5.4.1 Fast Mixing Environment 

The situation, when the time scale te on which the correlation function for the environment 
decays is much smaller than the time scale ic of the central system, is of considerable physical 
interest. This includes various "brownian" like baths, where the correlation times are smaller 
than the dynamical times of the central system in question. The correlation sums Sp^t), Sp{t) 
(t^:^ and 5R(t), 

Snit) := ^(S(i)(l - Pc) ® lS(t)), (5.58) 

giving the linear response decay of the fidelity F{t) = 1 — {6/h)'^tSF{t), the reduced fidelity 
F^{t) = l-{5/h)HSR{t) and of the purity (or generally purity fidelity) I{t) = l-2{5/n)HSp{t) 
can be significantly simplified in such situation. We will furthermore assume (14) = 0, with 
(A) = hmt^aot"^ Jq i^iO) dC denoting a time average. This assumption corresponds to an 
equilibrium situation where the average "force" Ve vanishes. The integration over the fast 
variable Vg can be carried out and we get for t ^ tc ^ tg 



i-e 



t^oo 



SF{t) = 


2cTe(K') 


Sn{t) = 


2ae{(y,2)-(y,)2| 


Sp{t) = 


2ae{(y,2)-(y,)2|, 


Se'(i))/2t, 


with Se(t) = / Ve{C)d^ 

Jo 



(5.59) 
with 



(5.60) 



being the integral of the autocorrelation function for the environmental part of the perturbation 
Ve alone. The result does not depend on the initial state of the environment. Some interesting 
conclusions can be drawn from these linear response results (|5.59|) . 

We can see that the decay time scale depends only on the time average diagonal correla- 
tions of the central system {V^{t)) and not on the full correlation function. This is a simple 
consequence of the separation of time scales and means that the decay of all three stability 
measures does not depend on the dynamics of the central system (e.g. being mixing or regular). 
Furthermore, reduced fidelity -FR(t) and the purity I{t) will decay on the same time scale (|5.59|) . 
meaning that the decay of the reduced fidelity is predominantly caused by the loss of coherence, 
i.e. entanglement between the two factor spaces. This in turn means that the reduced fidelity, 
which is a property of echo dynamics, i.e. of comparison of two slightly different Hamiltonian 
evolutions, is equivalent to the decay of the purity. 

If the initial state of the central system Pc(0) is a Gaussian wave packet (coherent state) 
then the dispersion (V^) — {Vc) is by factor of order h smaller than (V^'^). Thus for coherent 
initial states of the central system, irrespective of the initial state of the environment, the i^R(t) 
and I(t) are going to decay on a 1/h times longer time scale than F{t). We have therefore 
reached a general conclusion based on the assumption of chaotic fast environment (which is 
often the case), namely that the coherent states are most robust against decoherence (decay 
of purity), provided te *^ ^c and decoherence time is longer than the correlation time of the 
environment, tg <C tdec- If decoherence is even faster than the time scale of the environment, as 
is the case for macroscopic superpositions, then formulas ()5.59|) are not valid any more as one 
is effectively in a regular regime with Sp{t) oc t^. Decoherence time is then in dependent not 



just of systems dyna mics but also of environmenta l dynamics characterised by ag ()Braun et al. , 
|2Qa3;EmnZ£lil[2003; Stninz fc Haakfl. IF"" 



In the regime of fast chaotic env ironment one can derive a master equa.tion fo r a reduced 
density matrix of the central system ( Kolovskvl . 1 1994 : iMevstre &: Sargent IIJ . |l990|) . We take a 
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Figure 5.7: Various correlation sums from formulas H5.59() in a fast chaotic regime (solid curves, 
as indicated in the figure). Chain lines indicate corresponding theoretical time averages. For 
details see text. 



partial trace over the environment of the echo density matrix p (t) and write it for a small time 
step Ai. This time step Ai must be larger than the correlation time ie of the environment and at 
the same time smaller than the correlation time t^ of the central system. For the environmental 
part of the correlation function we assume fast exponential decay (particular exponential form 
is not essential) which is independent of the state p(0) 



ii,[v,{t)vs')p\ 



— exp{- 

6p 



|t-t'|Ac}trcM. 



(5.61) 



Assuming the perturbation to be a product V{t) = Vc{t) (8i Ve{t) and the average "force" 
tre[T4(t)/9] to vanish together with the exponential decay of environmental correlations of the 
form (|5.6H) for an arbitrary state, yields a master equation for the reduced density matrix 
it) := tre[p^(t)], 

{t) = -^ae[V,{t),[V,{t),p^m- (5.62) 



P^ 



P^' 



h' 



Master equation is strictly valid in the Markovian limit of the correlation function being a delta 
function in time or if considered on time scales larger than the correlation time. 

For numerical demonstration we use a double kicked top with: Vc^e = S^/S, S = 200, 
6 = 1.5 ■ 10~^, coherent initial state at {'d,(p)c,e = (7r/\/3, 7r/\/2) and parameters Oc = 0, 
7c = 7r/50 for the central system and Oe = 30, 7e = 7r/2.1 for the environment. Actually, we 
could take any value of Oc and would get qualitatively similar results. The only advantage of 
using regular central dynamics Qc = is that it is then possible to explicitly calculate averages 
{V^) and ly^ - {Vcf. Namely, if Oc = and 7c <1 we get 



iK') 



1 



:i 



yc) + ^(i + yc) 



{Vc') - {VcY 



45 



(1 + y; 



(5.63) 
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Figure 5.8: Decay of F'^{t),F^{t) and I{t) for fast chaotic environment. The dashed line is an 
exponential function with the exponent given by the values of o"e and (V^'^) (|5.59j) while the two 
dotted lines have slopes —2 and —1. For parameters see text. 



The values of these two quantities for our initial condition are (Vc) = 0.202 and (V^^) — (Vc) = 
0.399/5, and are shown in Figure 15.71 with two dotted lines (by pure coincidence we have 
o"e = 0.20 « (V^)), together with numerically calculated time dependent (not yet averaged) 
{V^{t)) and {V^{t)) — {Vc{t)) . This time dependent values oscillate on a time scale ^ 50, 
which is much longer than the time ~ 10 in which Ue converges and so the assumption t^ <C tc is 
justified. The values of all three quantities are then used in the linear response formulas (|5.59|) 
to give us the time scales on which F, Fpt and / decay. The results are shown in Figure 15.81 
We can see that the fidelity again decays exponentially as predicted, but the reduced fidelity 
and the purity have a power-law tails. The decay time can be estimated by the lowest order 
expansions ()5.59() . Using the values of o"e and (V^) — (Vc) (theoretical expression 1)5. 63() ). we 
get rp ~ 1/(0.165^5'^) and tr^p ^ 1/(0.325^5') (rp and tr are for F^ and F^) which agrees 
with numerics in Figure ESI The same general conclussion again holds: the more chaotic the 
environment is (smaller Ue), the slower the decay of all three quantities. Purity and reduced 
fidelity both decay on a 1/h longer time scale than the fidelity in accordance with expressions 
H5.63() for coherent initial states. 



5.4.2 Fast Regular Environment 

Here we will explore perhaps a less physical situation of a regular environmental dynamics, i.e. 
one with non-decaying correlation function. The double integral of environmental correlations 
grows as oc t^ and we can define the average correlation function 



Ce:= lim(S,^(t))/t^. 

t— »oo 



(5.64) 



If in addition the correlations of the system also do not decay then the correlation sum of 
the total system will grow as oc i^ which is just the regular regime already discussed before. 
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Figure 5.9: Times r at which F'^ (t) , F^{t) , I (t) fall to level 0.37 for different 5 and fast chaotic 
environment. Symbols give the numerics and lines give the theoretical dependence of r. All is 
for S = 100. For other parameters see text. 



Here, we will focus on a different situation where the integral of systems correlation function 
converges, i.e. the dynamics of the central system is mixing. We will additionally assume the 
average "position" Vq to be zero (V^) = 0. The transport coefficient of the central system CTc is 
then 

(Tc := lim (S2(t))/2t, Se(t) = / V^iOdC- (5.65) 

The expressions for SF{t), Sp{t) 1)5. 28(1 and S'R(f) (|5.58|) are independent of time and can be 
simplified to 

Sf = 2acCe 
5r = 2acCe 
Sp = 2ac{c,-JV^^]. (5.66) 

Note that now the reduced fidelity F^{t) decays on the same time scale as the fidelity F{t). 
This must be contrasted to the case of a fast mixing environment 1)5. 59() . where F^{t) decayed 
on the same time scale as the purity. If the initial state of the environment Pe{0) is a coherent 
state, then purity will decay on a 1/h times longer time scale than fidelity or reduced fidelity. 
On the other hand, for a random initial state of the environment, the average force (Ve) = 
vanishes, and all three quantities decay on the same time scale. 

For the purpose of numerical experiment we chose Vc = Sz/S, and Vc = S^/ S"^ in order to 
have a less trivial situation of non- vanishing average force (14). The initial condition is again 
(t?, (/7)c,e = (7i"/\/3) 7r/\/2) and the parameters are J = 200, Uc = 30, 7c = 7r/7 and Ue = 0, 
7e = 7r/2.1 and the perturbation strength (5 = 6- 10~^. By choosing the explicitly solvable case 

_ 2 

CKe = we can calculate Ce and (T4) , say for the simple case of a 7r/2 rotation, 7c = vr/2. 



C, 



\ii - ylf + }g{-^yt + ^vl + 1) + oiMS^ 
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Figure 5.10: Correlation sums occurring in 1)5. 66(1 (solid curves) for S = 200. The top chain line 
gives the best fit for o"c and the two lower chain lines give theoretical time averaged correlation 
functions for the environment ()5.67|1 . All is for a fast regular environment. See text for details. 



Ce - (Ve) 



h^'^' 



yl) 



1 



+ Y^(llye-ll2/e+2)+0(l/5' 



(5.67) 



For our parameters we have yc,e = 0.772 giving Ce = 0.0407 and C^ — (14) = 0.120/S'. The 
values of these coefficients are shown in Figure 15.101 (lower two dotted lines) and nicely agree 
with the numerics. In Figure 15. lll we can observe the exponential decay of fidelity and reduced 
fidelity on the same time scale (both curves almost overlap) and the decay of the purity on a 
1/fi longer time scale. For longer times purity decay is again algebraic. In Figure [5.12l we show 
dependence of the decay times on 5. The dependence for purity is quite interesting. If one looks 
at the time the purity falls to 0.99 one has agreement with linear response (by definition). But 
if one looks at the purity level 0.37, they don't agree as well, meaning that the shape of purity 
decay may change (not only the scale) as one varies 5 or fi. On the other hand, this may also 
be simply a finite size effect. 

5.5 Freeze in a Harmonic Oscillator 



Previously we have discussed the so called freeze of fidelity in regular systems having a non- 
singular derivative of the unperturbed frequencies, O = du/dj ^ 0. In the present section we 
will consider the case of a harmonic oscillator, for which fi = and the theory explained in 
Section 14.21 can not be used. For the sake of numerical demonstration we will use a Jaynes- 
Cummings model (see Section [2.1.2(1 with the perturbation in the G' parameter, i.e. the only 
case not explained in Section 15.3.21 The initial state will be a product of coherent state for a 
spin and an oscillator. The perturbation will be 



V 



h 



(a+S+ + a5_) = [6G']. 



(5.68) 
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Figure 5.11: Decay of F'^{t), F^^{t) and I{t) for fast regular environment. The dashed line is the 
exponential function with the exponent given by a product of cJc and Ce (|5.66j) . The straight 
dotted line has slope —1. See text for details. 



Provided we have G' = for the unperturbed system the correlation function has a zero time 
average regardless of other parameters. This is a simple consequence of the symmetries as 
the perturbation is a "counter-rotating" term, while the unperturbed Hamiltonian has only 
a "co-rotating" term. Due to the symmetry the perturbation is residual. The fact that the 
perturbation is residual has nothing to do with the perturbation G' breaking integrability. For 
instance, if we make a perturbation in G' and any other parameter at the same time, we will 
have a perturbation that breaks integrability but is not residual. 

Because we want to study the case of r^ = 0, we choose G = so the unperturbed Hamilto- 



nian, 



Hq = Tiujo^a + heSz, 



(5.69) 



is uncoupled. As a consequence purity fidelity equals purity. Residual perturbation can be 
written as a time derivative of another operator W 1)4. 6|) , which is in our case 



-4i««'-i-^' 



The echo operator can than be written as 



W 



[h 



2S{uj + e) 



{a~^ S+ — aS-) 



Ms{t) = exp \ -^ ({W{t) - W{0)}6 + ^r{t)6^ + 



(5.70) 



(5.71) 



with T(t) given by 



r(t) := j\{r)dT-^[W{0),Wit)], 



R 



^i-.^> 



i^ + e 



(S2 - S'i -S,- 2a+a5,) . 



(5.72) 
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Figure 5.12: Times r at which F'^{t) (pluses), F^^{t) (circles) and I{t) (triangles) fall to level 
0.37, and times r when I{t) falls to 0.99 (squares) for varying 6 and fast regular environment. 
Symbols are the numerics and lines give theoretical dependences of r. Everything is for S = 100. 

For the semiclassical calculations we will also need the classical limits of W and R. Introducing 
canonical action angle variables for the spin (i.e. "central") system. 



Tt-Sz -^ jc = cos ■&, 
and for the oscillator ("environment") 



ns± 



2% 



^2je e^ 



-i6lc 



l-J2e±'^^ 



I |2 /c 
Je = \a\ /S, 



^, 



the classical Hamiltonian reads 

ho = u- j, u = (e, w), j = (jcJe)- 

The classical limits w of W and r oi R are 



(5.73) 



(5.74) 



(5.75) 



w 



cos (^c + ' 



V2je(l-Jc^) 

UO + £ 

r = ^(l-j,^-2jej- 

5.5.1 The Plateau 

The height of the plateau is given by the expectation value of the echo operator 

i^piat = I (exp {-\5{W{t) - Wm/fi)) |2. 



(5.76) 



(5.77) 



For a harmonic oscillator one has to evaluate this expression explicitly and can not use time 
averaging as has been done for a general unperturbed dynamics in Section ^2 We immediately 
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Figure 5.13: The plateau 1 — F{t) (upper symbols) and 1 — I{t) (lower symbols) for S* = 50 and 
coherent initial state. Perturbation of strength 5 = 0.01 is in G' , lo = e = 0.3, G = 0, a = 1.15. 
Chain lines represent the theory 1)5. 8715. 88|1 and circles full numerical simulations. 



see, that the leading semiclassical order will vanish for coherent initial states. We would get 
I exp (— i5s(jf*, 0*, t)/?i)p = 1, where we denoted 



sij,e,t):=w{j,e + ut)-w{j,e). 



(5.78) 



Recall that in case of il 7^ this leading order already gave a nonzero contribution. Here 
though, we have to calculate the next order. For coherent initial states this is easily done using 
the stationary phase method, i.e. expanding the phase around the position of the packet to 
lowest order. The calculation is actually very similar to the calculation of the fidelity decay for 
coherent initial state and non-residual perturbations. Expansion of the phase gives 



sij,e,t) = sif,e*,t) + s' ■ ij-j*,e- e*) + 



with 



( 



ds{j*,9*,t) dsij*,e*,t). 



(5.79) 



(5.80) 



dj ' de 

Note that now s{j,9,t) depends also on the angles and therefore the ASI method cannot be 
used directly. Instead, we will use the classical averaging over the initial Gaussian distribution 
in a phase space. We will replace the quantum expectation value with the classical average. 
Using the compact notation x := {j,0) = {jcjc^c^c), we can write the classical density 
corresponding to the coherent initial state as 



p{x;x*) 



— VdeW exp{-{x-x*)-D{x-x*)/n), 
ttTl J 



(5.81) 



with x* being the position of the initial packet and a matrix D (of size 2d x 2d) determining 
the squeezing of the initial packet m d = dc + de degrees of freedom system. The above classical 
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density is normalised as Jpdx = 1. For our choice of the initial state being a product of a 
spin coherent state H2.28|) and an oscillator coherent state H2.31() the matrix D is diagonal with 
elements 



D 



11 



1 - J^ 



D. 



22 



1 

2^ 



D 



33 



1-Jc, 



D. 



44 



2je. 



The fidelity plateau is now calculated as 



(5.82) 



plat 



/ dxp{x; 0) exp {—\5s' ■ x/h) 



(5.83) 



resulting in a Gaussian function 



Fpi'^ = exp 



-t'ha 



l^ha 



-s'-D-'s'. 



(5.84) 



Note that the result is formally very similar to the expression for the fidelity decay for a non- 
residual perturbation (|3.,36j) . Linear response expansion of the plateau is of course F9^ = 

1 — -^i^har and SO the plateau for a harmonic oscillator is by factor 1/h higher than for general 
systems (|4.62|) . Beyond linear response though, the plateau will decay faster (i.e. as a Gaussian) 
for a harmonic oscillator than for systems with fi 7^ where the asymptotic decay of the plateau 
was a power law (e.g. 1/{5S) for the numerical model used in Section [4. 2|) . 
For our Jaynes-Cummings model we have 



s{j,0,t) = 2 sm I — - — t I ■ ^— sm 



Lo + e 



+ 0^ + t 



U) + £ 



(5.85) 



The coefficient fhar is then 



t'ha 



{u + e)- 



2 1^ + ^. 
■ sm 1 



:i-jc)+2jc 



jl + {I - jl) cos" [e, + e, + t 



Lo + e 



(5.86) 

The same t^ar would be obtained from the quantum calculation of the second moment of S(i). 
Results of the numerical calculation are shown in Figure 15.131 The initial packet was at 
(??, c/?) = (1,1) and a = 1.15, giving actions jc = 0.54 and je = 0.026 and angles ^c = 1 and 
9e = vr/2. Because je <C jc we can neglect the second term in i/har and get 



jiCIS 



1 - Fp^^^ = 3.93r5sin^ (0.3t). 



(5.87) 



The agreement of this theoretical prediction with numerics can be seen in Figure 15.131 

For the reduced fidelity and the purity one can go trough similar calculations. We will just 
list the linear response result obtained from the quantum expectation values, 



1-/ 



5^S 



(w + e) 



■smM^^t 2jc 



ic +(l-Jc)cos 



^'-c + ^e + t^ 



(w + e)' 



sm 



u + e 



i)(l-Jc)' 



(5.88) 



Theory for the purity shown in Figure IB. 131 also agrees well with numerics. The purity plateau 
for coherent initial states is 7l- independent, just as the whole decay of purity. 
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5.5.2 Beyond the Plateau 

After a sufficiently long time, when the second term r(i) becomes important in the echo operator 
H5.71() . the decay will be determined by the operator R. As our R is non-residual, we can use 
the theory for general perturbations, just using a "renormalised" perturbation strength 5^/2. 
The fidelity and the reduced fidelity will decay as Gaussians, with the decay times given by C 
H5.2U() and Cr (|5.21|) . respectively, while the decay of the purity in our 1 + 1 degrees of freedom 
system will be I{t) = 1/^1 + u{6t)'^ (|5.42jl . with u = ACp/h"^ (|5.24|) . For our simple example, 
the average f is equal to r and the derivatives occurring in the expressions for C's are 

df 2 df 2 (9^f 2 

ojc w + e dje uj + e ojcOje oj + e 

Taking into account that the squeezing parameters for coherent initial states are Ac = 1/(1 — j^) 
and for an oscillator Ag = l/(2je), we obtain 

^ = ^ (^ + g)2 {(Jc + Je)'(l - ic ) + 2jgje} 
^P = ^' (^^g)2 (l-Jc)ie- (5.90) 

This then immediately gives the fidelity decay (|4.()(i|) . the reduced fidelity (|5.44j) and the purity 
(E311), 



Fit) = exp i-d'^S'Ct'/A), FR(t) = exp (-5^5^C7Rt74), I{t) 



VT+CpSW^ 

(5.91) 
Note the ^-independent decay of purity. 

In Figure 15.141 we show numerical results for the Jaynes-Cummings model with spin size 
S = 1/h = 50 and a coherent initial state placed at {'&,(p) = (1, 1) and a = 1.15, giving actions 
jc = 0.54 and je = 0.026. The perturbation is 1^ = [6G'] of strength 6 = 0.01, while the 
parameters of the unperturbed Hamiltonian are uj = e = 0.3, G = 0. The average correlation 
functions (|5.90|) are CS = 1.35, CrS* = 1.26 and CpS^ = 0.10. For clarity we show in the figure 
only the fidelity and the purity as the reduced fidelity would almost overlap with the fidelity. 
This theoretical values are then used to compare with numerics. Agreement is excellent. We 
also show a numerical calculation of the classical fidelity (symbols). One can see that for a 
harmonic oscillator and a residual perturbation the classical fidelity agrees with the quantum 
one. Recall, that for a residual perturbation and a general unperturbed system, i.e. having 
rj 7^ 0, the classical fidelity followed the quantum only up to ti oc \/T/h time. For harmonic 
oscillators though, the quantum fidelity is equal to the classical despite the perturbation being 
residual. Purity can also be calculated purely classically. The classical purity can be defined as 

Id := Jdx.plit), pc{t) := Jdxepit), (5.92) 

with p{t) being the classical density at time t, i.e. p{t) = p{(f)~^{x);x*) with cf){x) being 
the Ha miltonian flow in phase space^^. S uch classical purity (linear entropy) has been used 
before ( Wehrll . 119781 : lAneelo et al. I . l2004h . For our system and perturbation the integrals in 



(|5.92|) are Gaussian and the result is the same as the one obtained from the quantum deflnition 
of purity (|^?nT|) . 

^Here the density has to be square normahsed, Jdxp^ = 1. 
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Figure 5.14: Long time fidelity and purity decay for the same data as in Figure 15. 131 The theo- 
retical decay for purity (chain curve) (|5.91|) . agrees with the numerical I{t) until the asymptotic 
plateau is reached. Similarly, theory for the fidelity overlaps with the numerics. Symbols show 
the numerical result for classical fidelity which in this case agrees with the quantum fidelity. 

5.6 Decoherence for Cat States 

In this section we want to study decoherence for macroscopic superpositions of states, the so- 
called Schrodinger cat states. We would like to demonstrate the accelerated decoherence for 
macroscopic superpositions, without resorting to any effective master equation description. The 
goal is to show that if we start with the initial state of the central system in the superposition 
of two macroscopically separated states |ti) and \t2), the decoherence (decay of purity) is much 
faster than if we start with a single state |ri) of the central system. 

We can explain the accelerated purity decay for cat states using our results for the reduced 
fidelity decay for coherent initial states. Decoherence of a cat state will cause the reduced 
density matrix to evolve from a coherent superposition of two packets at the beginning to a 
mixture of two packets after the decoherence time idee , 



1 



Pc(0) ~ - (|Tl)(ri| + |t2)(t2| + \ti){t2\ + |T2)(ri| 



tdc 



1 



(|ri)(ri| + |T2)(r2|) 



(5.93) 



The purity I{t) of the initial coherent state is /(O) = 1 while the purity of the final state is 
-^(idcc) = 1/2, if the states |ti) and \t2) are orthogonal. 

Let us suppose that the initial state is a product state of the coherent state of the environ- 
ment \a) and a cat state of the central system, i.e. a superposition of two coherent states |ri) 
and \t2), 

|V'(O)) = ^(|Ti) + |r2))0|a). (5.94) 



^/2 



We will assume that we are in a regular regime, for which the decay of purity of individual prod- 
uct coherent states has been derived before 1)5. 41|) and was seen do decay on a 7i-independent 
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time scale id ~ l/f^i where 5 is a coupling strength between the central system and the environ- 
ment. For times smaller than the decoherence time of individual coherent states id; the initial 
cat state (|5.94|) will evolve into a superposition of product states, 

iV'(i)) ~ Ixi) 55 |/3i) + 1X2) ® 1/32), (5.95) 

where we used the notation |xi) <8) |/3i) ~ C/(i)|ri) (8) |a) and |x2) '8 |/92) ~ U{t)Rc\Ti) |a), 
where for the second state we have written \t2) =: -Rc|''"i) with some unitary matrix Re. Note 
that the propagator U{t) has a coupling of strength S between the two subsystems. The second 
state can therefore be thought of to be obtained from the same initial state |ri) ® |a) under a 
slightly different, perturbed evolution Us{t) := U{t)Rc- The coherences in the reduced matrix 
(off diagonal elements in p^) of the evolved state (|5.95j) are proportional to the overlap (/3i|/32). 
The square of the overlap (/3i|/32) is in fact the overlap of two reduced density matrices, one 
obtained with the unperturbed evolution U{t) and the other with the perturbed one Us{t). 
Because the final states |/3i,2) as well as the perturbed evolution Us(t) depend on the initial 
states of the central system \ti^2), the square of the overlap is nothing but the reduced fidelity ^ 
of the environment § . The decay of coherences and therefore of the purity will be given by the 
reduced fidelity of a product coherent initial state |ri) |q). As we have seen before (|5.44|) . the 
Cr for a coherent state in a regular regime is Cr = ^fi{v'^ ■ A~^v'^) giving the reduced fidelity 
decay exp (— /^CRi^/^ ), which in turn determines the decay of coherences and therefore also 
purity decay. Here the perturbation strength / depends on the unitary matrix R^ i.e. on 
the "distance" between the initial states |ti^2) of the central system, as well as on U{t). This 
perturbation strength / is of course different than the coupling 6 of the evolution U{t). For 
instance, in the simple case of \t2) being a space shifted packet |ri), \t2) = e~'^'^/^|ri), the 
perturbation strength is proportional to the distance x between the two states. As the coupling 
6 causes decoherence of individual packets on a time scale t^ ~ 1/6 while the decoherence of 
the cat state happens on the decay time scale of the reduced fidelity tdec ~ Vh/l, it is sufficient 
to have / > 6\'h in order to obtain faster decoherence for cat states. The separation / of the 
two packets constituing the cat state has to scale only as / ~ Vh and therefore does not have 
to be macroscopically large in the semiclasical limit. For an environment with many degrees of 
freedom de, and if derivatives ■Dg are all nonzero, Cr will be proportional to d^. The reduced 
fidelity for the central system on the other hand will have Cr proportional to d^ and so if 
de S> dc ) the reduced fidelity of the environment will decay faster than the reduced fidelity of 
the central system. Therefore, at the decoherence time tdec we will still have |xi,2) ~ ki,2), 
justifying the decoherence scenario (|5.93|) . For cat states we therefore have Cp oc h instead of 
Cp (X h a.s for a single coherent state, resulting in an accelerated decoherence. 

Of course, the above argumentation is by no means a strict proof but is just a simple 
illustration o f how an accelerated decoh erence arises. An argumentation using the fidelity has 
been used bv JKarkuszewski et al. I ( 20021 ^. but without realizing the key role played by the time 



scales of purity and reduced fidelity decays. To repeat, several ingredients were needed: (i) the 
dynamics was assumed to be regular, (ii) for the reduced fidelity of product coherent states we 
have Cr oc h, (iii) the purity of constituent product coherent states decays on much longer h- 
independent time scale (i.e. Cp oc h'^) and approximation 1)5. 95() was possible, (iv) furthermore, 
if de ^ dc the resulting decohered states \x) are still approximately equal to the initial states 
\t) of the central system as the reduced fidelity decay is faster for the environment than for 
the central system. All this together enabled us to "derive" a faster decoherence decay for 
macroscopic cat states. 

■'■Actually we used a slightly different definition of the reduced fidelity as an overlap of the initial and the echo 
state. 

''In general the states |/3i,2) can not be obtained by a Hamiltonian evolution on the environmental subspace 
alone. 
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For chaotic dynamics and times longer than the Ehrenfest time ~ ln?i we have seen that 
fideUty, reduced fideUty and purity all decay on the same time scale irrespective of the initial 
state. A natural question arises: How does decoherence of macroscopic superposition in chaotic 
system behave? Actually, decoherence for macroscopic super positions is so fast, that it happens 
on tiine scales shorter than an v dynamical scale of the system (|Braun et al. 1 . 120011 : IStrunz et al. . 
20031 : IStrunz &: Haakd . l2003i ) . It therefore happens on an "instantaneous" time scale on which 
every system behaves as a regular one (i.e. no decay of correlation functions). The chaoticity 
therefore has no influence on decoherence of macroscopic superpositions. 



Numerical Example 
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Figure 5.15: Dependence of Cp (Eq. I5.99( exact expression) for the cat state (|5.96|) on the 
position angles of two constituent coherent packets. We have take ipi = ip2 = 1 and S* = 4 (left 
figure) and S = 64 (right figure). 

We take a Jaynes-Cummings model with the unperturbed Hamiltonian Hq = hcoa'a + ehSz 
and the perturbation in the coupling G of strenght 5, V = [dG]. Purity fidelity is then equal to 
the purity which will be used as a measure of decoherence. The initial state will be a product 
state of a coherent state |a) for the oscillator acting as an environment and the superposition 
of two coherent packets for the spin, 



|V'(0))cOc|Ti)+|r2). 

The Heisenberg evolution is simple, so that the integral of perturbation S(t) is 



AV25 



sin (l^] {e-"^/2a+5_ + e^^/^a^+j , 



(5.96) 



(5.97) 



where A := e — w is a detuning. From now on we set A = in order to have a non vanishing Cp. 
As we take a coherent oscillator initial state, the expectations over oscillator space in correlation 
function sum (|5.17j) can be explicitly made and one ends up with the very neat result 



Cp = ^"^('5+5'- 



{S+){S-)} 



(5.98) 



where now the expectation values are only over a spin state, e.g. the cat state (|5.96|) or any 
other. A similar result h as been obtained before trough a master equ ation approach to deco- 
herence, see for instance ( Benedict fc CziriakLll999l : lFoldi et al. I . l200lh . Our result though has 
quite different perspective as it was derived directly from the Hamiltonian dynamics without 



5.6. Decoherence for Cat States 



105 




20 40 60 80 100 120 140 

t 

Figure 5.16: The purity decay for a random spin initial state (fastest decaying curve), 
Schrodinger cat state (middle curve) and coherent state (slowest decaying curve) . The oscillator 
initial state is always a coherent state with a = 1.15. The parameters of the Jaynes-Cummings 
model are lv = e = 0.3, G = G' = and the perturbation is in G (i.e. in the coupling) with 
6 = 0.02, S = 20. The cat state consists of two coherent SU(2) packets at {'&,(p) = (1,1) and 
(0, 1). The theoretical decay times r = h/6\/2Cp are r = 46, 26 and 13 for the coherent state, 
the cat state and the random state, respectively. Cp for the coherent state has been calculated 
according to the theoretical formula in Figure EH for the cat state according to (|5.99|) and for 
a random state using (|5.100|) . 



any resort to a Markovian description of the reduced dynamics or to dissipation. The calculation 
of Cp (|5.98|) for a cat state (|5.96j) using formulas from Appendix El is straightforward although 
tedious, with the exact but considerably involved result. Rather than showing the whole result, 
we concentrate on the leading order in 1/5* (i.e. h) which will dictate the decay in the semiclas- 
sical limit and therefore will also govern the decoherence of macroscopic superpositions. The 
result is 



Cp 



85 



+ o{n'^ 



V 



sin 7?i + sin ■i?2 — 2sint?i sim?2Cos ((/?i — 932)- 



(5.99) 



The / in the above expression has an interesting geometrical interpretation: it is a distance 
between the projections of positions of the two coherent states to the x-y plane. There is 
indeed a faster decoherence for cat states than for a single coherent state. For a cat state 
Cp ccfi = 1/ S whereas we had Cp (xTi for a single coherent packet. The exact dependence of 
Cp on the position of both packets in a cat state can be seen in Figure 15.151 The accelerated 
decoherence can be seen as the increasing of Cp away from the diagonal ("!?! 7^ 'Q2) for 5 = 64, 
where the exact result for Cp can be approximated by the leading order expression (|5.99|) . 
Note that due to symmetry there is a slow decoherence for states with ■(?! + t?2 = tt despite of 
the fact that this leads to a macroscopic superpositions. Linear response decay time of purity 
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Figure 5.17: Dependence of numerically determined decay times of purity on S (= 1/h) for 
a random spin initial state, a Schrodinger cat state and a coherent state. The perturbation 
strength is (5 = 0.01 while all other parameters are the same as for Figure 15.161 Theoretical 
decay times r = h/6\/2Cp are shown with full lines. 



for macroscopic superpositions is therefore 1)5. 99() tdec = ^^/^/{Sl), while the decay time for 
individual coherent packets is (using Cp from Figure [5^ t^ = l/(-v/2cos^ {'d/2)5). We get an 
accelerated decoherence, tdec -^ t^, for / S> vS^cos^ {'^l'^)- 

For random initial state of a spin and a coherent packet for an oscillator we can calculate 
Cp using the expression (|5.98|) and get 



Cp 



S + l 



random spin i.e.. 



(5.100) 



Again we have Cp oc ?i, just like for cat states, but with a smaller prefactor. If we had a random 
state on the whole Hilbert space, i.e. also for an oscillator, we would have Cp oc li . Therefore, 
the decoherence rate for a cat state is between the decoherence rate of a single packet and the 
decoherence rate for a random state. 

The accelerated decoherence can be seen in Figure l^. IHI where we show numerically calculated 
purity decays for all three different initial spin states. Revivals of purity on a classical time scale 
^ 1/5 due to integrability of unperturbed Hamiltonian can also be seen. This revivals will be 
absent in a general non-integrable Hamiltonian. To further illustrate the scaling of decay times 
with S we show in Figure IB. 171 numericallv calculated decay times in which the purity decays 
to level 0.9 for all three different initial states, coherent, cat state and random state of a spin, 
while the oscillator is initially always in the same coherent state. 



Chapter 6 

Application: Quantum Computation 



In theory, there is no difference between theory and 
practice. In practice, there is no relationship between 
theory and practice. 

— Anonymous 



Quantum information theory i s relatively recent, for a review see ( Nielsen fc Chuand . 12001 : 



Steand . 119981 : lEkert fc Joszal . ll99(3 ). Its beginnings go back to '80 and in recent years theoretical 



concepts have been demonstrated in experiments. While quantum cryptography, a method of 
provably secure communication, is already commercially available, quantum computation on 
the other hand is still limited to small laboratory experiments. 

In order to perform quantum computation you obviously need a quantum computer. A 
quantum computer can be considered as a many-body system of n elementary two-level quantum 
systems — called qubits. The union of all n qubits is called a quantum register |r). The size of 
the Hilbert space M and therefore the number of different states of a register grows with the 
number of qubits as AA = 2". Quantum computation then consists of the following steps: 

• Load the data for the quantum computation in the initial state of the quantum register 
|r), resulting in a general superposition of AA basis states. 

• Then perform the actual computation, represented by a unitary transformation U. As 
U acts on an exponentially large space it is usually decomposed into simpler units U = 
Ut ■■ -112111. Such sequence of T elementary one-qubit and two-qubit quantum gates Ut, 
t = 1,2, . . . ,T is called a quantum algorithm. 

• Finally, we read out the result of the computation by performing measurements on the 
qubits of the final register state U\r). 

A quantum algorithm is called efficient if the number of needed elementary gates T grows with at 
most polynomial rate in n = log2 M, and only in this case it can generally be expected to outper- 
form the best classical algorithm. At present only few efficient quantum algorithms are k nown , 
perhaps the most genera llv useful being the Quantum Fourier Transformation (QFT) (|Shoii 
19941 : ICoppersmithL Il994l ^ which will also be the subject of our study. 



There are two major obstac les for performing pra ctical quantum computation. First, there 
is the problem of decoherence ( Chuang et al. V |l995|) resulting from an unavoidable generally 



time- dependent perturbation due to the coupling between the qubits and the environment. If 
the perturbation couples only a small number of qubits at a time then such errors can be 
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eliminated at the e x pense of extra qubits by quantum error correcting codes (ISteane , \l^ 
Calderbank fc Shoii Il99(il : 



iig, ^^, for another approach see Jxia. X. Lwi El^ 



Second, even if one knows an efficient error correcting code or assumes that quantum com- 
puter is ideally decoupled from the environment, there will typically exist small unknown or 
uncontrollable residual interaction among the qubits which one may describe by a general 
static perturbation. Therefore, understanding the stability of quantum algorithms with re- 
spec t to various types of perturbations i s an iniportant problem, for some results on this topic 



see 



dMiauel et al L Il996. 1 997: Gea-Banaclochel Il998l Il999l. Ignnnl: ISong fc ShepelvanskvL hVlOA 



Georgeot fc ShenelvauskvL fi^; JRermaTi pi »/. 1 120011. I2nn2allbl: ICelaTdo et nl ibooi) 



6.1 General Framework 

Let us perturbe the i-th quantum gate by a perturbation* generated by Vt, 



U!:=eKp{-[6Vt)Ut. 



(6.1) 



We set h = 1 and use the superscript 6 to denote a perturbed gate and the subscript t to denote 
a discrete time index, i.e. a gate number. We allow for different perturbations Vt at different 
gates. The perturbed algorithm is simply a product of perturbed gates, U = C/j. • • • t/j^. Fidelity 
will again serve as a measure of stability and we have 



F{T) = I {Ms{t)) 



Ms := U{-T)U\T), 



(6.2) 



with U{t) := Uf ■ -Ui and similarly for U (t). For our generally time dependent perturbation 
the echo operator equals to 



Ms 



-iSVT(T) 



,g-i5V2(2)g-i<5yi(l) 



(6.3) 



where Vt{t) := U'^{t)VtU{t) is the perturbation of f-th gate Vt propagated with the unperturbed 
gates U{t) = Uf-Ui. Beware that Vt{t) is time dependent due to two reasons, one is due to 
the interaction picture (time in parentheses) and the second is that the perturbation itself is 
explicitly time dependent, i.e. different perturbation for different gate (time as a subscript). In 
quantum computation one is usually interested in initial states containing a maximum amount 
of information, thereby being close to random states. With this in view, we take the initial 
state average to be a trace over the whole Hilbert space. 



(•) 



^tr(.). 



(6.4) 



By this prescription we study the average fidelity over the whole Hilbert space. Without sacrific- 
ing generality we furthermore assume the average perturbation to be traceless, J2t=i ^r^ = 0. 
Trace of the average perturbation only changes the phase of the fidelity amplitude and has 
therefore no influence on the fidelity (i.e. probability). The linear response expansion of the 
fidelity then reads 



F{T) = l-5^ J2 C{t,t') 



t,t'=i 



with the correlation function 



C{t,t'):=j^tv[Vt{t)V,{t% 



(6.5) 



(6.6) 



'It is expected that the computation is stable only when the evolution is close to unitary. 
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Here we used the fact that the average trace vanishes so the term (S(i)) is zero. Decomposition 
of a given quantum algorithm U into quantum gates is by no means unique. An interesting 
question then is, given perturbations Vt, which form of the algorithm has the highest fidelity, 
i.e. is the most stable? We will see that the standard QFT algorithm can be rewritten in a non 
trivial way so that it becomes more stable against static perturbations. The guiding principle 
in the construction of this new algorithm will be to study the correlation function 1)6. 5(1 and 
trying to minimise its sum. 

6.1.1 Time Dependent Perturbations 

If our perturbation Vt is time dependent, i.e. we have different perturbations on different gates, 
then the decay of correlation function will not only depend on the unperturbed dynamics, but 
also on how strong these perturbations are correlated at different gates. Let us have a closer 
look at one extreme example. If the perturbation Vt is an uncorrelated noise, as it would be in 
the case of coupling to an ideal heath bath, then the matrix elements of Vt may be assumed to 
be Gaussian random variables which are uncorrelated in time, 

(V;jfc(i)V^m(t'))noise = J^^jmhAt'- (6.7) 

Hence one finds {C{t,t'))^^^^^ = 5tt', where we have averaged over noise. In fact the average 
of the product in Mg 1)6. 3|) equals to the product of the average and yields the noise- averaged 
fidelity 

(^m)noise=exp(-52r), (6.8) 

which is independent of the quantum algorithm U. This result is completely general provided 
that the correlation time of the perturbation is sma//er than the duration of a single gate. On the 
other hand, for a static perturbation Vt = V one may expect slower correlation decay, depending 
on the 'regularity' of the evolution operator U, and hence faster decay of fidelity. Importantly, 
note that in a physical situation, where perturbation is expected to be a combination Vt = 
V^tatic + ynoise{t), the fidelity drop due to a static component is expected to dominate long-time 
quantum computation T ^ oo (i.e. large number of qubits n) over the noise component, as soon 
as the quantum algorithm exhibits long time correlations of the operator Kitatic- If Kitatic = 
is zero, the quantum computation can be stabilised by making "adiabatically" slow evolution 
of gates. In the following we will focus exclusively on the static perturbations being the worst 
ones. 

6.2 Quantum Fourier Transformation 

We will consider Quantum Fourier Transformation algorithm (QFT) and will consider its stabil- 
ity against static random perturbations. The perturbation Vt ^V will be a random hermitian 
matrix from a Gaussian unitary ensemble (GUE). Gaussian unitary ensamble is invariant under 
unitary transforma tions and the matrix elements in an arbitrary basis are independent random 
Gaussian variables ( Mehtal . ll99ll l. Due to hermitian symmetry they are real on the diagonal and 



complex off-diagonal. The GUE matrices can be used to model quantum statistical properties 
of classically chaotic Hamiltonians and have been first applied to studies of nuclear resonances. 
Second moments of a GUE matrix V are normalised as 

{yjkVlm)QUE = ^jmSkl/J^, (6.9) 

where the averaging is done over a GUE ensemble. 
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Let us briefly describe the QFT algorithm. Basis qubit states in a Hilbert space of dimension 
AA = 2" will be denoted by \k), fc = 0, . . . , 2" — 1. The unitary matrix J/qft performs the 
following transformation on a state with expansion coefficients Xk 

J\f-l M-l 

U(^F^{Y,Xk\k))=Y.Xk\k), (6.10) 

fc=0 fc=0 

where Xk = —tt?'}2j=q ex.p{2iTijk/M)xj. The resulting expansion coefficients Xk are Fourier 
transformed input coefficients Xk- The "dynamics" of the QFT is decomposed into three kinds 
of unitary gates: One-qubit gates Aj acting on j-th qubit 

where the basis is ordered as (|0), |1)), diagonal two-qubit gates Bj^ = diagjl, 1, l,exp (i^jfc)}; 
with 6jk = 7r/2 ■'j and transposition gates Tjfc which interchange the j-th and k-th qubit, 
Tjk\ . . . j . . .k . . .) = \ . . .k . . . j . . .). There are n A-gates, n{n — l)/2 B-gates and [n/2] transpo- 
sition gates, where [x] is the integer part of x. The total number of gates for the algorithm is 
therefore T = [n{n + 2)/2]. For instance, in the case of n = 4 we have a sequence of T = 12 
gates (time runs from right to left) 

f^QFT = T03T12A0B01B02B03A1B12B13A2B23A3. (6.12) 

For the GUE perturbation we can average the correlation function C{t,t') (|6.6j) over the 
GUE ensemble, resulting in 



(C(t,t')> 



GUE 



■^trC/(t,t') 



(6.13) 



where U{t,t') is the unperturbed propagator from gate t' + 1 to t, U{t,t') := Uf • • C/f+i with 
the convention U{t,t) = 1. Averaging over the GUE is done only to ease analytical calculation 
and to yield a quantity that is independent of a particular realization of the perturbation. 
Qualitatively similar (numerical) results are obtained without the averaging. 

We have calculated the correlator {C{t, t')) for the QFT 1)6. 13() which is shown in the right 
Figure lol One can clearly see square red plateaus on the diagonal due to blocks of successive 
B-gates. Similar square plateaus can also be seen off diagonal (from orange, yellow to green), 
so the correlation function has a staircase-like structure, with the A-gates responsible for the 
drops and B-gates responsible for the flat regions in between. This can be easily understood. 
For "distant" qubits A: — j » 1 the gates Bj^ are close to the identity and therefore cannot 
reduce the correlator. This slow correlation decay results in the correlation sum 

T 

v:=Y. Cit,t'), (6.14) 

t,t'=i 

being proportional to v oc n"^ (sum of the first n squares) as compared to the theoretical 
minimum u ccT ccn^. 

In view of this, we will now try to rewrite the QFT with a goal to achieve a smaller correlation 
sum, hopefully v cc v? . From (|6.1.S)) we learn that the gates that are traceless (e.g. A-gates) 
reduce the correlator very efficiently. In the plain QFT algorithm ()6.12|) we have n — 1 blocks 
of B-gates, where in each block all B-gates act on the same first qubit, say j. In each such 
block, we propose to replace B^^ with a new gate Gj^ = RL,Bjfc, where a unitary gate Rjk will 
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Figure 6.1: Correlation functions (C'(i, t'))GUE ^l^-l'^P ^or n = 10 qubits and static GUE per- 
turbation. The right figure shows the standard QFT (|H.12|) with T = 60, while the left figure 
shows the IQFT with T = 105 gates. Colour represents the size of elements in a log-scale from 
red (e^") to blue (e~^^ and less). 



be chosen so as to commute with all diagonal gates Bj; in the block, whereas at the end of the 
block we will insert Rj^ in order to "annihilate" RL so as to preserve the evolution matrix of a 



.t 



whole block. The unitarity condition RLRj/j = 1 and [Rjk, ^ji] = for all j, k, I leaves us with a 
6 parametric set of matrices Rjfc- By further enforcing trRjfc = in order to maximally reduce 
the correlator, we end up with 4 free real parameters in Rjk- One of the simplest choices, that 
turned out to be as suitable as any other, is the following 



R 



'jk 





1 
Vo 





1 







1 






-l/ 



(6.15) 



with the basis states {\jk)) ordered as (00,10,01,11). The R gate can be compactly written 

.Uj... {jTj (B bk) ■ ■ ■), where © is an addition modulo 2, bar 



as R 



'jk\ 



bk...) = {-ir 



denotes a negation and aj, b^ are or 1. Furthermore, we find that R- gates also commute among 
themselves, [Rjfc,Rj7] = 0, which enables us to write a sequence of R-gates whichever way we 
like, e.g. in the same order as a sequence of G's, so that pairs of gates Gjk, Rjfc operating 
on the same pair of qubits {j,k), whose product is a bad gate Bj^, are never neighbouring. 
This is best illustrated by an example. For instance, the block B01B02B03 will be replaced by 
RoiR'02Ro3R'OiBoiRo2Bo2Ro3Bo3 = RoiR'02Ro3GoiGo2Go3. This is how we construct an improved 
quantum Fourier transform algorithm (IQFT). For the IQFT we need one additional type of 
gate, instead of diagonal B-gates, we use nondiagonal ones R and G. To illustrate the obvious 
general procedure we write out the whole IQFT algorithm for n = 4 qubits (compare with 

UiQFT = T03T12A0R01R02R03G01G02G03A1R12R13G12G13A2R23G23A3. (6.16) 

Such IQFT algorithm consists in total of T = [n(2n+ 1)/2] gates (note that it does not pay of to 
replace a block with a single B gate as we have done, so we could safely leave B23 = R23G23). The 
correlation function for the IQFT algorithm is shown on the left of Figure 16.11 Almost all off- 
diagonal correlations are greatly reduced (to the level oc 1/AA^), leaving us only with a dominant 



112 



Chapter 6. Application: Quantum Computation 




Figure 6.2: Dependence of the fidehty F{T) on the number of qubits n for the QFT (pluses) 
and the IQFT algorithms (crosses), for fixed 6 = 0.04. Numerical averaging over 50 GUE 
realizations is performed. The full curve is exp (— (5^{0.47n^ — 0.76n^ + 2.90n}) and the dashed 
one exp (— (5^{1.22n^ + 1.78n}). For n = 10 the trace is approximated by an average over 200 
Gaussian random register states. 



diagonal. If we would have only diagonal elements, the fidelity would be F{T) = 1 — 6'^T, (as 
in the case of noisy perturbation or decoherence, however, with a different physical meaning 
of the strength scale 6) where the number of gates scales as T oc n^. From Figure l6?T] it is 
clear that we have a very fast correlation decay for the IQFT. Studying the scaling with n, 
the correlation sum v has decreased from v oc n^ to v cc n'^. To further illustrate this, we 
have numerically calculated fidelity by simulating the quantum algorithm and applying the 
perturbation exp(—[6V) at each gate. The results are shown in Figure W^ where one can see 
much faster decay of the fidelity for the QFT than for the IQFT. Note that the IQFT has 
higher fidelity despite the perturbation being applied ~ n? times for the IQFT and only ~ nP/l 
times for the QFT. As we have argued before, the sum of 2-point correlator (jHISI) gives us only 
the first nontrivial order in the (5-expansion. For dynamical systems, being either integrable or 
mixing and ergo die, we have shown in previous chapters that also higher orders of the fidelity 
can approximately be written as simple powers of the correlation sum z/, so that the fidelity has 
the simple form 

F{T) ^ eiip{-u6^). (6.17) 

Although quantum algorithm is quite inhomogeneous in time, we may still hope that exp {—vd'^) 
is a reasonable approximation of the fidelity also to higher orders in 6. This is in fact the case 
as can be seen in Figure l6.l-{[ Note also that the leading coefficient in the exponent for the 
IQFT fidelity, lim„^oo J^/n^ = 1.22, is close to the theoretical minimum of 1. As the definition 
of what is a fundamental single gate is somehow arbitrary, the problem of minimising the sum 
I' depends on a given technical realization o f gates and the nature of the perturbation V for 



an experimental setup. It has been shown (jGelardo et al. 1 1200,11 ). that the IQFT algorithm 
improves stability against GUE perturbations also for a more realistic model of a quantum 
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Figure 6.3: Dependence of the fidelity F{T) on 6 for the QFT (pluses) and the IQFT (crosses), 
for fixed n = 8. Solid curves are functions exp(— 1/5^) (see text) with u calculated analytically 
(j?rn|) and equal to i/ = 216 for the QFT and 93.2 for the IQFT. 



computer, namely for an Ising quantum computer. 

We should mention that the optimisation becomes harder if we consider few-body (e.g. two- 
body) random perturbation. This is connected with the fact that quantum gates are two-body 
operators and can perform only a very limited set of rotations on a full Hilbert space and 
consequently have a limited capability of reducing correlation function in a single step. For 
such errors the fidelity will typically decay with the square of the number of errors T (i.e . 
gates), like ~ exp(— J^T^), that is the same as for regular systems, see also (|Gea-Banaclochel . 
|2000|) for a similar result. This means that the very fact that the algorithm is efficient, having a 
polynomial number of gates, makes it very hard to reduce the correlation function and therefore 
causes a fast fidelity decay. However, our simple approach based on n-body random matrices 
seems reasonable, if the errors due to unwanted few-body qubit interactions can be eliminated 
by other methods. 



Chapter 7 

Conclusion 



I hate quotations. Tell me what you know. 

— Ralph Waldo Emerson 



We have studied the decay of quantum fidehty, of reduced fidehty and of purity in quantum 
systems. We considered two extreme cases of system's dynamics, mixing and regular. The 
dependence on the initial conditions, in particular for a random and coherent initial state, and 
the influence of the perturbation type on quantum stability has been analysed. 

For a general type of perturbation, having a nonzero diagonal matrix elements in the eigenba- 
sis of the unperturbed system, i.e. having a nonzero time average, the fidelity decay depends on 
the mixing properties of system dynamics. For mixing systems in the Fermi golden rule regime 
fidelity decay is exponential, with the decay time given by the transport coefficient, which is in 
turn the integral of the correlation function of the perturbation. The decay time scales with 
Planck's constant and with the perturbations strength as Tm ~ ?i /5^. In this regime quantum 
fidelity decay is much slower than classical fidelity decay and moreover, it will in general be 
the slower the more chaotic the corresponding classical system is. This surprising result does 
not violate the quantum-classical correspondence though, as for large perturbation strengths, 
in the so called Lyapunov regime, quantum fidelity agrees with classical fidelity. Whether we 
observe the quantum or the classical behaviour depends on the order of two noncommuting 
limits, namely the semiclassical limit of vanishing Planck constant h ^ and the limit of van- 
ishing perturbation strength 5^0. For sufficiently long times fidelity decay in mixing systems 
does not depend on the initial state. In regular systems and for perturbations with a nonzero 
time average the fidelity decay for wave packets is governed by a ballistic separation of the 
packets. For coherent initial states the resulting decay is Gaussian with the decay time scaling 
as Tr ~ yh/6, and can be for sufficiently small perturbations smaller than for chaotic systems. 
For random initial states in a regular regime quantum fidelity decays according to a power law, 
F ~ {h/6t) in d degrees of freedom system. In regular systems we also considered the decay 
averaged over random positions of the initial coherent state, resulting in an algebraic decay but 
with the power being system specific. By a semiclassical method we theoretically calculated all 
decay times in both cases of regular and mixing dynamics in terms of classical quantities only, 
despite the fact that the quantum and the classical fidelity do not agree for mixing systems. 

The quantum fidelity decay is markedly different for perturbations with a zero time average, 
which can be written as a time derivative of another operator. For such perturbations fidelity 
freezes at a constant plateau regardless of the dynamics and starts to again decay only after 
a much longer time, scaling as ~ 1/(5. This freezing is a pure quantum phenomenon as the 
correspondence with classical fidelity ends before the plateau starts and classical fidelity decays 
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much faster. The only exception where the classical fidelity also exhibits freezing and agrees 
with quantum fidelity is for a harmonic oscillator, for which the plateau is higher than for other 
systems. We explicitly calculate the plateau value for mixing and regular dynamics and coherent 
and random initial states. For mixing dynamics a universal relation holds between the plateau 
-^plat f°^ random and -FS^ for coherent states, (F^t^)'^ = F^f^. In the linear response regime 
the scaling of the plateau is 1 — -Fpiat ~ 5"^ /h'^ regardless of the dynamics. The asymptotic decay 
of fidelity after the plateau ends is also theoretically calculated in terms of a "renormalised" 
perturbation with strength (5^/2. For regular systems this long time Gaussian decay happens 
on a time scale ~ yh/d'^ for coherent initial states whereas it is power law with the prefactor 
scaling as ~ h/5'^ for random initial states. For mixing dynamics the long time decay does not 
depend on the initial state and is exponential with the decay time ~ h /5 or Gaussian with the 
decay time ~ ^ ~ ' /(5^. The crossover from the exponential to the Gaussian decay happens at 
the Heisenberg time. For one dimensional regular systems we also explain echo resonances, a 
sudden revivals of fidelity. 

We also study composite systems, composed of a central system and an environment, and 
connect the decay of purity with the decay of reduced fidelity. We prove a rigorous inequal- 
ity between the fidelity (characterising the stability of the whole system), the reduced fidelity 
(characterising the stability of the central system) and the purity (characterising the entan- 
glement). For mixing systems fidelity, reduced fidelity and purity all decay on the same time 
scale. For regular systems though, reduced fidelity has a Gaussian decay for coherent states 
whereas purity decays on an fi independent time scale. We explicitly calculate the purity decay 
and the power of the asymptotic algebraic decay depends on the perturbation and can range 
between 1 and dc, the dc being number of degrees of freedom of the central system. All decay 
constants are explicitly calculated. We also discuss an interesting case where the time scales of 
the central system and the environment are vastly different and one can use averaging over the 
faster system to simplify the theory. Decoherence for macroscopic superpositions of coherent 
states is derived and shown to be faster than for a single coherent state. 

Finally, we show an application of fidelity theory. By "randomising" the quantum Fourier 
transform algorithm we are able to make it more resistant against random perturbations from 
the environment. 



Appendix A 

Spin Wigner functions 



The Wigner function enables us to represent the quantum density matrices in a phase space and 
thereby compare it with the classical probability densities. If we have a one particle quantum 
system, described by a canonical pair [q,p\ = \h, the Wigner function Wp of a quantum state 
given by a density matrix p is 

Wp{q,p) := — / dx{q - x\p\q + x) exp {-i2px/h), (A.l) 

if |x) is an eigenstate of operator q with an eigenvalue x. In this appendix we will use a hat for 
quantum operators. For spin state such a definition can not be used as the Hilbert space has 
different structure due to SU(2) commutation relations of spin operators. 

We would like to obtain a Weyl symbol of an arbitrary operator A acting on a Hilbert space 
of size 25 + 1, i.e. on a state space of spin of size S. In the special case when the operator 
A is equal to a density matrix, the re sulting \yeyl sy mbol is called a Wigner function. Such 
functions have been first proposed bv lAgarwall (|l98ll ). The Weyl symbol Wa(^;V^) will be a 



function of coordinates on a sphere. Furthermore, we would like the standard trace dot product 
for operators to carry over to Weyl symbols, i.e. we demand that the following equality should 
hold 

tr(i,#) = fwAWsdn, (A.2) 

with dO = dipsmiDM. Spherical harmonic functions Yf^ constitute an orthonormal basis on a 
sphere, 

''YrYf*dn = 5u'Smm', (A.3) 



and so we can expand the Weyl symbol over Yf^. On the other hand, an arbitrary operator A 
can be in turn expanded over multipole operators T™ forming an orthogonal basis in the Hilbert 
space of operators, 

tr{frff^)=Sl,6,nrr.,. (A.4) 

Therefore we can write the operator A as 

2S I 

^ •= E E ^irnfr, ai^ = tr(if;'"t). (A.5) 

l=Qm=-l 

From this and the orthogonality property of Yj^ we immediately see that the Weyl symbol 

defined as 

25" ; 

Wa{^,v):=Y^ E ai^Yn^,^), (A.6) 

1=0 m=-l 
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will satisfy dot property ()A.2|) we demanded. The coefficients aim can be calculated using the 
explicit form of the multipole operators 

T}'" := V (-l)'^~™-^V2rTT f ^ ^ ^)\m + q){q\, (A.7) 

q— — S 

( S I s\ 

with I I being the Wiener 3? symbol. 

\-m- q m q J ^ ^ •' -^ 

The spin Wigner function for a state represented by a density matrix p is obtained by taking 
A = p m. the above formulas. As the density matrix is a Hermitian operator, the resulting 
Wigner function is real. In the case of two pure states, pi = |'(/'i)(V'i|) P2 = |V'2)(V'2| the dot 
condition gives simply 

\{iJi\^2)? = I Wp,Wp,dn. (A.8) 



For examples of Wigner functions of some simple states see e.g. ( Dowling et al. Lll994^ . 



Appendix B 

Coherent State Expectation Values 



Boson coherent states 

Expectation values of expressions involving creation and annihilation operators a"^ and a for a 
harmonic oscillator coherent state |a) are frequently needed. Using the definition of |q) (|2..Sn|) 
it is easy to show the following equality 

{a\g{a+)f{a)\a) = g{a*)f{a), (B.l) 

with two polynomials g and /. Let p{a, a^) be some polynomial function and by : p{a, a^) : we 
will denote a polynomial in a normal order which can be obtained from p{a,a~^) by using the 
commutation relation [a, a+] = 1 to bring all terms involving a"*" to the left of terms with a. For 
instance, ■.aa~^ : = a^a+ 1. Then we can write the expectation value of an arbitrary polynomial 
as 

{a\p{a,a~^)\a) = {a\ ■.p{a,a'^) : \a) =: p{a,a*) : . (B.2) 

Spin coherent states 

For spin coherent states formulas are a bit more complicated due to different group struc- 
ture. The easiest systematic method for the calculation of expectation values of polynomi- 
als i n operators ^y^ in coherent state \'d*,ip*) (|2.2H|) is using generating function formal- 
ism ( Arecchi et al. V Il972h . For brevity let us denote a spin coherent state with a complex 



parameter r := e"'' tan(??*/2), i.e. \t) := \'&*,ip*). The following expression holds, 

where 

X(^ r, C) ■= (ri le^^+e^'^^e^^- Iro) = ^ — -— ^ g'^^i"^^)^ fB 4) 

This formula, together with the commutation relations, can be used to calculate some of the 
lowest powers 



{r\S.\r) -- 


= Sz 


{r\S'jr) - 


- 5V + |(i-.2) 


{t\S!\t) - 


= 5V + ^.(l-z2)-|.(l-.2) 


{r\S^\r) -- 


4 
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14z2 + 3) + |(-3z^ + 4z2_l), (B.5) 
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where z = cos 'd* . The expressions for other spin operators Sy and Sx are obtained by replacing 
z with y or X on the right hand side, respectively. Some other useful expectation values are 

{t\S+S^\t) = 52(l-z2) + |(i + ^)2 

{t\S^S+\t) = S\l-z') + ^{l-zf. (B.6) 
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